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^ ■ Abstract. We present a method to obtain congruences modulo powers of 2 for se- 

quences given by recurrences of finite depth with polynomial coefficients. We apply 
this method to Catalan numbers, Fuf5-Catalan numbers, and to subgroup counting 
functions associated with Hecke groups and their lifts. This leads to numerous new 
results, including many extensions of known results to higher powers of 2. 
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1. Introduction 



Ever since the work of Sylow [35j, Frobenius [HI |T3], and P. Hall [T7], the study 
of congruences for subgroup numbers and related numerical quantities of groups has 
(N \ played an important role in group theory. 

Divisibility properties of subgroup numbers of (finitely generated) infinite groups 

may to some extent be viewed as some kind of analogue to these classical results for 

^ \ finite groups. To the best of our knowledge, the first significant result in this direction 

'' was obtained by Stothers [M]: the number of index-n- subgroups in the inhomogeneous 

Q ■ modular group PSL2{Z) is odd if, and only if n is of the form 2^ — 3 or 2^^^ — 6, for 

some positive integer k > 2. A different proof of this result was given by Godsil, Imrich, 

and Razen [T3]. 

The systematic study of divisibility properties of subgroup counting functions for 
S^ • infinite groups begins with [27j. There, the parity of subgroup numbers and the number 

of free subgroups of given finite index are determined for arbitrary Hecke groups Sj{q) = 
C2*Cq with g > 3. Subsequently, the results of [27] were generalised to larger classes of 
groups and arbitrary prime modulus in [3l |20l |25l |26l [28]. A first attempt at obtaining 
congruences modulo higher prime powers was made in [29], where the behaviour of 
subgroup numbers in PSL2{Z) = S^{3) is investigated modulo 8 and a congruence 
modulo 16 is derived for the number of free subgroups of given index in PSL2{'L). 

A common feature of all the above listed sequences of subgroup numbers is that 
they obey recurrences of finite depth with polynomial coefficients. The purpose of this 
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paper is to present a new method for determining congruences modulo arbitrarily large 
powers of 2 for sequences described by such recurrences. Our method is inspired by the 
observation that many of the aforementioned results say in essence that the generating 
function for the subgroup numbers under consideration, when reduced modulo a 2- 
power, can be expressed as a polynomial in the basic series 

$w = E^" (1.1) 

n>0 

with coefficients that are Laurent polynomials in z. What our method affords is an 
algorithmic procedure to find such polynomial expressions, provided they exist. By 
applying our method to Catalan numbers, to (certain) FuB-Catalan numbers, and to 
various subgroup counting problems in Hecke groups and their lifts, we obtain far- 
reaching generalisations and extensions of the previously mentioned results. In order 
to give a concrete illustration, the recent result [221 Theorems 6.1-6.6] of Liu and 
Yeh determining the behaviour of Catalan numbers Cat„ modulo 64 can be compactly 
written in the form 

oo 

^ Cat„ z" = 32^^ + 16^^ + 6^2 + 13^ + 1 + (32z^ + 32^^ ^ 20^^ + 44z + 40) $(z) 

n=0 

+ [ 16^^ + 56;^^ + 30z + 52 + — J $^(z) + (32^^ + 60^ + 60 + — j $^(^) 
+ i'i'lz^ + 16^2 + 48z + 18 + — J ^^{z) + (32^2 + 44) ^^{z) 

+ (48^ + 8 + — J $^(^)+ (32^ + 32 + - j <l>^(z) modulo 64, (1.2) 

as may be seen by a straightforward (but rather tedious) computation. Our method 
can not only find this result, but it produces as well corresponding formulae modulo 
any given power of 2 in a completely automatic fashion, see Theorems [13] and [TH in 
Section O 

In a sense, which is made precise in Section HJ our method is very much in the spirit of 
Doron Zeilberger's philosophy that mathematicians should train computers to automat- 
ically produce theorems. Indeed, Theorems [T3| [T9l [211 [33|[36[ imply that our algorithm is 
able to produce a theorem on the behaviour modulo any given 2-power of the subgroup 
counting functions featuring in these theorems, and, if fed with a concrete 2-power, 
our implementation will diligently output the corresponding result (provided the input 
does not cause the available computer resources to be exceeded . . . ). Moreover, when 
discussing subgroup numbers of lifts of PSL2{1), (such as the homogeneous modular 
group S'L2(Z)), a crucial role is also played by an application of the holonomic func- 
tions approach to finding recurrences for multi-variate hypergeometric sums, pioneered 
by Wilf and Zeilberger ^371 [39] , and further developed in [U [5l [19] . 

The rest of this introduction is devoted to a more detailed description of the contents 
of this paper. In Section [2] we discuss our main character, the formal power series ^{z) 
defined in (ll.ip . While ^{z) is transcendental over Q[2;] (or, equivalently, over Z[z]), 
it is easy to see that it is algebraic modulo powers of 2. The focus in that section is 
on polynomial identities for ^{z) modulo a given 2-power which are of minimal degree. 
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Then, in Section HI we describe our method of expressing the generating function of 
a recursive sequence, when reduced modulo a given 2-power, as a polynomial in $(2) 
with coefficients that are Laurent polynomials in z. The method relies in an essential 
way on the polynomial identities from Section |2l The problem how to extract the 
explicit value of a concrete coefficient in a polynomial expression in $(2;) (such as (11. 2p ) 
modulo a given 2-power is discussed in Section [3l where we present an efficient algorithm 
performing this task. This algorithm is of theoretical value (minimal length relations 
between powers of $(-2) such as the ones in Proposition |2] are established by applying 
this algorithm to the powers of $(z); see also Appendix \^ as well as of practical 
significance, as is demonstrated by the derivations of Theorems [27] and EH 

As a first illustration of our method, we apply it to Catalan numbers, thereby sig- 
nificantly improving numerous earlier results in the literature; see Section [51 This is 
contrasted in Section [6] with an example (concerning particular Fufi-Catalan numbers) 
where our method is bound to fail. The reason is spelled out in Theorem [151 which, 
at the same time, also gives a new description for the parity pattern of the numbers of 
free subgroups of given index in the Hecke group i3(7). 

The subsequent sections contain several applications of our method to the problem 
of determining congruences modulo a given 2-power for numbers of subgroups of Hecke 
groups ^[fi) and their lifts 

r„(g) = C2™ * C,„ = (x,y I x2- = y-?- = 1, x2 = y-?), m > 1. (1.3) 

Ubiquitous in these applications is — explicitly or implicitly — the intimate relation 
between subgroup numbers of a group F and numbers of permutation representations 
of r, in the form of identities between the corresponding generating functions. This is 
directly visible in the folklore result (lll.ip (which is not only used in Sections [Til ESI 
and [HI but also lies behind the crucial differential equation (19. ip in Section [9l cf. its 
derivation in [H]), and also indirectly in Lemma [T71 via the A-invariants yl^(io(g)), see 
\2n\ Sec. 2.2]. In the cases relevant here, the numbers of permutation representations of 
r satisfy linear recurrences with polynomial coefficients — to make these explicit may 
require the algorithmic machinery around the "holonomic paradigm" (cf. [H El [13 [231 
[311 [371 [39] ) , see Sections [TTI and [T2l for corresponding examples. Via the aforementioned 
generating function relation, such a recurrence can be translated into a Riccati-type 
differential equation for the generating function of the subgroup numbers that we are 
interested in. It is here, where our method comes in: it is tailor-made for being applied 
to formal power series F(^z) satisfying this type of differential equation, and it affords 
an algorithmic procedure to find a polynomial in $(2;) which agrees, after reduction of 
the coefficients of F(^z) modulo a given power of 2, with the power series Fi^z). 

We start in Section [7] with free subgroup numbers of lifts Vmisi)-, for primes g > 3, 
where we prepare the ground for application of our method. More specifically, in 
Proposition [18] we present a lower bound for the 2-adic valuation of the number of 
free subgroups of index n in Yra{(l)i where g is a Fermat prime. In particular, this result 
implies that the sequence of free subgroup numbers under consideration is essentially 
zero modulo a given 2-power in the case when m is even. In Section [8] we show that our 
method provides an algorithm for determining these numbers of free subgroups of Fm(3) 
modulo any given 2-power in the case when m is odd. The corresponding results (see 
Theorems [T^ and [2U|) go far beyond the previous result j29[ Theorem 1] on the behaviour 
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of the number of free subgroups of PSL2 (Z) modulo 16. Our method provides as well an 
algorithm for determining the number of all subgroups of index n in PSL2{'L) modulo 
powers of 2, as we demonstrate in Section [9l Not only are we able to provide a new 
proof of Stothers' result [31] (which was stated in the second paragraph above), but our 
method leads as well to refinements modulo arbitrary powers of 2 of Stothers' result 
and of the mod-8 results in [29l Theorem 2] mentioned earlier; see Theorems [21] and 
[2H For the homogeneous modular group S'L2(Z) (being isomorphic to the lift r2(3)) 
and for the lift r3(3), however, our method from Section H] fails already for the modulus 
8 = 2^. We overcome this obstacle by instead tuning our computations with the target 
of obtaining results modulo 16 = 2^. Indeed, this leads to the determination of the 
number of subgroups of index n in SL2{'L) and in r3(3) modulo 8 (see Theorems [2^ 
and [30]), but direct apphcation of our method does not produce corresponding results 
modulo 16. Only by an enhancement of the method, which we outline in Appendix [D] 
we are able to produce descriptions of the subgroup numbers of SL2{1) modulo 16, see 
Theorem [28] For the subgroup numbers of Fa (3) even this enhancement fails, and this 
shows that the generating function for these subgroup numbers, when coefficients are 
reduced modulo 16, cannot be represented as a polynomial in ^{z) with coefficients 
that are Laurent polynomials in z. Still, the results in Theorems [22] [2H1 and [30] go 
significantly beyond the earlier parity results [2D1 Eq. (6.3) with \H\ = 1] for these 
groups. This is explained in Sections [H] and [T^J with Section [TD] preparing the ground 
by providing formulae for the number of permutation representations of 6*^2 (Z) as well 
as other lifts of PSL2{1^). A further example where our method works for any 2-power 
is the subject of Section [13] there we apply the method to a functional equation (see 
(113. ip ) extending the functional equation for Catalan numbers (producing Fufi-Catalan 
numbers), and show that it works for any given 2-power; see Theorem [33] In fact, we 
apply a variation of the method here, in that the basic series ^{z) gets replaced by 
a slightly different series, which we denote by ^h{z) (see (113. 3p ). If Theorem [33 is 
combined with results from [27j, then it turns out that our method provides as well an 
algorithm for determining the number of free subgroups of index n in a Hecke group S) (q) 
and its lifts, where g is a Fermat prime, modulo any given 2-power; see Corollary [M] The 
same assertion holds as well for the problem of determining the number of subgroups 
of index n in the Hecke group S^{5), again modulo any given 2-power (see Theorems [36] 
and [37] in Section [H]). We conjecture that the same is true for Hecke groups S){q), with 
q a Fermat prime (see Conjecture [38]). The results of Sections [T3] and [T^ discussed above 
largely generalise the parity results ^7\ Cor. A', respectively Theorem B] for subgroup 
numbers of Hecke groups, although our results are not independent, in the sense that 
we base our analyses on prior results from [27j. 

Concluding the introduction, we remark that there is no principal obstacle to gen- 
eralising our method to other basic series and moduli. For example, one may think of 
analysing the behaviour of recursive sequences modulo powers of any prime p in terms 
of the obvious generalisation of ^{z), i.e., the series J2n>o^^"- ^^ ^^ i^ ^^'^^ ^^^ difficult 
to see that our results from Sections [5]and[T3]for FuB-Catalan numbers characterised by 
the functional equation (113. ip for their generating function have rather straightforward 
analogues for FuB-Catalan numbers whose generating function satisfies the functional 
equation 

zfP\z)-f{z) + 1 = 0. (1.4) 
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However we are not aware of any applications of this (or of variants) to congruence 
properties of subgroup numbers modulo powers of primes p different from 2. In fact, 
the known results (cf. [22] or f2^ Theorem 3]) strongly point to the fact that the 
phenomena that appear modulo primes different from 2 cannot be captured by series 
of the type ^„>o z^" ■ So, currently, we do not know of interesting applications in this 
direction, but we hope to be able to return to this circle of ideas in future pubhcations. 

Note. This paper is accompanied by several Mathematica files and a Mathematica 
notebook so that an interested reader is able to redo (most of) the computations that 
are presented in this article. Files and notebook are available at the article's website 



http : //www . mat . univie . ac . at/~kratt/artikel/modlif ts . html 



2. The 2-power series $(z) 

Here we consider the formal power series ^{z) defined in 01. ip . This series is the 
principal character in the method for determining congruences of recursive sequences 
modulo 2-powers, which we describe in Section |H It is well known that this series 
is transcendental over Zlz] (this follows for instance from the density argument used 
in the proof of Lemma [1] below) . However, if the coefficients of ^{z) are considered 
modulo a 2-power 2"^, then ^{z) obeys a polynomial relation with coefficients that are 
polynomials in z. The focus of this section is on what may be said concerning such 
polynomial relations, and, in particular, about those of minimal length. 

Here and in the sequel, given power series (or Laurent series) f{z) and g{z), we write 

f{z) = g{z) modulo 2"^ 

to mean that the coefficients of z^ in f{z) and g{z) agree modulo 2'^ for all i. 

We say that a polynomial A{z, t) in z and t is minimal for the m,odulus 2'*', if it is monic 
(as a polynomial in t), has integral coefficients, satisfies A{z., ^{z)) = modulo 2"*, and 
there is no monic polynomial B{z, t) with integral coefficients of t-degree less than that 
of A{z,t) with 5(2;, $(2;)) = modulo 2'^. (Minimal polynomials are not unique; see 
Remark O) Furthermore, we let ^2(0) denote the 2-adic valuation of the integer a, that 
is, the maximal exponent e such that 2^ divides a. 

The lemma below provides a lower bound for the degree of a polynomial that is 
minimal for the modulus 2"' . 

Lemma 1. If A{z,t) is m,inim,al for the modulus 2"' , then the degree d of A{z,t) in t 
satisfies f2((i!) > 7. In particular, the series $(-2) is transcendental over TL\z\. 

Proof. We introduce the following density function with respect to a given modulus 2'^ 
for a power series f{z) in z: 

D(/,2^;n) := |{e: 2"-i <e<2" and {z^) f{z) ^ modn\o 2^ }\, n = l,2,..., 

(2.1) 
where {z'^) f{z) denotes the coefficient of z'^ in f{z). Setting 

EUz):= Y. z2"^+2"^+-+2"'", (2.2) 

simple counting yields that 

D{Em, 2^; n) = J ~ -TTn""-', as n ^ 00. (2.3) 

\m — 1/ im — 1)1 
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(The modulus T^ does not play a role here.) Furthermore, by considering the binary 
representations of possible exponents e such that the coefficient of z^ in $"^(2;) does not 
vanish, we have 

D($'^, V; n) = {n'^-^) , as n ^ 00. (2.4) 

Indeed, by direct expansion, we see that 

^''\z)=m\E^{z) + R^{z), (2.5) 

where all monomials z'^ which appear with non- vanishing coefficient in Rm{z) have a 
binary expansion with at most m — 1 digits 1. Consequently, again by elementary 
counting, we have 

and hence D{Rm, 2'^; n) = 0(ra™~^) as n — )■ 00. 

Let us, by way of contradiction, suppose that the degree d of A(z, t) satisfies f 2(0?!) < 
7. Considering (12. 5 p with m = d, we see that Ed{z) appears in ^'^{z) with a non-zero 
coefficient modulo 2^ . Furthermore, the other terms in (^'^{z) (denoted by Rd{z) in 
(12. 5p ) have a density function modulo 2^ which is asymptotically strictly smaller than 
the density function of Ed{z). Consequently, if we remember (12. 3p . we have 

D($'^, 2^; n) ~ , , ' , n"'"^ = dn'^'\ as n ->cx). 

[d-iy. 

Moreover, by (12. 4p . all powers ^"^{z) with m < d have a density function modulo 2'^ 
which is asymptotically strictly smaller than n'^~^. Altogether, it is impossible that a 
linear combination of powers $™'((i), m = 0,1,. . . ,d, with coefficients that are poly- 
nomials in z sums up to zero modulo 2"', a contradiction to our assumption that d is 
the degree of a minimal polynomial A{z,t). The particular statement is an immediate 
consequence of the inequality just proven. D 

Proposition 2. Minimal polynomials for the moduli 2, 4, 8, 16, 32, 64, 128 are 
t^ + t + z modulo 2, 

if + t + zf modulo 4, 

t^ + 6t^ + {2z + 3)t2 + {2z + 6)t + 2z + bz^ modulo 8, 

{t^ + t + z) {t^ + 6t^ + {2z + 3)t^ + {2z + 6)t + 2z + 5^^) modulo 16, 

(t^ + 6t^ + (2z + 3)t^ + (2^ + 6)t + 2z + 5^2)2 modulo 32, 

(t^ + 6t^ + (2z + 3)t^ + (2^ + 6)t + 2z + 5^2)2 modulo 64, 

t* + 124t^ + t^(682 + 18) + t^(124z + 24) + t^ (622^ ^ 54^ ^ g^j 
+ t^ (20^2 + 76z + 28) + t^ {\\<oz^ + 114^2 + V2z + 92) 
+ t (I16;z^ + 28^2 + 8;z + 16) + 9z^ + 124z^ + 12^^ + 112z modulo 128. 

Proof. In order to be consistent with Section [31 let us change notation and write 

2"i+2"2-| 1-2"'"- 



i^M,...,iW := E 



z 

ni>--->nm>0 
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(with m occurrences of 1 in ifi i ... 1(2;)). Note that the above series is identical with 
the series which we earher denoted by Em{z). Straightforward calculations yield that 

<l>\z) = <l>{z)+2Hi^i{z)-z, (2.6) 

$3(2) = -2 ^ ^3.2" ^ 3^^ _ ^)^(^) ^ 6i/i,i(^) + 6i/i,i,i(^) - 3z, (2.7) 

n>0 

$4(2) = -12^^3.2'^ -8 Yl z'-^"'+'"'-8 J2 ^'"'^'■'"' + (13-18^)$(^) 

n>0 ni>n2>0 ni>n2>0 

+ (30 - l2z)Hi^i{z) + 36ifi,i,i(2) + 24i7i,i,i,i(2) + 5^^ - 13^. (2.8) 

In particular, relation (12.61) . together with Lemma [1], immediately implies the claims 
about minimal polynomials for the moduli 2 and 4. Moreover, a simple computation 
using (I2S])-([2S]) shows that 

$4(2) + 6$3(z) + (2z + 3)$2(z) + (22 + 6)$(z)+2z + 522 = modulo 8. (2.9) 

Together with Lemma [H this establishes the claims about minimal polynomials for the 
moduli 8, 16, 32, and 64. In order to prove the claim for the modulus 128, one uses the 
expressions for $^(2), i = 2, 3, . . . , 8, given above and in Appendix [Aj D 

Remark 3. Minimal polynomials are highly non- unique: for example, the polynomial 

it^ + t + zf + 2{f + t + z) 
is obviously also a minimal polynomial for the modulus 4. 

Based on the observations in Proposition [2] and Lemma [H we propose the following 
conjecture. 

Conjecture 4. The degree of a minimal polynomial for the modulus 2^ , 7 > 1, is the 
least d such that f2((i!) > 7. 

Remark 5. (1) Given the binary expansion of d, say 

d = do + di-2 + d2-4: + --- + dr-2'', < 7^ < 1, 
by the well-known formula of Legendre [2J4 p. 10], we have 



00 , 00 

£=1 L J £^]_ 



5^f/.2^" 



j=0 



EE*2- 



= 1 i=C 



= EE^*2^~' = 5^rf. {T - 1) =d-s{d), (2.10) 

i=l 1=1 i=l 

where s{d) denotes the sum of digits of d in its binary expansion. Consequently, an 
equivalent way of phrasing Conjecture H] is to say that the degree of a minimal polyno- 
mial for the modulus 2'^ is the least d with d — s{d) > 7. 

(2) We claim that, in order to establish ConjecturelU it suffices to prove the conjecture 
for 7 = 2*^ — 1, 5 = 1, 2, ... . If we take into account Lemma [1] plus the above remark, 
this means that it is sufficient to prove that, for each 6 > 1, there is a polynomial 
As{z,t) of degree 2^ such that 

As{z,^{z)) = modulo 22'-!. (2.11) 
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For, arguing by induction, let us suppose that we have aheady constructed Ai{z,t), 
A2{z,t), . . . , Am{z, t) satisfying (12.111) . Let 

a = ai •2 + a2-4 + --- + a^-2'", < a^ < 1, 

be the binary expansion of the even positive integer a. In this situation, we have 

rra mm 

YlA'^'{z,<l>{z)) = modulo ]^2"^(2*-i) _ "Qg^^Ca^-i) = 2"-^(°). (2.12) 

5=1 <5=1 (5=0 

On the other hand, the degree of the left-hand side of (12.121) as a polynomial in $(-2) is 

E^i«<52^ = a- ^_^ 

Let us put these observations together. In view of (I2.10p . Lemma [T] says that the 
degree of a minimal polynomial for the modulus 2"^ cannot be smaller than the least 
integer, S'^^ say, for which d^'^^—s{S'^^) > 7. (We remark that S'^^ must be automatically 
even.) If we take into account that the quantity a — s{a), as a function in a, is weakly 
monotone increasing in a, then (12.121) tells us that, as long as S"'^ <2 + 4 + --- + 2'" = 
2^+1 _ 2^ -^g have found a monic polynomial of degree d^'^\ B^(z,t) say, for which 
B.y(z, ^(z)) =0 modulo 2'^, namely the left-hand side of (I2.12p with a replaced by d^"'\ 
to wit 

m 

B,{z,t) = l[Af iz,t), 
5=1 

where d'-'^^ = d^ ■ 2 + dg ■ 4 + ■ ■ ■ + dm ■ 2"^ is the binary expansion of d^"'\ Hence, it 
must necessarily be a minimal polynomial for the modulus 2"'. 

Since (2™+^— 2) — s(2™"*"^— 2) = 2™+^— 2— m, we have thus found minimal polynomials 
for all moduli 2'^ with 7 < 2™+^ — -m — 2. Now we should note that the quantity a — s{a) 
makes a jump from 2*"+^ — m — 2 to 2™+^ — 1 when we move from a = 2™+^ — 2 to 
a = 2™+^ (the reader should recall that it suffices to consider even a). If we take 
A'^[z,t), which has degree 2 ■ 2"^ = 2'""''^, then, by (12. lip , we also have a minimal 
polynomial for the modulus (2^™"^) = 2^™ ~^ and, in view of the preceding remark, 
as well for all moduli 2'^ with 7 between 2™+^ — m — 1 and 2"^"*"^ — 2. 

So, indeed, the ffist modulus for which we do not have a minimal polynomial is the 
modulus 2^™ ~^. This is the role which Am+i{z,t) (see (12. lip with m + 1 in place of 
S) would have to play. 

The arguments above show at the same time that, supposing that we have already 
constructed Ai{z, t), A2{z, t), . . . , Am{z, t), the polynomial A'^{z, t) is a very close "ap- 
proximation" to the polynomial Am+i{z,t) that we are actually looking for next, which 
is only "off" by a factor of 2. In practice, one can recursively compute polynomials 
As{z,t) satisfying (12. lip by following the procedure outlined in the next-to-last para- 
graph before Lemma [6] in the next section. It is these computations (part of which are 
reported in Proposition |2]) which have led us to believe in the truth of Conjecture HI 

3. Coefficient extraction from powers of $(2;) 

In the next section we are going to describe a method for expressing formal power 
series satisfying certain differential equations, after the coefficients of the series have 
been reduced modulo 2^, as polynomials in the 2-power series ^{z) (which has been 
discussed in the previous section; for the definition see (II. ip ). the coefficients being 
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Laurent polynomials in z. Such a method would be without value if we could not, at 
the same time, provide a procedure for extracting coefficients from powers of $(2;). The 
description of such a procedure is the topic of this section. 

Clearly, a brute force expansion of a power $'^(z), where i^ is a given positive integer, 
yields 

r=l ai,...,ar>l 
ctlH \-a.i=K 

where 

ni>n2>--->nr>0 

The expansion (13. ip is not (yet) suited for our purpose, since, when ai, a2, . . . , a^ vary 
over all possible choices such that their sum is K, the series Ha^^a2,...,ar{^) are not linearly 
independent over the ring Z[z, z~^] of Laurent polynomials in z over the integer^, and, 
second, coefficient extraction from a series -ffai,a2,...,ar(-2^) can be a hairy task if some of 
the Oj's are even. 

However, we shall show (see Corollary [T]) that, if we restrict to odd a^'s, then the 
corresponding series -ffai,a2,...,ar(-2^); together with the (trivial) series 1, are linearly in- 
dependent over Z[z, 2;"^], and there is an efficient algorithm to express all other series 
Hbx,b2,...,bs{z)i where we do not make any restriction on the 6j's, as a linear combination 
over Z[2;, 2;"^] of 1 and the former series (see Lemma [9]). Since coefficient extraction 
from a series -ffai,a2,.--,ar(^) with all a^'s odd is straightforward (see Remark [8]), this 
solves the problem of coefficient extraction from powers of $(2;). 

As a side result, the procedure which we described in the previous paragraph, and 
which will be substantiated below, provides all the means for determining minimal 
polynomials in the sense of Section [2) as explained in Item (2) of Remark [5] at the end of 
that section, it suffices to find a minimal polynomial for the modulus 2^ ~^, 5 = 1, 2, . . . . 
For doing this, we would take a minimal polynomial As-i{z, t) for the modulus 2^ ~^, 
expand the square A^^_^{z,t), and replace each coefficient Ca^p of a monomial z^t^ in 
^5_i(2;, t) by Ca^/3 + 2^ ""^Xa^p, where x^^p is a variable, thereby obtaining a modified 
polynomial, Bs-i{z,t) say. Now we would substitute $(2;) for t, so that we obtain 
Bs-i{z, ^{z)). Here, we express powers of $(2;) in terms of the series -ffai,a2,...,a^(-2) with 
all Oj's being odd, and collect terms. By reading the coefficients of z^Ha^^a2,...,ar{z) in 
this expansion of Bs-i{z,^{z)) and equating them to zero modulo 2^ ~^, we produce 
a system of linear equations modulo 2^ ~^ in the unknowns Xa,i3- By the definition of 
As-i{z,t), after division by 2^ ~^, this system reduces to a system modulo 2, that is, 
to a linear system of equations over the field with two elements. A priori, this system 
need not have a solution, but experience seems to indicate that it always does; see 
Conjecture IH 

We start with an auxiliary result pertaining to the uniqueness of representations of 
integers as sums of powers of 2 with multiplicities, tailor-made for application to the 
series Ha^^a2,...,ar{z). 

The same is true for an arbitrary ring in place of the ring Z of integers. 
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Lemma 6. Let d,r,s be positive integers with r > s, c an integer with \c\ < d, and let 
ai, 02, . . . , flr respectively 6i, 62, ■ ■ • , ^s be two sequences of odd integers, with 1 < ai < d 
for 1 < i < r, and 1 < hi < d for 1 < i < s. If 

a{l^^'^ + a222("-^)'^ + ■ ■ ■ + ari^'^ = 6i2"^ + 622"^ + ■ ■ ■ + 6^2"^ + c (3.2) 

for integers ni,n2, . . . ,ns with ui > n2 > ■ ■ ■ > Us > 0, then r = s, c = 0, a, = hi, and 
Hi = 2d{r + 1 — 2) for i = 1,2, . . . ,r. 

Proof. We use induction on r. 

First, let r = 1. Then s = 1 as well, and (13. 2p becomes 

ai2^^ = 6i2'^i + c. (3.3) 

If ni > 2d, then the above equation, together with the assumption that ai is odd, 
implies 

22^ = c modulo 2^^^+^ 

However, by assumption, we have \c\ < d < 2^"^, which is absurd. 

If d < ni < 2d, then it follows from (13.31) that c must be divisible by 2"^ Again 
by assumption, we have \c\ < d < 2'^ < 2"^, so that c = 0. But then (13.31) cannot be 
satisfied since 61 is assumed to be odd. 

li < ni < d, then we estimate 

6i2"i +c<d{2'^ + l) <{2'^- 1)(2'^ + 1) < 2^"^, 

which is again a contradiction to (13.31) . 

The only remaining possibility is ni = 2d. If this is substituted in (13. 3p and the 
resulting equation is combined with \c\ < d < 2'^'^, then the conclusion is that the 
equation can only be satisfied if c = and oi =61, in accordance with the assertion of 
the lemma. 

We now perform the induction step. We assume that the assertion of the lemma is 
established for all r < R, and we want to show that this implies its validity for r = R. 
Let t be maximal such that nt > 2d. Then reduction of (13. 2p modulo 2^*^ yields 

6i+i2"*+^ + 6i+22"*+^ + ■ ■ ■ + &s2"^ + c = modulo 2^*^. (3.4) 

Let us write b ■ 2^"^ for the left-hand side in (13. 4p . Then, by dividing (13. 2p (with R 
instead of r) by 2^"^, we obtain 

^^22(«-i)'i + a22^(^-^^'^ + - ■ ■ + aij_i22'^ = 6i2"i-2'^+622''^-^'^ + - ■ ■+bt2'''-^'^+b-aR. (3.5) 
We have 

< 6 < 2-2'^rf (2^^-^ + 22'^-2 + ■ ■ ■ + 22^-"+* + 1) 

< 2-^'^d (2^^ - 2^'^-"+* + 1) <d. 

Consequently, we also have \b — aji\ < d. This means that we are in a position to apply 
the induction hypothesis to (13. 5p . The conclusion is that t = R — l,b — aR = 0,ai = bi, 
and Hi = 2d{R + 1 — z) for i = 1,2, . . . ,R— 1. If this is used in (13.21) with r = R, then 
we obtain 

aR2'"' = c 
or 

an2'"' = 6ij2"« + c. 
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depending on whether s = R — 1 or s = R. The first case is absurd since c < d < 2^"^ < 
0^2^'^. On the other hand, the second case has aheady been considered in (13.31) . and 
we have seen there that it follows that c = 0, a^j = fo^j, and n^ = 2d. 

This completes the proof of the lemma. D 

The announced independence of the series -ffai,a2,...,ar(-^) with all Oj's odd is now an 
easy consequence. 

Corollary 7. The series Ha^^a2,...,ar{z), with all ai 's odd, together with the series 1 are 
linearly independent over (Z/2Z)[z, z~^], and consequently as well over {Z/2'^Z)[z, z~^] 
for an arbitrary positive integer 7, and over %\z^ ^~^\- 

Proof. Let us suppose that 

N 

Po{z) + X]Pi(^)^aW,40^,,,^4«)(^) = 0, (3.6) 

where the Pi{z)'s are non-zero Laurent polynomials in z over Z/2Z (respectively over 
Z/2^Z or over Z), the Tj's are positive integers, and a^ , j = 1, 2, . . . , Tj, i = 1, 2, . . . , A^, 

are odd integers. We may also assume that the tuples (a^ , 03 , . . . , ar]), i = 1,2, . . . , N, 
are pairwise distinct. Choose iq such that rj^ is maximal among the Tj's. Without loss 
of generality, we may assume that the coefficient of z^ in pi^ [z) is non-zero (otherwise we 
could multiply both sides of (13. 6p by an appropriate power of z). Let d be the maximum 
of all Oj s and the absolute values of exponents of z appearing in monomials with non- 
zero coefficient in the Laurent polynomials Pi{z), i = 0, 1, . . . , A^. Then, according to 
Lemma |6] with r = Vi^, aj = clj , j = 1, 2, . . . , Tj^, the coefficient of 

(»0)22rd_^„(»0)22(r-l)d_^..._^_4>o)22d 



is 1 in if (ig) (io) iio){z), while it is zero in series z^H (i^) a^) (io){z), where e is a 
non-zero integer with lei < d, and in all other series z'^H (o (o (i)(z), i = 1, . . . ,io — 

'=' ' ' — ' a) 'aX ,...,ar- ^ '' 7 7 u 

l,io + 1,...,A^, where e is a (not necessarily non-zero) integer with |e| < d. This 
contradiction to (13. 6p establishes the assertion of the corollary. D 

Remark 8. Coefficient extraction from a series ifai,a2,...,a^(-2) with all a^'s odd is straight- 
forward: if we want to know whether z^^ appears in -ffai,a2v,ar(^)i that is, whether we 
can represent M as 

M = ai2"i + 022"" + ■ ■ ■ + a^2"'- 
for some ni,n2, . . . ,nr with ni > n2 > ■ ■ ■ > n,, > 0, then necessarily n,. = V2{M), 
Ur-i = V2{M — ar2""'), ctc. The term 2;*^ appears in ifai,a2,.--,ar(-2) if; ^^^ only if, the 
above process terminates after exactly r steps. This means, that, with n^.,n^_i, . . . ,ni 
constructed as above, we have 

M - (a,2"^ + ■ ■ ■ + a^_i2"'-^ + a^2"'-) > 

for s > 1, and 

M - (ai2"i + ■ ■ ■ + a,_i2"'-i + a,2"') = 0. 

It should be noted that, given ai, 02, . . . , a^, this procedure of coefficient extraction 
needs at most O(logM) operations, that is, its computational complexity is linear. 



12 M. KAUERS, C. KRATTENTHALER, AND T. W. MULLER 

Our next goal is to show that a series i^fei,62,.--,f's(-2^) can be expressed as a hnear 
combination over Z[z, z~^] of the series 1 and the series -ffai,a2,...,a^(^)) where all a^'s are 
odd. In doing this, we are forced to consider the more general series 



<t:t(^) - E 



I3i,l32,-,l3s{^\ \;^ ^fei2"i+622"2+...+b^2"'' 



ni+/3l>n2+,92>--->"s+/3s>0 

where, as before, 61, 62, • • • , &s are positive integers, and (3i, /32, ■ ■ ■ , /3s are integers. 

Lemma 9. For positive integers bi,b2, ■ ■ ■ ,bs and integers /3i, (32, ■ ■ ■ , f3s, the series 
HIj^I^'"1^'{z) can be expressed as a linear combination over 1\z^l'^''\ {for a suitable 
integer e) of the series 1 and series of the form -ffai,a2,.--,ar(-2)5 where all ai's are odd. 
Moreover, in the above expansion of the series Hb-^i,.2,...,bs{^) = -^fei 62 ' bs(^) ^^ have 
e = 0; that is, in that case all coefficients are in TL\z\. 

Proof. We describe an algorithmic procedure for expressing H^^f'^'"'fj'{z) in terms of 
series H2l'2'^'''''2^{z), where either r < s, or r = s and 

max{i : Oj is even or 7^ 7^ 0} < max{i : 6j is even or /3j 7^ 0}. 

In words, in the second case the length of the string of consecutive O's at the tail of the 
upper parameters respectively the length of the string of consecutive odd numbers at 
the tail of the lower parameters has been increased. 

Our algorithmic procedure consists of four recurrence relations, f l3.7p - fl3.10p below. 
For the first two of these, let bg = b'g2'^% where Cs = ^2(65). By definition, the number 
b'g is odd. Then we have 

^bi,b2,-,bs \^) 



E 

ni+l3i-l3s+es>--->ns-i+l3s-i-l3s+es>ns+es>-l3s+es 



J)i2"l+622"2+...+6^_i2"s-i+6;,2"s+'='> 



In the above sum on the right-hand side, let n'^ = ns + Cs be a new summation index. 
Then, for e^ < (3s, one sees that 



'"'^"-'{z). (3.7) 



'^bi,b2,...,bs y^) - 


Tjl3i-l3s+e.s,l32-l3s+es,...,l3s-i-l3s+es,0 /\ 

■ ^bi,b2,-,bs-i,b', y^) 






Ps-es 

+ 2^ ^ Hbub2,...,b.-1 

k=l 


..A-i 


1 the other hand, for Cg > (3s, one has 




'^bi,b2,...,bs y^) - 


tr/3i-/3^+es,/32-/3s+es,.-,/9s-i-/9s+es,0/ X 

■ ^bi,b2,-,bs-i,b's y-^) 






es-l3s-l 

\ ^ b'^2'' TTl3i-/3s+e,-k-l,P2-l3s+es-k-l,. 

k=0 


..,A-i 


Now consider 


rr/3l,...A,0,-,0 f^ 

^bu-,bh,bh+l,-,bsy^)' 





-Ps+Cs-k-l, 



{3.t 
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where 1 < h < s and all of &/i+i, ■ ■ ■ ,bs are odd. Similar to the proceedings above, let 
bh = b'^^^ ^ where th = V2{bh)- Again, by definition, the number b'f^ is odd. Then we 
have 

r7-/3i,...A,0,...,0 ( ^ 
^bu-,bh,b^+r,...M^) 

^ V^ ^6i2"i+---+6h_i2"h-i+6;2"h+'=h+feh+i2"ft+i+-+fe,2"^ 

ni+f5i-/3h+eh>--->nh-i+l3h-i-l3h+eh>nh+eh 
nh+lSh>nh+i>--->ns>0 

In the above sum on the right-hand side, let n'^ = n/^ + e/^ be a new summation index. 
Then, for eh < (3h, one sees that 

ttI3i,...,I3h,0,...,0 / X _ rr'Si-'Sh+eh,.-,/3h-i-/3h.+eh,0,0,...,0/ N 
-"bi,...,6h,6h+i,...,6«l^J ^bi,...,bh-iK'^h+i,-,bs ^^) 

, V^ rr/3i-/3h+eh+fc,.--,^h-i-/3h+eh+fc,fc,0,...,0/ \ /^ q\ 
A^ bi,...,bh-i,b'i^+bh+i2'',bh+2,---,bs ^ '' V • J 

k=0 

On the other hand, for et > (3h, one has 

rr/3i,...,/3h,0,...,0 / X _ rTl3i-l3h+eh,-,(3h-i-l^h+eh,0,0,--;0 ( \ 
^bi,...,bhK+l,-,bsy^' ~ ^bi,...,bh-i,b'^,bh+i,-,bs ^^f 

eh-Ph 

b-^,...,bh-l,b'^^'2''+bh+l,bh+2,■■■,bs ^ '' ^ ' ' 

k=l 

It is clear that, if we recursively apply fl3.7p - fl3.10p to a given series H^^f^'"'fj{z), and 
use H^{z) = 1 as an initial condition, we will eventually arrive at a linear combination of 
1 and series H^f'a^''^.. ar(^) ~ -^01,02,. ..,ar(^) with all a^'s being odd, where the coefficients 
are polynomials in z^^'^" for a suitable e. (Potential fractional exponents come from the 
relation 03.71) .) This proves the first assertion of the lemma. 

Now let us consider the case where all the /3j's are zero. Suppose that we have an 
expansion as described in the first part of the lemma for -^61,62,.. .,6s (-2^); 

Hh,M,...,bsi^) = $^c(a)z^(-)i7,(z), (3.11) 

a 

where the sum is taken over all finite tuples a = (ai, 02, ... ) with all a^'s being odd, 
and where only finitely many coefficients c(a) are non-zero. We also allow the tuple a 
to be the empty tuple () and make the convention that Hq{z) = 1, so that the series 1 
is as well included in the linear combination on the right-hand side of fl3.11l) . 

Let us now consider exponents e(a) that are not integral. Let e be a real number 
strictly between and 1, and concentrate on exponents e(a) with fractional part e; in 
symbols {e(a)} = e. Then we isolate these exponents e(a) in the relation (13. lip , and 
since there are no fractional exponents on the left-hand side, we obtain 

0= Y. c{si)z<''^H,{z). 

{e(a)}=. 
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After dividing both sides through by 2;'^, an apphcation of Corollary [7] shows that c(a) = 
for all a with {e(a)} = e. Thus, all exponents e(a) actually occurring in f l3.1ip with 
non-zero coefficients c(a) are in fact integral. This completes the proof of the lemma. D 

Computer computations suggest that, if we restrict our attention to the series 
Hbi,b2,...,bs{z): which are the ones that we are actually interested in, there is a strength- 
ening of Lemma [9] (see also Appendix Rll. 

Conjecture 10. For any positive integers &i,62, ■ ■ ■ ,bs, the series H},^^b2,...,bs{z) can he 
expressed as a linear combination over Z[z, z^^] of the series 1 and series of the form 
Hai,a2,...,ar{z)i where all ai 's are odd, r < s, and ai + 02 + ■ ■ ■ + a,. < 61 + 62 + ■ ■ ■ + &s- 

To conclude this section, let us provide an illustration of the above discussion. We 
set ourselves the task of determining the coefficient of 2;^°^^^^^^^°^^^ in $^(2;). In order 
to accomplish this task, we first express $^(2;) in terms of series i^f^j ...^^(2;) with all Oj's 
being odd. This is done by means of the expansion (13.11) and the algorithm described 
in the proof of Lemma [91 The resulting expansion is displayed in Appendix El 

Now we have to answer the question, in which of the series -f^ai,...,a^(^) that appear 
in this expansion of $^(z) do we find the monomial 2;^"^^^^^^^"^^^. Using the algorithm 
described in Remark [HI we see that 

1099511640192 = 5 ■ 2^ + 1099511639552, 

= 3 . 2^3 + 2^2 + 2^ + 1099511611392, 

= 2^° + 3 ■ 2^2 + 2^ 

= 2^ + 2® + 3 • 2^ + 1099511639040, 

= 3 . 2^2 + 2^ + 1099511627776, 

= 2^ + 3 • 2^ + 1099511639552, 

= 2^° + 2^3 + 2^2 ^ 27, 

= 1 + 3 -2° + 1099511640188, 

= 3 ■ 2^ + 1099511639808, 

= 1 + 2^ + 2^ + 2° + 1099511640184, 

= 2^^ + 2^2 + 2^ + 1099511627776, 

= 1 + 2^ + 2° + 1099511640188, 

= 2^2 + 2^ + 1099511635968, 

= 2 + 2^ + 1099511640188, 

= 1 + 2° + 1099511640190, 

= 2^ + 1099511640064. 

Here, the third line shows that ^I0995ii640i92 ^ppg^rg j^ Hi^^^i{z), and the seventh 
line shows that it appears in if 1^1. 1^1 (2;) (thereby making it impossible to appear in 
Hi^i.i,i,i{z)) , while the remaining lines show that it does not appear in any other term 
in the expansion of $^(z) displayed in Appendix [XI Hence, by taking into account the 
coefficients with which the series ifi 3 1(2;) and if 1^1^1^1(2;) appear in this expansion, the 
coefficient of ^1099511640192 -^ ^5(^^) jg gggj^ ^Q gq^^l '^43 + 240 = 200. 
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4. The method 
We consider a (formal) differential equation 

V{z-F{z),F\z),F"{z),...,F^^\z)) = Q, (4.1) 

where P is a polynomial with integer coefficients, which has a power series solution 
F{z) with integer coefficients. In this situation, we propose the following algorithmic 
approach to determining the series F{z) modulo a 2-power 2^'^", for some positive 
integer a. We make the Ansatz 

F{z) = Y^ ai{z)^\z) modulo 2=^-^", (4.2) 

i=0 

with ^{z) as given in (II. ip . and where the aj(z)'s are (at this point) undetermined 
Laurent polynomials in z. Now we substitute (14. 2 p into (I4.ip . and we shall gradually 
determine approximations ai^p{z) to ai{z) such that (14. ip holds modulo 2^, for (3 = 

1, 2, . . . , 3-2". To start the procedure, we consider the differential equation (14. ip modulo 

2, with 

F{z) = Y^ ai^i{z)^\z) modulo 2. (4.3) 

1=0 

Using the elementary fact that ^'{z) = 1 modulo 2, we see that the left-hand side of 
(14. ip is a polynomial in ^{z) with coefficients that are Laurent polynomials in z. We 
reduce powers $'^(z) with k > 2""*"^ using the relation (which is implied by the minimal 
polynomial for the modulus 8 given in Proposition [2]) □ 

{^'^{z) + 6^^{z) + {2z + 3)^^{z) + {2z + 6)^{z)+2z + 5z'^Y" =0 modulo 23-2". (4.4) 

Since, at this point, we are only interested in finding a solution to (14. ip modulo 2, the 
above relation simplifies to 

^^"''" (z) + $2"+^ (z) + z^"^' = modulo 2. (4.5) 

Now we compare coefficients of powers $^(z). A; = 0, 1, ... , 2°""*"^ — 1 (see Remark [TTIl . 
This yields a system of 2""'"2 (differential) equations (modulo 2) for the unknown Laurent 
polynomials 0,^1(2;), i = 0,1, . . . , 2°''^'^ — 1, which may or may not have a solution. 



Actually, if we would like to obtain an optimal result, we should use the relation implied by a 
minimal polynomial for the modulus 2^'^ in the sense of Section [51 But since we have no general 
formula available for such a minimal polynomial (cf. Item (2) of Remark [5] in that section) , and since 
we wish to prove results for arbitrary moduli, choosing instead powers of a minimal polynomial for 
the modulus 8 is the best compromise. In principle, it may happen that there exists a polynomial in 
$(2) with coefficients that are Laurent polynomials in z, which is identical with F(z) after reduction 
of its coefficients modulo 2^'^° , but the Ansatz (|4.2p combined with the reduction (|4.4I) fails because 
it is too restrictive. We are not aware of a concrete example where this obstruction occurs. The 
subgroup numbers of 51/2 (Z) (which we treat modulo 8 in Section[TT]by the method described here, and 
modulo 16 by an enhancement of the method outlined in Appendix |D]) are a potential candidate when 
considered modulo 2^'^ for a > 1. On the other hand, once we are successful using this (potentially 
problematic) Ansatz, then the result can easily be converted into an optimal one by further reducing 
the polynomial thus obtained, using the relation implied by a minimal polynomial for the modulus 

23-2° _ 
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Provided we have already found Laurent polynomials ai^j3{z), i = 0, 1, ... , 2"+^ — 1, 
for some /3 with l</3<3-2" — 1, such that 

na+2_i 



E ^^A^)'^'i^) (4-6) 



i=0 



solves (14. ip modulo 2^, we put 

ai,p+iiz):=ai,p{z) + 2%f,+,{z), ^ = 0, 1, . . . , 2"+^ _ i, (4.7) 

where the 6j,/3+i(z)'s are (at this point) undetermined Laurent polynomials in z. Next 
we substitute 

E a,^+i(^)<l>^(^) (4.8) 

j=0 

instead of F{z) in (^^. Using the fact that $'(2) = Ef=o2"^^""^ modulo 2^+\ we 
expand the left-hand side as a polynomial in ^{z) (with coefficients being Laurent 
polynomials in z), we apply again the reduction using relation (14. 4p . we compare co- 
efficients of powers ^''{z), k = 0, 1, . . . , 2""*"^ — 1 (again, see Remark [TT]) . and, as a 
result, we obtain a system of 2""*"^ (differential) equations (modulo 2^+^) for the un- 
known Laurent polynomials 64^^+1(2;), i = 0,1, . . . , 2°"'"^ — 1, which may or may not have 
a solution. If we manage to push this procedure through until /3 = 3 • 2" — 1, then, 
setting ai{z) = 0,^3. 2^(2;), i = 0,1, . . . , 2°+^ — 1, the right-hand side of (14.21) is a solution 
to (14.11) modulo 2^'^°, as required. 

Remark 11. As the reader will have noticed, each comparison of coefficients of powers 
of ^{z) is based on the "hope" that, if a polynomial in $(2) is zero modulo a 2- 
power 2^^ (as a formal Laurent series), then already all coefficients of powers of $(-2) 
in this polynomial vanish modulo 2^. However, this implication is false in general (see 
Lemma [39] below for the case of modulus 2^ = 16). It may thus happen that the method 
described in this section fails to find a solution modulo 2^ to a given differential equation 
in the form of a polynomial in ^{z) with coefficients that are Laurent polynomials in 
z over the integers, while such a solution does in fact exist. As a matter of fact, this 
situation occurs in the analysis modulo 16 of the subgroup numbers of SL2{1>), see 
Theorem [23 In Appendix [Dl we outline an enhancement of the method, which (at 
least in principle; Appendix [D] treats only the case of the modulus 16 explicitly) allows 
us to decide whether or not a solution modulo a given power in terms of a polynomial 
in ^{z) with coefficients that are Laurent polynomials in z over the integers exists, and, 
if so, to explicitly find such a solution. 

It is not difficult to see that performing the iterative step (14. 7p amounts to solving 
a system of linear differential equations in the unknown functions hi^p+i{z) modulo 2, 
where all of them are Laurent polynomials in z, and where only first derivatives of the 
6i_/3+i(z)'s occur. Solving such a system is equivalent to solving an ordinary system of 
linear equations, as is shown by the lemma below. 

Given a Laurent polynomial p{z) over the integers, we write p''°\z) for the odd part 
\{p{z) — p{—z)) and p^^\z) for the even part h{p{z) +p{—z)) oi p{z), respectively. 
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Lemma 12. Let Cij{z) and dij{z), 1 < i,j < N, andri{z), 1 < i < N, be given Laurent 
polynomials in z with integer coefficients. Then the system of differential equations 

N N 

J]Q,,(z)/,(2;) + ^d,,,(z)/;(2;)=r,(2;) modulo 2, 1 < ^ < iV, (4.9) 

j=i i=i 

has solutions fj{z), I < j < N, that are Laurent polynomials in z over the integers if, 
and only if, the system of linear equations 

N N N 

J2^iz)f^{z) + Y.c^{z)ff{z) + Y.z-'d^{z)ff{z) = rf\z) modulo 2, 
j=i i=i i=i 

N N N 

E^S(^)/f ^w + E^Sw/f (^) + E^"'<^5w/f w = ^?^w ^°duio 2, 
j=i j=i j=i 

l<i<N, 
(4.10) 

has a solution in Laurent polynomials /• (z), /• (z) in z over the integers for I < j < 

N. 

Proof. We write fj{z) = fjiz) + fj{z), and observe that 

f'^{z) = z-^ff\z) modulo 2. 

If this is used in 04. 9 1) . and if we separate the even and odd parts on both sides of the 
equations, then JKWf with f!j^\z) = ff\z) and ff\z) = f^°\z), j = 1,2,...,N, 
results after little manipulation. 

Conversely, suppose that gj (z), g^ (z), j = 1,2,..., iV, is a solution to the system 
fHTTU]) . that is, 

N N N 

Y.^^9piz) + J2c^]{z)gf\z) + Y,z-'d!t}{z)gf\z) = rf\z) modulo 2, 
i=i i=i i=i 

N N N 

Y.^f}{z)gf\z) + Y,^}{z)gf\z) + Y.z-'d!f^{z)gf\z) = rf\z) modulo 2, 
i=i i=i i=i 

l<i<N. 
(4.11) 

At this point, the g^ (2;)'s need not be even Laurent polynomials, and the g): {zYs 
need not be odd Laurent polynomials. We have to prove that there exists a solution 
fj {z), fj (z), j = 1,2,..., A^, such that all /j (z)'s are even Laurent polynomials and 

(2) 

all /• {zYs are odd Laurent polynomials. 
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By separating even and odd parts of the g^ (2;)'s and the g^ (-2)'s, we obtain the 
equations 

N N N 

Y.^{^^r^^ +j:^h^)(9fr{z) +j:z-'d^{z){gfriz) = rr\z) 
j=i j=i j=i 

modulo 2, 

(4.12) 

N N N 

j=i j=i j=i 

modulo 2, 

TV iV TV 

T.^i^KafY'^iz) +Y.^]{z){gfr\^) +T.'~'di:H')(9fY"'i^ = -1°^^) 
j=i j=i j=i 

modulo 2, 

(4.13) 

TV TV AT 

i=i i=i i=i 

modulo 2, 

l<i<N. 

Combining KV^ and (|4A3|) . we see that {gf^Y^^\z),{gfY°\z), j = 1,2,...,N, is 

a solution to (14.101) . and now the {g, y'^\zys are indeed even polynomials while the 

i9j Y°\zys are indeed odd polynomials. Addition of both sides of (14.121) and (I4.13P 
then yields that 

/,(z) := ig^'^Y^\z) + {gft\z), 1 < J < AT, 

is a solution to (14.91) in Laurent polynomials in z over the integers. D 

In general, it is difficult to characterise when the system (I4.10p has a solution. What 
one has to do is to solve the system over the field of rational functions in z over Z/2Z, 
and then to see whether possibly occurring denominators cancel out or, in the case of a 
parametric solution, whether denominators can be made to cancel by a suitable choice 
of the parameters. One simple case, where a characterisation is possible, is given in 
Lemma [221 which is crucial for the proof that the generating function for the number 
of subgroups of index n of PS'L2(Z), when these are reduced modulo a given power of 
2, can always be expressed as a polynomial in ^{z) with coefficients that are Laurent 
polynomials in z. 

We remark that the idea of the method that we have described in this section has 
certainly further potential. For example, the fact that the series ^{z) remains invariant 
under the substitution z ^ z"^ (or, more generally, under the substitution z — )■ z^ , 
where h is some positive integer) — up to a simple additive correction — can be 
exploited in order to extend the range of applicability of our method to equations 
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where we not only allow differentiation but also this kind of substitution. This is actually 
already used in a very hidden way in Section [Tl](cf. f27\ Theorem 12], setting in relation 
subgroup numbers of the Hecke group S){q) with subgroup numbers of Cg*Cq modulo 2; 
in terms of generating functions, the meaning of this theorem is that the generating 
function for the former numbers can be expressed in terms of the generating function for 
the latter numbers by a relation which involves a substitution z — )■ z'^). Furthermore, as 
we already mentioned in the introduction, there is no obstacle to modifying the method 
presented here to work for recursive sequences which are reduced modulo powers of p, 
in connection with the series ^„>o ^^"; although at present we are not able to offer any 
interesting applications in this direction. 



5. A SAMPLE APPLICATION: CATALAN NUMBERS 

The Catalan numbers, defined by Cat„ = 7^(^^), n = 0, 1, ... , are ubiquitous in 
enumerative combinatorics. (Stanley provides a list of 66 sequences of sets enumerated 
by Catalan numbers in [5^ Ex. 6.19], with many more in the addendum J23J-) Recently, 
there have been several papers on the congruence properties of Catalan numbers modulo 
powers of 2, see [HI [221 [301 [38]. In particular, in |22] the Catalan numbers are deter- 
mined modulo 64. As we already mentioned in the introduction, the corresponding 
result (cf. [22l Theorems 6.1-6.6]) can be compactly written in the form (11.21) . Clearly, 
once we know the right-hand side of (11.21) . the validity of the congruence (11.21) can be 
routinely verified by substituting the right-hand side into the well-known functional 
equation (cf. ^ (2.3.8)]) 

zC\z)-C{z) + 1 = 0, (5.1) 

where C{z) = Xl^o Cat„ z^ denotes the generating function for the Catalan numbers, 
and reducing powers of ^{z) whose exponent exceeds 7 by means of the relation (14. 4p 
with a = 1. We shall now demonstrate that the method from Section [H allows one not 
only to find the congruence (ll.2p algorithmically, but also to find analogous congruences 
modulo arbitrary powers of 2[3 

Theorem 13. Let ^{z) = J^nyo^"^" ' ^'^'^ ^^^ ^ ^^ some positive integer. Then the 
generating function C{z) for Catalan numbers, reduced modulo 2^'^", can be expressed 
as a polynomial in $(z) of degree at most 2"+^ — 1 with coefficients that are Laurent 
polynomials in z over the integers. 

Proof. We apply the method from Section [H We start by substituting the Ansatz (14. 3 p 
in (15.11) and reducing the result modulo 2. In this way, we obtain 



z 

i=0 i=0 



J2 «JiW'^^'W+ Yl ai,i{z)¥{z) + l = modulo 2. (5.2) 



^In principle, one could use the generalisations of Lucas' theorem due to Davis and Webb 0, and 
to Granville [15], respectively, to analyse the classical expression ;j^(^^) for the Catalan numbers 
modulo a given 2-power, or, more generally, the right-hand side of (|13.2p . But this approach would 
be rather cumbersome in comparison with our method, and it is doubtful that one would be able to 
derive results which are of the same level of generality as Theorems [131 HH or |331 
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We may reduce $^*(z) further using the relation fl4.5p . This leads to 

2a+l_i 2"-l 

1=0 j=0 

2"-l 2°+2-l 

+z Y^ {a^i+2"+i^i{z) + 0^+3.2^^1 (-2)) $^'+^°^'(z)+ ^ ai,i(2;)<l>X-2;) + l = modulo 2. 

i=0 j=0 

(5.3) 

Now we compare coefficients of ^^{z), for i = 0, 1, . . . ,2""*"^ — 1. For i odd, we see 
immediately that this implies that ai^i{z) = modulo 2. Proceeding inductively, we 
now suppose that a2fiu,i{z) = modulo 2 for odd u and some positive integer /3, /3 < a. 
Reading off coefficients of ^^ *, where i is odd, we then obtain 



2 { \ [ 2"+Hl 2 f \ [ 2"+^+l 2 / \ 

^^2/5i,ll^J + ^ ^2''j+2«+l,ll^J + ^ '^2'3j+3-2«,ll^J 

+ -20-2/9 i+2«,l(^) + -2«2'3i+2«+l,l(-2^) + '^2/3+ii,l(-2) = modulo 2. 

However, due to our inductive assumption, all squared terms on the left-hand side of 
this congruence vanish, and we conclude that a2i^+ii^i{z) = modulo 2. 

So far, we have found that all coefficient Laurent polynomials a^ 1(2;) vanish modulo 2 
except possibly 00,1(2;) and 0211+1 i(z). The corresponding congruences that we obtain 
from extracting coefficients of ^^{z) and <l>^° (z), respectively, in (15.31) . are 

za^i{z) + z'^ ^^ a\a+i i{z) + aQ^i{z) + 1 = {i modulo 2, (5.4) 



^-0,1 

2 



zago+i i{z) + 02^+1, i(;z) = modulo 2. (5.5) 



-1 



The only solutions to (15.51) are 020+1, i(z) = modulo 2, respectively a2a+\^i{z) = z 
modulo 2. The ffist option is impossible, since it would imply that, modulo 2, the series 
C(^z) reduces to a polynomial; a contradiction to the well-known fact (easily derivable 
from Legendre's formula ^21, p. 10] for the p-adic valuation of factorials; cf. (12.101) ) that 
the Catalan number Cat^ is odd if, and only if, n = 2^^ — 1 for some h. Thus, 

02^+1, i(-2) = z~^ modulo 2. 

Use of this result in (15. 4p yields the congruence 

zal-y{z) + ao,i(-2) + 22"+'-i + 1 = modulo 2 (5.6) 

for ao.i(z). We let 



ao,i(z) =00,1(2:) + Xl^^ ^ 

fc=0 

and substitute this in (15. 6p . Thereby, we get 

a a 

zal^^{z) + ^z^'"^'"^ + ao,i(-2) + X]'^^''~^ + ^2-+i-i _^ 1 = modulo 2, 

fc=0 k=0 

or, after simplification, 

za^i{z) +ao,i(z) = modulo 2. 
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Again, either 00,1(2;) = modulo 2, or ao,i(-2) — ^^^ modulo 2. Here, the second option 
is impossible, since it would imply that C{z) contains a negative 2;-power, which is 
absurd. 

In summary, we have found that 

a 

o-o^iz) = y^z"^ ~^ modulo 2, 

fc=0 

a2a+i^i{z) = z~^ modulo 2, 

with all other ai^i{z) vanishing, forms the unique solution modulo 2 in Laurent polyno- 
mials ai^i{z) to the system of congruences resulting from fl5.3p . 

After we have completed the "base step," we now proceed with the iterative steps 
described in Section m We consider the Ansatz (I4.6l) - fl4.8l) . where the coefficients ai^^i^z) 

are supposed to provide a solution Cpi^z) = ^j^q ~ '^i,i3{z)^\z) to (15. ip modulo 2^. 
This Ansatz, substituted in (IS.ip . produces the congruence 

2"+2-l 

zCliz) - Cp{z) + 2^ J2 &i,/3+i(^)$'(^) + l = modulo 2^+\ 

i=0 

By our assumption on Cp{z), we may divide by 2^. Comparison of powers of $(2) then 
yields a system of congruences of the form 

hp+i{z) + ^o\i{z) = modulo 2, z = 0, 1, . . . , 2°+^ - 1, 

where Polj(2;), i = 0,1,..., 2"+^ — 1, are certain Laurent polynomials with integer 
coefficients. This system being trivially uniquely solvable, we have proved that, for an 
arbitrary positive integer a, the algorithm of Section H] will produce a solution Cs.2c{z) 
to (15.11) modulo 2^'"^" which is a polynomial in $(-2) with coefficients that are Laurent 
polynomials in z. D 

For example, our computer program needs only about 30 seconds to come up with 
the corresponding congruence modulo 2^'^ = 4096. 

Theorem 14. Let $(2) = Xln>0'^^"- '^hen we have 

Y^ Catn ^" = 2048^^^ + 3072^^3 + 2048^^^ + 3584^^^ + 640^^° + 2240^^ + 32/ 

+ 832^^ + 2412/ + 1042/ + 2702/ + 53z^ + 2z^ + z + 1 
+ (2048z^^ + 3840z^° + 2112/ + 2112/ + 552/ 

+3128/ + 2512/ + 4000/ + 3904/) $(2) 
+ (2048^^^ + 3072^^^ + 1536z^° + 1152/ + 1024/ + 4000z^ + 3440^^ 
+3788/ + 3096/ + 3416/ + 2368/ + 288^) ^^{z) 
+ (2048z^^ + 2048z^° + 2304/ + 512/ + 2752z^ + 3072/ + 728 



/ 
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+3528^^ + 1032^3 + 3168^2 + 3456^ + 3904) <^^{z) 

+ (2048^^2 + 3072^1^ + 1024^^° + 2048z^ + 1152z® + 1728^^ + 2272/ + 2464z^ 

48 \ 
+3452/ + 3154/ + 2136/ + 3896^ + 1600 + — 1 ^*{z) 

+ (2048^^° + 2048/ + 1792z^ + 1792z^ + 1088/ + 1536/ 

2272 \ 
+1704/ + 3648/ + 3288/ + 200z + 3728 + J <l>^{z) 

+ (2048^^^1024/ + 1536/ + 3200/ + 2816/ + 1312/ + 3824/ 

+140/ + 592/ + 3692^ + 488 + ) ^''(z) 

z 



+ (2048/ + 2304^^ + 2304/ + 3520^^ + 960/ + 2456z 

4024' 
+2128/ + 2936^ + 1784 + 1 (!>'(z) 



3 



+ (2048z^° + 1024/ + 2048/ + 512z^ + 3968/ + 1088/ + 1888/ 

339\ 

+832/ + 1444^2 + 2646^ + 3258 + ) $^(2) 

z 

+ (2048z^ + 3328/ + 1536/ + 3008/ 

3152' 
+320/ + 2168/ + 1144^ + 3992 + ) $^(2) 

z 

j/ 



+ (2048/ + 3072z^ + 512/ + 1408/ + 2560 

2380 



+3424/ + 3408/ + 1316^ + 3608 + J $'^(z) 

+ (2048^^ + 2048/ + 2816/ + 3072/ + 1856/ 

, 3904\ -n, , 

+2688/ + 1288^ + 3880 + | $^^(2) 

z 



/ 



+ (2048/ + 1024^^ + 3072/ + 2048/ + 1408 

358 



+2624/ + 1440z^ + 224z + 948 + 1 <^'^{z) 

2384' 
+ ( 2048/ + 2048^^ + 3328/ + 2816/ + 1984/ + 384z + 2488 + ) $^^(z) 
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+ ( 2048z^ + 1024/ + 512z^ + 2432/ + 1792/ + 3040^ + 336 + ) ^^^(z) 

2fiQfi\ 

+ ( 2048/ + 768/ + 256^" + 64z + 2752 + $^^(z) 



/ 26Q6\ 
( 2048/ + 768/ + 256^2 ^ 54^ ^ 2752 + 1 



modulo 4096. (5.7) 



The reader is reminded that coefficient extraction from an expression such as the 
one on the right-hand side of (15. 7p is straightforward, via the algorithm described in 
Section [3] (see (13. ip and the proof of Lemma [9]). 

6. A NON-EXAMPLE 

Consider the equation 

zF%z)-F{z) + 1 = 0, (6.1) 

which has a unique formal power series solution F{z). We note that the coefficients 
in the series are special instances of numbers that are now commonly known as Fufi- 
Catalan numbers, which have numerous combinatorial interpretations; cf. ||2|, pp. 59-60]. 
It was shown in [271 Eq. (36)] that the coefficient of z^ in the series F{z) has the same 
parity as the number of free subgroups of index 14A in the Hecke group io(7) = C2* C-j. 
If we try our method from Section HI then already at the mod-2 level we fail: let 
F[z) = ai{z)^{z) + ao{z) modulo 2, for some Laurent polynomials ao{z) and 01(2;). 
Upon substitution in (16. ip and simplification using (cf. Proposition [2]) 

$^(2) + $(2) + 1 = modulo 2, 

we obtain 

zaQ{z) + zaQ{z)al{z) + ao{z) + zaf{z) + 1 

+ $(2;) [zaQ{z)al{z) + zal{z)ai{z) + ai(z)) = modulo 2. 
or, equivalent ly, 

za^i^z) + zaQ{z)al{z) + zal{z) + ao{z) + 1 = modulo 2 (6.2) 

zaQ{z)al{z) + zal{z)a'l{z) + ai{z) = modulo 2. (6.3) 

However, this congruence has no solution in Laurent polynomials ao{z) and 01(2;). For, 
the Laurent polynomials al{z), aQ{z)a1{z), a1{z), all of them being squares, contain only 
even powers of z when the coefficients are reduced modulo 2. Consequently, the term 
ao{z) + 1 on the left-hand side of (16. 2p can only contain even powers (modulo 2). In 
particular, ao{z) must contain the term 1. If we now suppose that ao{z) and/or ai{z) 
contain negative powers of z, then we obtain a contradiction regardless whether the 
orders (the minimal e such that z"^ appears in a Laurent polynomial) of ao(^) and ai{z) 
are the same or not. This implies that both ao{z) and 01(2;) are actually polynomials 
in z, with aQ{z) being of the form 09(2;) = 1 + 00(2;), where aQ{z) is a polynomial without 
constant term. If we now multiply both sides of (16.21) by al{z) and both sides of (16. 3p 
by al{z), and subsequently add the two congruences, then we obtain 

zal{z) + ao{z)al{z) + al{z)ai{z) + 0.1(2;) = modulo 2. 
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Dividing by ai{z) and replacing aQ{z) by l + aQ{z), we obtain the equivalent congruence 

za[{z) +'dQ{z)ai{z) +'a^{z) + 1 = modulo 2. 

This congruence has no solution since ao(-2) has no constant term modulo 2. 

In the next theorem, we reveal the deeper reason why our method must fail for F{z). 
Namely, it shows that exponents e of terms ^^ which survive in F{z) after reduction of 
its coefficients modulo 2 may have an arbitrary number of blocks of consecutive I's. In 
contrast, a polynomial in $(^) of degree d with coefficients that are Laurent polynomials 
in z can only have terms z^, where e contains at most d blocks of consecutive I's, apart 
from a right-most block of bounded lengthjj 

Theorem 15. Let F{z) be the unique formal power series solution to the functional 
equation (16. ip . Then the coefficient of z"' in F{z) is odd if, and only if, the sequence 
of binary digits of n is built by concatenating (in any order) blocks o/OOll and 01. In 
particular, the number of free subgroups of index n in S){7) = C2*Cj is odd if, and only 
if, the above condition holds. 

Proof. By replacing F[z) by 1 + G{z) in (16. ip . we obtain 

z {I + G{z)f - G{z) = 0, 

or, equivalently, 

G{z) 



{i+G{z)r 

so that G{z) is the compositional inverse of the series z/{l + z)^. By the Lagrange 
inversion formula (cf. [32| Theorem 5.4.2 with /c = 1]), we obtain for n>l that 

{z-)F{z) = {z-)G{z) = -{z-')^-^-^ 

= ^-(z--^){l + zr = U ^" U^f6n + 1 
n ^ ' n\n — IJ Qn + 1\ n 

By the well-known theorem of Legendre [2T| p. 10] (cf. (I2.10p ). we see that the coefficient 
of 2;"" in F{z) is odd if, and only if, 

s{Qn + 1) - s{bn + 1) - s{n) = 0, (6.4) 

where, as in Section [2l s{m) denotes the binary digit sum of m. Another way to phrase 
(16. 4p is to say that, whenever we find a 1 in the binary expansion of n, then there must 
also be a 1 in the binary expansion of 6n -(- 1 at the same digit place. 

We are now ready to establish the claim of the theorem. In view of the above 
considerations, it suffices to show that the condition on n in the statement of the 
theorem is equivalent to (16. 4p . 

Let n be a positive integer with the property that its binary expansion is formed by 
concatenating blocks of the form 0011 and 01. We prove that (16.41) holds in this case by 
induction on n. It is routine to check that our assertion holds true for ra = 1, 2, . . . , 15. 
Now, let n = 4ni + 1, with some positive integer rii. In other words, the right-most 



The length of this right-most block is in fact bounded by the maximal modulus of an exponent of 
z occurring in a Laurent polynomial coefficient. 
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digits in the binary expansion of n are 01 and the binary expansion of rii is formed by 
concatenating blocks of the form 0011 and 01. In that case, we have 

s(6?7, + 1) — s{bn + 1) — s{n) 

= s (4(6ni + 1) + 2 + 1) - s (4(5ni + 1) + 2) - s(4ni + 1) 
= s(6ni + 1) + 2 - s(5ni + 1) - 1 - s{ni) -1 = 0, (6.5) 

by the induction hypothesis apphed to rii. On the other hand, ii n = 16ni + 3 for some 
positive integer rii, that is, if the right-most digits in the binary expansion of n are 0011 
and the binary expansion of ni is formed by concatenating blocks of the form 0011 and 
01, then 

s{6n + 1) — s(5n + 1) — s{n) 

= s (16(6^1 + 1) + 2 + 1) - s (16(5^1 + 1)) - s(16r2i + 2 + 1) 
= s(6ni + 1) + 2 - s(5ni + 1) - s{ni) -2 = 0, (6.6) 

establishing again the truth of (16. 4p . 

In order to prove the converse, let us suppose that n satisfies (16. 4p . We start by 
showing that the binary expansion of n cannot contain any of the substrings 111, 000, 
1011. We call an occurrence of any of these substrings a "violation." 

Assuming that the right-most violation is a substring of the form 111, we have 

n= ...1110..., 

2n= ...1110...0, 

4n= ...1110... 00, 

6n + l = ...101 ...1, 

since to the right of the substring 111 in n there are only blocks of the form 0011 
and 01 according to our assumption, which implies that there cannot be any carries 
"destroying" the substring 101 in 6n+l. However, this means that at the place where we 
find the bold-face 1 in (the binary expansion of) n we find a in (the binary expansion 
of) 6n + 1, a contradiction to (16. 4p . 

Now we assume that the right-most violation is a substring of the form 000. In that 
case, we have 



n = 


. . . 1000 . . . , 


2n = 


...1000...0, 


An = . 


..1000... 00, 


6n + l = 


...10 ...1, 



a contradiction to (16. 4 p for the same reason. 
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Finally we assume that the right-most violation is a substring of the form 1011. Then 
we have 



n = 


...10110..., 


2n = 


... 10110... 0, 


An = . 


..10110... 00, 


M = 


...0001 ...1, 



since to the right of the substring 1011 in n there are only blocks of the form 0011 
and 01 according to our assumption, which implies that there cannot be any carries 
"destroying" the substring 0001 in 6n + 1. However, this means that at the place where 
we find the bold-face 1 in ra we find a in 6ra -|- 1, again a contradiction to (16. 4p . 

Now let us suppose that n = 2"ni+no, where Hq is an a-digit (binary) number formed 
by concatenating blocks of the form 0011 and 01. By applying the computations (16. 5p 
and (16. 6p (possibly several times), one sees that n satisfies (16. 4p if, and only if, rii does. 

We claim that rii cannot be even. For, if it were, say rii = 2^n2, then rii has a 1 at 
the digit place /3 while 6ni + 1 has not, a contradiction to (16.40 . But then the already 
established fact that n, and hence rii, cannot contain any of the substrings 111, 000, 
1011 implies that the right-most digits in the binary expansion of rii form either a block 
01 or a block 0011. This provides an inductive argument that the binary expansion of 
n is formed by concatenating blocks of the form 0011 and 01, and thus completes the 
proof of the theorem. D 

7. Free subgroups in lifts of Hecke groups 

For integers m, g with m > 1, g > 3, and q prime, we consider the group Tm{q) as 
defined in (II. 3p . Denote by f^ (m) the number of free subgroups of index 2qm\ in 
Tmiq)- The purpose of this section is to estimate the 2-adic valuation of f^\ni) in the 
case when g is a Fermat prime. This estimate is based on a recurrence relation for these 
numbers, which, in turn, results from a specialisation of a differential equation in [271 
Sec. 2]. Moreover, this differential equation for the generating function of free subgroup 
numbers in Tm{q) will become of crucial importance in Sections [HI and [T3l 

In order to present the aforementioned differential equation, we first need to compute 
several important invariants of Tm{q)- Using notation and definitions from [2?i Sec. 2], 
we have mv^{q) = 2gm, xO^miq)) = -|^, and thus, 

/i(r„(g)) = 1 - mr^(g)x(r„(g)) =q-l = fi{Sj{q)). 

Moreover, for the family of zeta-invariants {CK{^m{(l)) '■ i^ \ 2qim} oiVmiq)-, we find that 



U'^miq)) = < -1, ^^ = 2gm, (7.1) 




Our first result in this section compares the A- invariants of Tm{q)i as defined in 
Eq. (14)], to those of the underlying Hecke group i3(g), and provides an estimate for 
their 2-adic valuation. 
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Lemma 16. (i) For an integer m > 1 and a prime q > 3, we have 

A^(r^(g)) = m''-%(io(g)), < fi < q - 1. 

(ii) We have 

V2{A^{Sj{q)) > fi, 0<fi<q-l. 
Proof, (i) In view of ( 17. ip . we have, for < // < g — 1, that 

A^iTUq)) = ^ X^(-l)^~^' (^) 2gm(j + 1) [ J] (2gmj + A;) (2gm(j + 1))"' 

^' j=0 ^■' '' \ l<k<2qm J 

{k,2qm)=m 



i=0 ^-^ ^ l<k<2gm 

{k,2qm)=m 






771 



i=0 ^-^ ^ l<fc'<2g 

(fc',2g)=l 



^' 1=0 ^"'^ l<fc'<2(; 



i=0 ^-^ ^ l<fc'<2ij 

(A:',2g)=l 



m'?-iA^(i3(g)), 



as claimed. 

(ii) This follows from [271 Lemma 1]. D 

Our next result is a recurrence relation for the subgroup numbers f^ (m). 
Lemma 17. For m, q as in Lemma [TBI and X > 1, we have 

/SiM=i:i: e e C '^ jH(2gr)-w— ^A,(i,(g)) 



/i=l !^=1 /^i,...,/^i/>0 Ai,...,A^>0 

MiH hMi.=MAiH (-Ay=A-/j 



wi/i initial value f{^ (m) = m'^~'^Ao{Sj{q)). 
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Proof. Setting = Tm{q) in |17l Eq. (18)], and using Part (i) of Lemma [TBI to compute 
-^fj.{^m{(l)) in terms of A^{S){q)), leads to the differential equation 



(7.3) 



for the generating function Gm{q', z) := X1a>o f\+ii^^^- Comparing the coefficient of 
z^ in (I7.3p for A > 1 yields (17. 2p : while, for A = 0, we obtain the required initial value 

f['\m). U 



Given these preparations, we can now show the following estimate for the 2-adic 
valuation of f^y^) i^ the case when g is a Fermat prime. 

Proposition 18. (i) Let m > 1 be an integer, and let q > 3 be a Fermat prime. Then 
we have 

V2{f?{m)) > V2{m){X + g - 2), A > 1. (7.4) 



In particular, if m is even, then f^ [m) is zero modulo any given 2-power for all suffi- 
ciently large values of X. 

(ii) For q = 3, equality occurs in Inequality (17.41) if, and only if , X + 1 is a 2-power. 

Proof, (i) Since (I7.4p is trivially true for m odd, we may suppose that f 2(771) > 0. We 
use induction on A. For A = 1, we have 



V2Ui'\rn)) = V2{m'^-'Ao{Sj{q))) > (g - l)v2{m) = (A + g - 2)v2{m), 



as desired. Now suppose that our claim (I7.4p holds for all f-f (m) such that 7 < L, 
with some integer L > 2, and consider an arbitrary summand 



S = S{n, u,ni,...,n^,Xi,..., X^) 
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in the recurrence relation (17. 2p with A = L — 1. We find that 
V2iS) > ■^2 ( ) - V2ii^\) -v+{q-v- l)v2{m) 



/il, . . . ,/i, 






> f2 ( ) - t;2('^!) - z/ + (g - z/ - l)v2{m) + ix 

^^i, • • • , /^t^y 



i=i 

>V2[ ) - t^2(i^0 - // + (g - // - l)t;2("^) + /i 

\/il, . . . ,/iI,/ 

+ (L - l)f2(m) + (g - 2)vv2{m) 

>{L + q- 2)f2(m) + (g - 3)z/f2(m) + f 2 ( ) + /i - z/ - ^2(2^!), 

V/ii, . . . ,/ii./ 

where we have used Part (ii) of Lemma dn] plus the induction hypothesis in the second 
step. Since z/ < yU, the desired inequality for the 2-adic valuation of fj^{m) will follow, 
if we can show that 

V2{iy\)<{q-3)uv2{m)+V2( ^ Y (7.5) 

Since g is a Fermat prime, we have g — 1 = 2" for some a > 1. Thus, if z/ < g — 1, then, 
by Legendre's formula for the p-adic valuation of factorials (cf. [21, p. 10]), we get 

j>0 i>l 

and (17.51) holds, the left-hand side already being compensated by the term 

(g — 3)z/t>2(m) > g — 3. 

On the other hand, for z/ = g — 1, we have /i = z/ = g — 1, ;Ui = --- = /ig_i = 1, 
f2(z^!) = g — 2, and 

V2( ^ )=V2i{q-l)\) = q-2, 

and the desired conclusion holds again. We have thus shown that every summand S on 
the right-hand side of (17. 2p satisfies V2{S) > {L + q — 2)^2 (m), which implies that 

V2{fi:\m))>{L + q-2)v2{m), 

completing the induction. 

(ii) For g = 3, the recurrence relation (17. 2p . with A replaced by A — 1, takes the form 

/f(m) = 6mA/f_\(m)+ ^ /f(m)/f(m), A > 2, (7.6) 

fl+l/=X—l 



30 M. KAUERS, C. KRATTENTHALER, AND T. W. MULLER 

with initial value /} (m) = 5m?. In order to establish our second claim, we rewrite 
Equation fl7.6p as 



d3)f^\ _fi_\/(3) 



f^>{m) = 6mXf^\{m) + 



2ES/i'H^),/lVi("^)' A = 0(2), 

2E3 /f M/fViM + (/^M)'' A = 1 (2), 

for A > 2, (7.7) 

and argue again by induction on A. For A = 1, Inequality f l7.4p is sharp, as required. 
Now suppose that, for A < L with some L > 2, Inequality fl7.4p with g = 3 is sharp if, 

/ON 

and only if, A + 1 is a 2-power, and consider fl (m) as given by (17.70 . Setting m = 2"m' 
with m' odd, we have 

V2iQmLfl^l,{m)) = l + a + v^^L) + v^UflM)) > a{L + 1) + 1. 

Consequently, if L is even then, by what we have already shown, 

t^2(/f (m))>a(L + l) + l. 

For L odd, all terms except possibly [fi-i i'm)) are divisible by 2"^^+^)+^ or a higher 
2-power, while this exceptional term satisfies 

V2{f^^{m)Y>aiL + l), 

2 

with equality occurring (according to our induction hypothesis) if, and only if, -^yi + l = 
2"^ for some 7 > 1; that is, if, and only if , L + 1 = 2'^'^'^ is a 2-power. This completes 
the induction, and the proof. D 

8. Free subgroup numbers for lifts of the inhomogeneous modular 

GROUP 

In this section, we investigate the behaviour of the numbers /| (m) of free subgroups 
in lifts of the inhomogeneous modular group PSL2{'L) = Sj{3) modulo powers of 2. As 
mentioned in the introduction, the best previous result available in the literature is |2S1 
Theorem 1], which determines the behaviour of /| (1) modulo 16. The results in this 
section solve the problem of determining f^ (m) modulo powers of 2 not only for m = 1 
and the 2-power 2^ = 16, but for all m and modulo any power of 2. 

Let Fm{z) := 1 + X]a>i f\ ("^) ^^ be the generating function for these numbers. (In 
the notation of the previous section, Fm{z) = 1 + zGm{3] z).) By speciahsing g = 3 in 
(17. 3p . one obtains the differential equation 

(1 - (6m - 2)z)F„Xz) - 6mz^F^{z) - zF^{z) - 1 - (1 - 6m + 5m^)z = 0. (8.1) 

Theorem 19. Let ^{z) = X]n>0'^^"' ^^^ ^^^ '^ ^^ some positive integer. Then, for 
every positive integer m, the generating function Fm{z), when reduced modulo 2^'^", can 
be expressed as a polynomial in $(2;) of degree at most 2"+^ — 1, with coefficients that 
are Laurent polynomials in z over the integers. 
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Proof. In view of Proposition [THl the assertion is trivially true for even m, the polyno- 
mial in ^{z) being a polynomial of degree zero in this case. We may thus assume from 
now on that m is odd. 

We apply the method from Section HI We start by substituting the Ansatz (14. 3 p in 
( 18. ip and reducing the result modulo 2. In this way, we obtain 

Y^ ai,i{z)^'{z) + z Y^ al-^{z)^^\z) + 1 = modulo 2. 

i=0 i=0 

This congruence is identical with the congruence (15.21) . Hence, we can copy the resulting 
solution from there. Namely, the unique solution to (15. 2 p (and, hence, to the above 
congruence) is given by 

a 

O'0,i(-2) = /.-^^ ~^ modulo 2, 

fc=0 

a2«+i,i(-2) = z^-^ modulo 2, 

with all other aj^i(z) vanishing. 

After we have completed the "base step," we now proceed with the iterative steps 
described in Section HI We consider the Ansatz (I4.6p - (l4.8p . where the coefficients aj,/3(-2) 

are supposed to provide a solution Fm,i3{z) = ^j^q ~ ctj,/3(-2)'^*(-2) to (18. ip modulo 2^. 
This Ansatz, substituted in (18. ip . produces the congruence 

2°+2-l 



2^ J2 bi,p+i{z)¥{z) + {l-{6m-2)z)F,nA^) 



j=0 



- Gmz^F^^z) - zF^^z) - 1 - (1 - 6m + 5m'^)z = modulo 2^+\ 

By our assumption on Fm^isiz), we may divide by 2^. Comparison of powers of $(-2) 
then yields a system of congruences of the form 

bi,(3+i{z) + Po\i{z) = modulo 2, ^ = 0, 1, . . . , 2°+^ - 1, 

where Polj(2;), i = 0, 1, . . . , 2°"*"^ — 1, are certain Laurent polynomials with integer 
coefficients. This system being trivially (uniquely) solvable, we have proved that, for 
an arbitrary positive integer a, the algorithm of Section H] will produce a solution 
F^ 23-2" (z) to (18. ip modulo 2^'^" which is a polynomial in ^{z) with coefficients that 
are Laurent polynomials in z. D 

We have implemented this algorithm. As an illustration, the next theorem contains 
the result for the modulus 640 

Theorem 20. Let $(-2) = ^„>o^^"- Then, for all positive odd integers m, we have 
1 + V /f ^ (m) z^ = 32z^ + 48/ + 32z^ + (16m + 8)/ + (16m + 8)/ 



A>1 



+ (2m^ + 34) / + (4m^ - 4m + 24) z^ + {5m'^ + 12) z + 1 



To be precise, our implementation finds an expression for each fixed m. These particular results 
can then be put together "manually" into the uniform expression displayed in (j8.2p . 
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+ (48^^ + 2Az^ + 12^2 + mz + 40) ^{z) 

12\ 
16z^ + (16m + 32);^^ + {Am^ - 32m + 68)^^ + ^Qz^ + 22;^ + 12 + — J ^"^{z) 

+ (322^ + 32z^ + (16m - 16)^^ + 40^^ + 4z + 52 + — j $^(2) 

+ (32/ + 32/ + 32z^ + (16m + 24)z^ + (16m + 40)z^ 

35 \ 

+ {2m^ + 16m + 38) 2 + 24 + — 1 ^*{z) 

+ (32;z^ + 16/ + (16m - 8)z + 44) $^(2) 
+ ( 16/ + 16m/ + [Am^ - 16m + 20) ;z + 44 + — J $^(2) 

+ (32^3 + 32^2 + (16m + 16)2 + 40 + - J ^\z) modulo 64. (8.2) 

9. Subgroup numbers for the inhomogeneous modular group 

For a finitely generated group F, let s„(F) denote the number of subgroups of index 
n in F, and write Sy{z) for the (shifted) generating function YlnXd^n+iiX) z"'. 

In this section, we focus on the sequence (s„(PS'L2(Z)))^^^ and its generating func- 
tion S{z) := SpsL2iz)(yZ). We shall show that our method solves the problem of deter- 
mining these subgroup numbers modulo any given power of 2, thus refining the parity 
result of Stothers [35 and the mod-8 result from [29l Theorem 2] mentioned in the 
introduction. 

By the first displayed equation on top of p. 276 in [T^ (cf. also [201 Eq. (5.29)] with 
H = {1} and a = b = h = 1), the series S{z) obeys the differential equation 

(-1 + 4/ + 2/ + 4/ - 2/ - Az^)S{z) + {z' - z^''){S'{z) + S\z)) 

+ l + z + Az^ + Az^ -z^ + Az^ - 2/ - 2/ = 0. (9.1) 



The differential equation (19.11) has a unique solution since comparison of coefficients of 
z^ fixes the initial values, and yields a recurrence for the sequence (s„(PS'L2(Z)))^^_^ 
which computes s„+i(PS'L2(Z)) from terms involving only Si{PSL2{'L)) with i < n. 



Theorem 21. Let ^{z) = Yln>o^'^"^ ^'^^ ^^^ '^ ^^ some positive integer. Then the 
generating function S{z) = SpsL2{z){^): when reduced modulo 2^'^" , can be expressed 
as a polynomial in $(z) of degree at most 2"+^ — 1 with coefficients that are Laurent 
polynomials in z over the integers. 
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Proof. We apply the method from Section m We start by substituting the Ansatz fl4.3p 
in fl9.ll) and reducing the result modulo 2. In this way, we obtain 

2«+2_l /2°+2-l 2"+2-l 

i=0 \ i=0 i=0 

2"+2-l 



Y^ al^{z)^^\z) ] + 1 + z + z^ = modulo 2. 



j=0 



We may reduce $^*(z) further using Relation (14.51) . This leads to 

Y a,i(^)$'(^) + (/ + ^i°) Yl ^(^^A^W-\^)^'i^)+ E <iW*^W 

i=0 \ i=0 i=0 

i=0 i=0 

2"-l \ 

+ E («?+2«+m(-2) + «'+3.2m(-2)) <I>''+'"^' (z) + 1 + z + z^ = modulo 2. (9.2) 



i=0 



In the same way as in the proof of Theorem [19], one sees that all coefficients aj,i(-2) 
vanish modulo 2, except possibly 00,1(2;) and 02^+1,1 (-2). The corresponding congruences 
obtained by extracting coefficients of $°(z) and $^° {z), respectively, in (19.21) . are 

ao,i(-2) + (/ + -2^°) (a'oi{z) + al^z) + z^"^\lo.+iAz)^+l + z+z'^ = modulo 2 (9.3) 

and 

02^+1, i(z) + {z"^ + z-^^) [a2a+i i{z) + a^a+i i{z)) = modulo 2. (9.4) 

The only solutions to (19.41) are 02^+1,1 (,2) = modulo 2, respectively a2o+i^i{z) = z~^ + 
z~'^ modulo 2. The ffist option is impossible, since there is no Laurent polynomial 
ao,i(^) solving the equation resulting from (19. 3p . Thus, we have 

a2Q+i,i(-2) = z'"^ + z~'^ modulo 2. (9.5) 

Use of this result in (19. 3p yields the congruence 

ao,i(z) + /(l + z^) {a'f^Az)+alAz))+z'^"^'-^l+zY + 'i- + z+z^ = modulo 2. (9.6) 

for ao,i(-2). We let 



a 

k 



ao,i{z) = ao^z) + z^^ + 2;"^ + z'^ + (1 + z^) E ^^'"^ 



fc=2 
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and substitute this in f l9.6p . Thereby, we get 

a 

So,i(z) + z\l + z^) {a'^^^{z) + Zl^{z)) + Z-' + z-^ + z-^ + (1 + z^) ^ z""-' 

k=2 

+ z\l + z") iz-^ + z~' + Y^ ^''"' + ^~'' + ^"' + ^'' + (1 + -2')' E '^''^'"'' 

\ fc=2 fc=2 

+ ^2°^'-^(l + z^)^ + l + z + / = modulo 2, 
or, after simphfication, 

o-o^i{z) + 2;'^(1 + 2;'^) (oq 1(2;) + Oq 1(2;)) = modulo 2. 



Again, either ao,i(-2) = modulo 2 or ao^i{z) = z~'^ + z~^ modulo 2. Here, the second 
option is impossible, since it would imply that S{z) contains a negative 2;-power, which 
is absurd. 

In summary, we have found that 



00,1(2) = z~'^ + z~'^ + z~^ + (1 + -2^) ^ 2^'°"^ modulo 2, 

fc=2 

a2c+i 1(2) = z~'^ + 2;"^ modulo 2, 



with all other ai^i{z) vanishing, forms the unique solution modulo 2 to the system of 
congruences resulting from (19. 2p in Laurent polynomials Oj 1(2;). It should be noted that 
all aj.i(2;)'s, 1 <i < 2^" — 1, are divisible by 1 — 2;^ modulo 2, as is ao.i(-2) — 1- 



After we have completed the "base step," we now proceed with the iterative steps 
described in Section m The arguments turn out to be slightly more delicate here than 
in the proof of Theorem [13 To be more precise, when considering the Ansatz (14.61) - 
(14. 8p . where, inductively, the coefficients ai^js^z) are supposed to provide a solution 

Si3{z) = J2i=o ~ '2i,/3(-2)'^*(^) to (19. ip modulo 2^, we must also assume that ai^js^z), 
1 < i < 2^° — 1, and ao^^lz) — 1 — 2z^, are all divisible by 1 — z^. The reader should 
note that the divisibility assumptions do indeed hold for /3 = 1, the term — 2z^ being 
negligible, since for /3 = 1 we are computing modulo 2^ = 2. 
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The above Ansatz, substituted in f l9.ip . produces the congruence 



4 = 



i=0 4=0 

2"+2-l 



,Q + 2_ 



4 = 




2'»+2-l 



+ Y <A^Wi^) + 



4=0 




+ 1 + 2 + 422^4^3-2^ + 4^^-22^-2^^ = modulo 2^+^ (9.7) 



By our inductive construction, we know that the terms contained in hues 3-7 of fl9.7p 
are divisible by 2^. Hence, if we were to divide by 2^ and compare coefficients of $*(2), 
for i = 0, 1, . . . , 2°"^2 — 1, we would obtain the modular differential equations 

bi,^+i{z) + (2^ - 2^°) {b[p+i{z) + {t + l)6i+i,;3+i$'(^)) + Pol4(;2) = modulo 2, 

z = 0,l,...,2°+2_i, (9.8) 

where Polj(2), i = 0, 1, . . . , 2"'*'2 — 1, are certain Laurent polynomials with integer 
coefficients. If i is odd, then the term {i + l)6j+i ^3+1 in (19.81) vanishes modulo 2. Hence, 
in this case, the differential equation (19. 8p is of the form appearing in Lemma [221 
The lemma then says that such a differential equation has a solution if, and only if, 
the Laurent polynomial Polj(2) satisfies the condition given there. We must therefore 
verify this condition for our Laurent polynomials Polj(2), arising through division of 
lines 3-7 of (19. 7p by 2^^. We shall actually prove (see the following paragraph) that 
Poli(2) is divisible by (1 - 2^)2, for all i with 0<i< 2"+^ _ 1. por odd i, Corollary [23] 
thus implies not only unique existence of solutions bi^p+i{z) but also divisibility of these 
solutions by 1 — 2^. If we now consider Equation (19. 8p for even i, that is, 

6^,^+1(2) + (2^ - 2i°)6:^^+i(2) + 2^(1 - 23)6,+i_^+i$'(2) + Poli(2) = modulo 2, 

i = 2,4,...,2°+2-2, 

then we see that divisibility of &4+i,^+i(2) by 1 — 2'^ guarantees that we may again apply 
Corollary [23] to obtain that there is also a unique solution 64/3+1(2) for even i, and that 
this solution is divisible by 1 — 2^. In summary, we would have proved that, for an 
arbitrary positive integer a, the algorithm of Section [H produces a solution 823-2" (2) 
to (19. ip modulo 2^'^°' , which is a polynomial in $(2) with coefficients that are Laurent 
polynomials in 2. This would establish the claim of the theorem. 
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It remains to prove that lines 3-7 of fl9.7p are indeed divisible by (1 — z^)^. In order 
to see this conveniently, we write 

ai^p{z) = di^p{z){l - z^) modulo (1 - ^^)^ z = 1, 2, . . . , 2"+^ _ i, 

and 

ao,/3(^) = -1 - 2z^ + (io,/3(2;)(l - z^) modulo (1 - z^Y, 

where the di^p{zys are polynomials of the form po+Pi -2+^2-2^, for some integers po, Pi, P2- 
(It should be noted that it is at this point where we use our inductive hypothesis on 
divisibility of the coefficients ai^p{z).) Then, reduction of lines 3-7 of (19. 7p modulo 
{1 — z^Y leads to the remainder 

2"+2-l 



(3 + 2^(1-23)) _i .2^2 ^(1-^3) Y^ d,,p{z)^\z) 



j=0 



2Q + 2_]^ 



+ z\l-z'^){ ~4z-3z^ J2 d,^^{z)¥{z) + l + 4z'^ + 4z^ 

\ j=0 

+ 3 + 6^2 + z{l - z^) modulo (1 - z^. (9.9) 

Using the fact that 

z\l - z^) = z{l - z^) modulo (1 - z^f 

and similar reductions modulo (1 — z^Y, one sees that the expression in (19. 9 p is in fact 
divisible by (1 — z^)'^, as claimed. This completes the proof of the theorem. D 

Recall that, given a Laurent polynomial p{z) over the integers, we write p'^°\z) for 
the odd part ^{p{z) — p(—z)) and p^^\z) for the even part ^{p{z) + p{—z)) oi p{z), 
respectively. 

Lemma 22. The differential equation 

a{z) + {z'^ ~z^^)a\z) + Fo\{z) = modulo 2, (9.10) 

where Pol(z) is a given Laurent polynomial in z with integer coefficients, has a solution 
that is a Laurent polynomial if, and only if, Pol'°-*(z) is divisible by 1 + z^ (modulo 2). 
In the latter case, the unique solution is given by 

a(z) = Pol^'^) (z) + l±i-Pol(°) (z) = Pol(^) (z) + i±^^t^Pol(") (z) modulo 2. 
1 + z° 1 — z-^ 

Proof. Let ao{z) = z^ . Then it is obvious that 

ao(^) + {zJ — z}'^)o!q{z) = ao{z) modulo 2 

if m is even, and that 

ao{z) + (z'^ — 2^°)ag(2;) = (1 + z^)ao(z) + ^^00(2;) modulo 2 

if m is odd. The assertion of the lemma follows now immediately. D 

Corollary 23. //, in the differential equation ( ]9.10p . the Laurent polynomial Pol{z) is 
divisible by (1 — z^)"^ (modulo 2), then the uniquely determined solution a{z) is divisible 
by 1 — z^ (modulo 2). 
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We have implemented the algorithm described in the proof of Theorem |2I1 As an 
illustration, we present the result for the modulus 64. 

Theorem 24. Let $(2;) = J2n>o^'^" ■ Then we have 

n>0 

= 32^^° + 48/^ + 48/^ + 32z^*^ + 32^^^ + 32^^^ + 48^^^ + IQz^^ + 40^^^ + IQz^^ 
+ 32^2^ + 32^2^ + 16;z^^ + IQz^^ + 52^^° + 32;z^^ + 40z^^ + 60^^^ 

+ 48z^*^ + Az^^ + 322^^ + 42^2 + 362" + 16;z^° + 60^^ + 2/ + 16^^ + Az^ 

. 4 o , 56 36 51 33 52 1 

+ 60^^ + 44^^ + 16^2 + 54^2 + 60;z + 32 + — + — + — + — + — + — 

z 2;^ Z"^ z^ z^ z' 

+ I 32^^ + 32^2^ + 32^2^ + 32^2^ + 16^^^ + 32^^^ + 32^^° + 32z^^ + 32^^^ 



48^^^ + IQz^^ + 16z^^ + I62" + 32^^° + 32/ + 48/ + 8z^ + 8/ + 48z^ + 24z + 32 

20 12 8 36 4 24\^, , 

+ ( 32/^ + 32^2^ + 32^2^ + 32^2^ + 32^^^ + 32^^^ + 48z^^ + 32^^^ + 48z^^ + 32^^^ 

+ 48^^^ + 32z^2 + 56^^° + 8/ + 16/ + 48z^ + 24z^ + 56^^ + 44/ + 16z^ 

2 60 50 48 8 50 52 52\^2, , 

+ 48/ + 4O2; + 44 + — + — + — + — + — + — + — I <l>^(z) 

Z Z'^ z-^ z^ z° z^ z' 

+ ( 32/^ + 32/^ + 32/^ + 32/° + 32;z^^ + 48z^^ + 32;z^^ + 32^^^ + 32^^^ 

32 20 

+ 32;z^^ + 16^^° + 48/ + 8/ + 482^ + 56/ + 82^ + 16/ + 48^ + 56 + — + — 

z z^ 

52 4 36 12 36\^n, , 

z"^ z^ z° z° z' 

+ ( 32/^ + 32z^^ + 32/^ + 32z^2 + 32/^ + 32z^° + 32^^*^ + 32/^ + 16^^^ + 32/^ 

+ 32/3 ^ 4g^22 ^ ^g^2i ^ 4Q^2o ^ 32^19 ^ 32^18 ^ 24^^^ + IQz^^ + 48^^^ + 32z^^ 

+ IQz^^ + 82^^ + 32^" + 56^^° + 56/ + 44^^ + 40^^ + 48/ + IQz^ + 20/ + 56^^ + 30^^ 

40 34 52 17 26 40 29\ ^4, 

+ 32^ + 28 + — + ^ + ^ + ^ + ^ + -r + ^| ^\z 
Z Z'^ Z"^ z^ z° z° z' 

+ I 32/2 ^ 32_230 ^ 32^26 ^ 32^24 ^ 32^23 _^ 32^22 _^ 32^21 ^ ^3^20 ^ ^g^is ^ 32^16 ^ 43^14 
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+ 32z^^ + 48z^^ + 48z^^ + 32/ + 16^^ + 56/ + 48z^ + 48/ + 40/ + 16z^ 

24 24 20 24 40 20\^,, , 

+ 322; + 56 + — + — + — + — + — + — ]^'iz) 

z z^ z-^ z^ z^ z^ 

+ (32/2 + 32/^ + 32/° + 32/^ + 32/^ + 32/^ + 48;^^^ + 16;^^^ + 48;^^ 

+ 16^^^ + A&z^^ + 322^2 ^ 32^11 + 56/ + 49/ + 55/ + 16/ 

, . n 8 52 60 30 20 20 14\ ., , 

+ 8/ + 56/ + 4/ + 56z + 32 + - + — + — + — + — + — + — ) $6(2) 

z z^ z-^ z^ z^ z° z' 

32/° + 32/^ + 32^21 ^ 32/0 ^ 43^18 ^ 32^16 ^ 43^14 ^ 32^13 ^ 43^10 ^ ^g^9 ^ g/ 

+ 32/ + 16/ + 16/ + 8/ + 48^ + 40 + - + 4 + ^ + ^ + ^ + ^ + ^V'W 

z / / / / / // ^ ' 

modulo 64. (9.11) 

10. Counting permutation representations of 1^(3) for m prime 

Let m be a prime and, for a finitely generated group F, let /ir(^) := |Ifom(r, S'„)|. 
We want to determine the function hY^{2,){n) counting the permutation representations 
of degree n of the lift T„i{'i) of the inhomogeneous modular group PSL2{1,) = S){3). 
Suppose first that m > 5, and classify the representations Tm{3) — > 5'„ by the image 
p E Sn oi the central element a;^ = y^. The permutation p must be of the form 

r 

P = Y\_^i^ < r < [n/m\, 

i=l 

with pairwise disjoint m-cycles (Tj. 

For fixed r, the symmetric group Sn contains exactly 

U "" )^rn-l)V = — ""■ ^, (10.1) 

r! \m, . . . ,m,n — mrj r! [n — mr)\ mJ' 

such elements p. 

Next, given such p, the image of the generator x will contain a certain number s 
of 2?72-cycles in its disjoint cycle decomposition, < s < \r/2\ , each of which breaks 
into two m-cycles when squared. Thus, in order to construct a square root of p (i.e., a 
possible image of x) , we need to 

(i) fix s in the range < s < \r/2\ , 

(ii) select s 2-element subsets from the r ?72-cycles of p, which can be done in 

1 f r\ f 2s \ r\ 



s\\2sj \2,...,2j 2^s!(r-2s)! 
different ways, 

(iii) lift each of these s pairs of m-cycles to a 2m-cycle (whose square is the product 
of the two given m-cycles), which can be done in m ways. 
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(iv) compute a^ for each of the r — 2s remaining ?7i-cycles a, where u is the multi- 
phcative inverse of 2 modulo m, and 

(v) select a permutation vr subject only to the condition that vr^ = 1 on the n — mr 
letters not involved in any of the r m-cycles of p. 
Hence, there are precisely 

r! he, (n - mr) J^ ^ ,.^_^ x, (10-2) 

„_n ■ V /■ 



s=0 

distinct square roots for each p involving r m-cycles. 

Similarly, the number of cubic roots of such p is given by 

l'-/3J 2i 

rl hc.(n — mr) y ^ . , — (10.3) 

fr^ 3*t!(r-3t)! ^ ' 

(classify the cubic roots by the number t of 3?Ti-cycles, and use the fact that each 
product of three disjoint m-cycles is the cube of precisely 2m^ 3m-cycles). Multiplying 
(llU.ip - (llU.3p and summing over r = 0, 1, . . . , [n/m\ , we find that 

[n/m\ Lr/2J [r/Sj 

h { ] = \ \^ \^ \^ r! hcAn - mr)hc-,{n - mr) 

^r„(3)W ^- 2^ 2^ 2^ 2-3*m'^-«-2*s!t!(n-mr)!(r-2s)!(r-3t)!' "^ " 

(10.4) 

The cases where m = 2, 3 have to be treated separately. By arguments similar to the 

ones above, one finds that 

, ,, /^^^^^ {2r)\hc,{n~Ar)hc,{n-Ar) 

/.r.(3)(n) =nll^l^ 22(-^)3V!.!(n-4r)!(2r-3.)!' ^'^'^^ 

and that 

/.r3(3) in) = n\ V v' ^ hc{n - 9r)hc{n - 9r) 

r=0 s=0 ^ ' ^ ' 

Furthermore, it is well-known that, for a prime p, 

^^'^""^ " ^ p'=A;!(n-p/t)!' 
which allows us to make Formulae (I10.4p -( 110.6p completely explicit. 

11. Subgroup numbers for the homogeneous modular group 

In this section we consider the problem of determining the behaviour of the number 
of index-n-subgroups in the homogeneous modular group S'L2(Z) modulo powers of 2. 
By a folklore result that goes back at least to Dey [9J, these subgroup numbers are 
in a direct relation to numbers of permutation representations of 5*^2 (Z). Our start- 
ing point is a recurrence with polynomial coefficients for the latter numbers, which 
is then translated into a Riccati-type differential equation for the generating function 
X^„>QSrt+i(S'L2(Z)) ^'^ of the subgroup numbers. (Equation (Ill.Sp displays this equa- 
tion when reduced modulo 16.) Our method from Section H] is then applied to this 
differential equation. Direct application already fails for the modulus 8. Interestingly, 
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if we instead apply our method with the minimal polynomial for the modulus 16 given in 
Proposition O then our algorithm produces a result for modulus 8 (see Theorems 1261 and 
[271) . but then fails at the level of modulus 16. By the enhancement of the method out- 
lined in Appendixini we are nevertheless able to treat the subgroup numbers Sn{SL2{Z)) 
modulo 16 as well (see Theorem |28|) . In view of the already substantial computational 
effort involved in the case of modulus 16, we did not try to push our analysis further to 
higher powers of 2. In particular, as opposed to the case of the subgroup numbers of 
PSL2{1j), it remains unclear whether it is possible to express the generating function 
^„>Q s„+i(5'L2(Z)) z"', when the coefficients are reduced modulo 2'^, as a polynomial in 
$(z) with coefficients that are Laurent polynomials in z over the integers for all 7 > 1. 
We feel, however, that this should be the case; see Conjecture | 



Let us start with the aforementioned result (cf. [HI Theorem 6.10], see pHl Prop. 1] 
for a conceptual proof, plus generalisations) relating the numbers of subgroups of a 
finitely generated group to the numbers of its permutation representations. 

Proposition 25. Let T be a finitely generated group. Then we have 

00 ^ / 00 ^\ 

J2 I Hom(r, Sn)\^ = exp ^ .„(r)- . (11.1) 






We take P = P2(3) = SLiiZ) and combine (TTTTl) with (TT03D . the latter giving 
an explicit formula for the homomorphism numbers hn '■= | Hom(S'L2(Z), S^)]. Us- 
ing the Guessing package [T6], we found a recurrence of order 30 for the sequence 
{hn/n\)n>o, with coefficients that are polynomials in n over Z. The validity of the recur- 
rence was verified by computing a certificate using Koutschan's Mathematica package 
HolonomicFunctions [19j|j However, this recurrence is not suitable for our purpose, 
for which we require a recurrence with coefficients that are polynomials in n over Z, and 
with leading coefficient n. A recurrence of this form, if it exists, must be a left multiple 
of the recurrence operator corresponding to the minimal order recurrence. The con- 
struction of such left multiples is known as desingularisation, and algorithms are known 
for this purpose [T]. This technique can be used to eliminate factors from the leading 
coefficient of the recurrence (whenever possible), but it cannot be used to ensure that 
the leading coefficient be a monic polynomial. The recurrence of order 32 mentioned in 
Footnote [6l with leading coefficient 1, could indeed be used for our purpose, by simply 
multiplying it by n. However, since this recurrence has high-degree polynomials as coef- 
ficients, we preferred to work with a different recurrence with lower degree polynomials 
as coefficients. The price to pay is that the order of such a recurrence will be higher. So, 
by an indeterminate Ansatz, we computed a candidate for a recurrence of the desired 
form of order 50, with polynomial coefficients of degree at most 5|j To be precise, it is 



The certificate has 4 megabytes, and, to obtain it, required about 30 hours of computation time. 
The coefficients of this recurrence are polynomials in n of degree 34. Interestingly, there is a recurrence 
of order 32 with leading coefficient 1. Although we did not try to prove it, it is likely that the recurrence 
of order 30 is the recurrence of minimal order. 

We used the function LinSolveQ of the Guessing package [TB], which uses modular arithmetic, in 
order to solve the arising system of linear equations. Mathematica' s built-in linear system solver is not 
capable of solving it on current hardware due to the huge numerators and denominators of rational 
numbers which arise during the computation. 
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the uniquely determined recurrence of the form 

^\^a{k,i)n \ --— = 0, n>50, 

fc=0 ^ i=0 ^ ^ '' 

where 

a(0, 0) = a(0, 2) = a(0, 3) = a(0, 4) = a(0, 5) = 0, 

a(0,l) = l, 

a(50, 5) = 47323476536606893277939021 129424044201294092725261226600745838\ 

993897087202045010603943040012232525, 

a(50, 4) = -853333370519051585059335896571817612918194491041969759097679\ 

3078743106989966250706985019403282594096, 

a(49, 5) = 2507660784286104701612089471873568042396155618028516886767837\ 

559764248217845308468763736164634176, 

a(49,4) = a(48,5) = a(48,4) = a(47, 5) = a(47,4) 

= a(46, 5) = a(46, 4) = a(45, 5) = a(45, 4) = 0. 

Subsequently, we checked that this recurrence is a left-multiple of the certified recurrence 
of order 30, thereby establishing validity of this candidate recurrence of order 50. This 
last recurrence was then converted into a linear differential equation with polynomial 
coefficients for the series 

i7(z) := V/i„^ = V|Hom(^L2(Z),S'„)|^. 

n=0 n=0 

Finally, this last mentioned differential equation can be translated into a Riccati-type 
differential equation for the generating function 



S{z) = Y,Sn+i{SL2mz^ (11.2) 



n>0 

for the subgroup numbers of SL2{Z). This is done by differentiating the relation flll.ip . 
with r = SL2{Z), several times and by dividing by H{z). This leads to relations of the 
form 

^^M = p^^s{z), S'{z), . . . ), A; = 1, 2, ... , (11.3) 

where Pk{S{z), S'{z), . . .) is a polynomial in S{z) and its derivatives, which can be 
determined explicitly using the Faa di Bruno formula for derivatives of composite func- 
tions (cf. [HI Sec. 3.4]; but see also [H [H]). Substituting these relations in the linear 
differential equation for H[z), one obtains the announced Riccati-type differential equa- 
tion for S{z). It turns out that this differential equation has integral coefficients, so 
that it is amenable to our method from Section |H The differential equation cannot 
be displayed here since this would require about ten pages. Its reduction modulo 16 
is written out in (lll.5p . Our method from Section H] with a = applied to (I11.5P 
does in fact not produce a result modulo 8 = 2^'^ (it stops at the level of modulus 4). 
If, however, we use the method from Section H] with the minimal polynomial for the 
modulus 16 in place of the one for the modulus 8, then the method goes through up to 
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modulus 8 (but fails for modulus 16). This yields the following theorem. It refines the 
parity result [Sni Eq. (6.3) with \H\ = 1, g = 3, m = 2]. 

Theorem 26. Let $(-2) = Xln>o'^^"- Then we have 

n>0 

= 4:Z^° + 4:Z^^ + 4:Z^^ + 4:Z^^ + 4:Z^^ + iZ^ + 6Z^ + 4:Z^ + 6Z^ +4:Z^ +4:Z+6 

736641 3 6 

^2 ^3 ^6 ^7 -;8 ^9 -;11 -;12 



46662444 6 2 



+ (4^3 + 4^2 + - + ^ + ^ + ^ + -^ + - + - + — + -^ + -^)$(z) 



, s 4 -^ 9 4 6 2 5 2 

+ { Az^ + Az'^ + Az^ + 6z^ + A + - + — + — + — + — 

z z^ z-^ z^ z° 

6 1 4 6 4 6 6 5\2/N 
( . 2 ^ 4 2 4 4 2 4 4 4 2X3,, 

+ 4^2+ + + + + + + + + + $3 2 

\ 2;^ z^ z^ z^ z^ z^ z^ 2;" z^^ 2"y 

modulo 8. (11.4) 

Proof. The Riccati-type differential equation for S{z) (as defined in (lll.2p ) modulo 16 
i^ 

Po{z) + Piiz)S{z) + P2{z)Sizf + Psiz)Sizf + p,{z)S{z)' + p5{z)S{zf + pe{z)S'{z) 

+ Pr{z)S\zf +Psiz)S{z)S\z) + pg{z)S{zfS\z) + pio{z) S {zf S' (z) 

+ p^^iz)Siz)S'{zf+Puiz)S"iz)+p,s{z)Siz)S"iz) + p^.iz) S (zf S" (z) 

+ p^,{z)S\z)S"{z)+p^,{z)S"'{z)+p^r{z)S{z)S"\z)+p^s{z)S""{z) = 

modulo 16, (11.5) 

with coefficients Pj{z), j = 0, 1, . . . , 18 as displayed in Appendix [Bl 

The differential equation (111.51) has a unique solution since comparison of coefficients 
of z^ fixes the initial values, and yields a recurrence for the sequence (s„(S'L2(Z)))^^_^ 
which computes s„+i(S'L2(Z)) from terms involving only Sj(S'L2(Z)) with i < n. 
Now we apply the method from Section H] with the polynomial 

(<l>2(z) + $(2) + z){^^{z) + 6^\z) + {2z + 3)^^{z) + {2z + 6)<^{z) + 2z + 5z^) (11.6) 

in place of the polynomial on the left-hand side of (14.41) to the differential equation (I11.5P 
(that is, in view of Proposition [2] we are aiming at determining the subgroup numbers 
of Sn{SL2{'L)) modulo 16). This yields the above result by means of a straightforward 
computer calculation. D 



We display the differential equation modulo 16 in order to prepare for Theorem [ 
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If we want to know explicitly for which n the subgroup number Sn{SL2{'L)) is con- 
gruent to a particular value modulo 8, then we should first apply the algorithm from 
Section E] (see (13. ip and the proof of Lemma [HD in order to express powers of $(z) on 
the right-hand side of flll.4p in terms of the series -f^ai,...,a,.(-2)- (The corresponding 
expansions are in fact listed in (I2.6P and (12 .yp .) The result is 

■5-^ /444\ /444\ 



n>0 



+ |4--^ + ;i + ;3 +71 + ^5 +;j + T7 + 7o+7T + IS + 7l3l«.aW 



+ 4^^ + 4^3 + 6^2 + 4 + 4/ + - + — + — + 



z z^ z^ z^ 

2 2 1 4 6 4 6 2 5\ 

+ iF + ^6 + ir + is + ^9 + ilo + ilT + il^ + ^ J ^1 ( ^) 

61267225 6 1 3 

+ i + i^ + i^ + i^ + i^ + i^ + i^ + i^ + i^ + i^ + i^ moduloS. (11.7) 
From this expression, it is a routine (albeit tedious) task to extract an explicit de- 
scription of the behaviour of the subgroup numbers of SL^iJ^ modulo 8. Since the 
corresponding result can be stated within moderate amount of space, we present it in 
the next theorem. 

Theorem 27. The subgroup numbers Sn{SL2{Z)) obey the following congruences mod- 
ulo 8 : 

(i) Sn{SL2{1)) = 1 (mod 8) if, and only if, n = 1,2,4,10, or if n is of the form 
2°" — 3 for some cr > 4; 

(ii) Sn{SL2{'L)) = 2 (mod 8) if, and only if, n = 7, 12, 17, or if n is of one of the 
forms 

3 ■ 2'" - 3, 3 ■ 2'" - 6, 3 ■ 2'" - 12, for some a > A; 

(iii) Sn{SL2{Z)) = 4 (mod 8) if and only if, n = 3,22,23,27,46,47,51, or if n is 
of one of the forms 

2" + 6, 2" + 7, 2" + 11, 2" + 12, 2" + 18, 2" + 21, for some a > 5, (11.8) 

2'" + 2^-2, 2'' + 2^ + l, 2'" + 2^ + 3, 

for some cr, r with a > 6 and 4 < r < o" — 1, (11-9) 

2<^ + 2^ + 2*^ - 12, 2^" + 2^ + 2^ - 6, 2^" + 2^ + 2^ - 3, 

for some a, r, u with a > 6, 5<z/<cr — 1, and 3 < r < z/ — 1; 

(11.10) 

(iv) s„(S'L2(Z)) = 5 (mod 8) if, and only if, n = 5, or if n is of one of the forms 

2" - 6, 2'" - 12, for some a > 5; 
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(v) SniSL2{Z)) = 6 (mod 8) if, and only if, n = 6,11,14,18,19,21,33,34,35,37, 
or if n is of one of the forms 

2<^_2, 2''-4, for some a > 5, (H-H) 

2'" + l, 2'' + 2, 2'' + 3, 2'" + 4, 2^" + 5, 2^" + 10, 2"" + 13, 

for some a > 6, (11.12) 

2'' + 2^-3, 2'' + 2^ -6, 2'" + 2^ -12, 

for some a, r with cr > 7 and 5 < r < a — 2; (11.13) 

(vi) in the cases not covered by items (i)-(v), s„(S'L2(Z)) is divisible by 8; in par- 
ticular, Sn{SL2{'L)) ^ 3, 7 (mod 8) for all n. 

As we already said earlier, the method from Section H] with the polynomial in (I11.6P 
in place of the polynomial on the left-hand side of (14. 4p applied to the differential 
equation (Ill.Sp does not actually produce a result modulo 16 (although this is what 
it would be designed to). It only produces the result modulo 8 given in Theorem [26] 
since, at the mod-16-level, the arising system of equations has no polynomial solutions. 
Nevertheless, by applying the enhanced method from Appendix [D] to this last system 
of equations, a solution modulo 16 can still be found, the result being displayed in our 
next theorem. 

Theorem 28. Let $(z) = J2n>o^'^" ■ Then we have 

J2sn+l{SL2{Z))z^ 
n>0 

= 8z^^ + 8z'^^ + 8z^^ + 8z^^ + 8/2 + gz^^ + 8z" + 8^^^ + 8z^^ + 8z^° + 8^^ + 8z^^ 
+ 8/^ + 8z^^ + 8z^^ + 8/2 + 8z^^ + 8z^° + 8z^^ + 8z^^ + 8/^ + 82^6 ^ 3^24 ^ g_^2i 
+ 12^20 + 12^1^ + 8z^^ + 8z^^ + 4:Z^^ + Az^^ + 12/ + Uz^ + 8/ + 12/ + 16z^ 
, ,0 10 12 3 14 9 4 12 4 3 6 

+ 12/ + 12/ + 8/ + 8+ — + — + — + — + — + — + — + — + — + —r 

z z^ z-^ z^ z^ z^ z^ z^ z^ z^^ 

+ ( 8;Z^^ + 8Z^2 ^ 3^71 ^ g^69 ^ g^68 ^ g^67 ^ g^66 ^ g^64 ^ g^63 ^ g^62 ^ g^61 

+ 8/° + 8z^^ + 8z^^ + 8/^ + 8/^ + 8/^ + 8/^ + 8/^ + 8/*^ + 8z^^ + 8/^ + 8/° 

+ 8^29 + 8;228 ^ g^27 ^ g^24 ^ g^23 ^ g^21 ^ g^l7 ^ g^l4 ^ g^l3 ^ g^l2 ^ g^9 ^ g^8 

7 -^ 4 ^ 16 6 12 8 6 4 4 8 14 \^, , 

+ 8/ + 8/ + 12z^ + 8/ + 4+ — + — + — + — + — + — + — + —r + —- $ ^ 

z z'^ z-^ z^ z^ z' z^ z^^ Z^^ J 

+ ( 8Z^2 ^ g^69 ^ g^68 ^ g^66 ^ g^65 ^ g^62 ^ g^61 ^ g^58 ^ g^57 ^ g^56 ^ g^55 



+ 8z^=^ + 8z^i + 8z^° + 8/^ + 8/« + 8z=^^ + 8;^=^^ + 8/^ + 8z'^ + 8z^'' + 8^24 + g^22 
+ 8z^^ + 8z^'^ + 82^^ + 82" + 8z^^ + 8/ + 8/ + 12/ + 8/ + 8 + - + -| + ^ 

2; Z^ Z'^ 
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4 10 12 12 8 2 \ 2/ 

+ 7 + ^ + 1^ + ;^ + ;^"^;^'^'' 

+ {8z'' + 8/^ + 8/^ + 8/^ + 8/^ + 8/° + 8/^ + 8/^ + 8/° + 8z'' + 8z'' 
+ 8z^' + 8z'^ + 8z^° + 8^3^ + 8^3=^ + 8z^^ + 8;^='^ + 8^30 + s^^s ^ §^25 ^ g^20 ^ g^i9 

+ 8z^^ + 8z^^ + 8^^^ + 8z^2 + 8^^° + 8^^ + 8z^ + 8^^ + 12^^ + 8^ + _+ + + 

z z^ z-^ z^ 

12 12 18 8 12 8 4 12 10 > o, , 

+ — + — + — + - + — + ^ + — + ^ + — I '^^^) 
z^ z° z' z^ z^ z^^ z^^ z^^ z^-^ 

+ ( 8z^2 + 8z''^ + 8^™ + 8z^'^ + 8/^ + 8/^ + 8/^ + 8/° + 8z^^ + 8z^^ + 8^^^ 

+ 8z^^ + 8z^^ + 8z^^ + 8^^^ + 8z^^ + 8z^^ + 8^^^ + 8z^° + 8^^^ + 8^^^ + 8z^^ + 8^^^ 
+ 8z^2 + 8z^° + 8^2^ + 8z^^ + 8;z^^ + 8;^^^ + 8z^^ + 8z^^ + 8z^^ + 8z^^ + 8;z^^ + 8z^ 

+ 12/ + 8z^ + 8/ + 8z^ + 12/ + 12/ + 6^2 + 8;z + 12 + - + — + — + — 

Z Z'^ z"^ z^ 

14 10 5 12 14 12 2 2 1 ,^^, , 

+ I 8/^ + 8z^^ + 8^*^^ + 8/2 + 8z^^ + 8z^^ + 8/^ + 8z^^ + 8z^° + 8^^^ 

+ 8/*^ + 8/^ + 8/3 + 8/° + 8/9 + 8/^ + 8/^ + 8/=^ + 8/« + 8/*^ 
+ 8^23 + 8^2° + 8;z^^ + 8;z^^ + 8^^^ + 8z^^ + 8/ + 8/ + 8/ + 8^^ 

4.9 4 8 14 12 14 4 4 6 \^,, , 

+ 8/ + 4/ + 8z + - + - + — + — + — + - + — + — U^ 2; 
2;^ ^r'^ 2;* z^ z' z^ z^^ z^-^ J 

modulo 16. (11.14) 

We did not attempt to push this analysis further to moduli 32, 64, etc., since the 
computational effort seemed immodest. With the (not very substantial) evidence of 
Theorems [2B1 and [2S] (but see Remark [52]) . we still expect the enhanced method to be 
successful for any given 2-power. 

Conjecture 29. Let $(z) = ^„>o^^"; ^^^ ^^^ 'j be a positive integer. Then the 
generating function X]n>o '^"+i('^-^2(^)) -z"; reduced modulo 2'^, can be expressed as a 
polynomial in $(z) with coefficients that are Laurent polynomials in z over the integers. 

12. Subgroup numbers for the lift r3(3) 

Continuing in the spirit of the previous section, we now consider the number of index- 
n-subgroups in the lift r3(3) (of the Hecke group i?(3) = PSL^iJ^^ modulo powers of 2. 
We shall see that, again, our method from Section |4] already fails for modulus 8. While 
this can again be overcome by, instead, designing the computation so that the target is 
modulus 16, the method then fails at the level of modulus 16. Moreover, for modulus 16, 
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even the enhancement of the method described in Appendix [D] fails (see Remark | 
This means that a new phenomenon, not covered by our Ansatz, arises in the behaviour 
of the subgroup numbers at the level of modulus 16. It would be of great interest to 
find an explicit description of the hidden scheme behind the mod-16 behaviour of the 
number of subgroups of index n in r3(3), and, more generally, of the behaviour modulo 
any power of 2. 

We take F = r3(3) in flll.ip and combine the resulting formula with fll0.6p . the lat- 
ter giving an explicit formula for the homomorphism numbers hn '■=\ Hom(r3(3), 5'„)|. 
Using the Guessing package [IS], we found a recurrence of order 42 for the sequence 
{hn/n\)n>o, with coefficients that are polynomials in n over Z. The validity of the recur- 
rence was verified by computing a certificate using Koutschan's Mathematica package 
HolonomicFunctions [19j|j However, again, this recurrence is not suitable for our pur- 
pose, for which we require a recurrence with coefficients that are polynomials in n over 
Z, and with leading coefficient n. By an indeterminate Ansatz, we computed a can- 
didate for a recurrence of the desired form of order 60, with polynomial coefficients of 
degree at most lOo To be precise, it is the uniquely determined recurrence of the form 

2^[2^Kk.t)n ) -— = 0, n>60, 

A;=0 ^ j=0 ^ ^ '' 

where 

6(0, 0) = 6(0, 2) = 6(0, 3) = 6(0, 4) = 6(0, 5) 

= 6(0, 6) = 6(0, 7) = 6(0, 8) = 6(0, 9) = 6(0, 10) = 0, 
6(0,1) = 1, 
6(60, 8) = 9649124343496238177846526221678676069879148435557456840677\ 
68567400990643919180258204664996863270960793634431477\ 
96875828563496094243333614632539311543926582958877938\ 
09887854513738722474642524334737161421912431106592005\ 
22984304410147101964876864298627928130880022459406799\ 
539461032349694733915947489297372243661012, 
6(60,10) =6(60,9) = 6(60,4) 

= 6(59, 10) = 6(59, 9) = 6(59, 8) = 6(59, 7) = 6(59, 6) = 6(59, 5) 

= 6(59, 4) = 6(59, 3) = 6(59, 2) = 6(59, 1) = 6(59, 0) 

= 6(58, 10) = 6(58, 9) = 6(58, 8) = 6(58, 7) = 6(58, 6) 

= 6(58,5) = 6(58,4) 

= 6(57, 10) = 6(57, 9) = 6(57, 8) = 6(57, 7) = 6(57, 6) = 6(57, 5) = 6(57, 4) 



The computation took about one week, producing a certificate of 28 megabytes. The coefficients 
of this recurrence are polynomials in n of degree up to 105. In this case, we were not able to find a 
recurrence with leading coefficient 1. (It may still exist.) The best that we found in this direction was 
a recurrence with leading coefficient a 1828-digit number. Again, although we did not try to prove it, 
it is likely that the recurrence of order 42 is the recurrence of minimal order. 

Again, we used the function LinSolveQ of the Guessing package fTB' in order to solve the arising 
system of linear equations. 
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= 6(56, 10) = 6(56, 9) = 6(56, 8) = 6(56, 7) = 6(56, 6) = 6(56, 5) = 6(56, 4) 
= 6(55, 10) = 6(55, 9) = 6(55, 8) = 6(55, 7) = 6(55, 6) = 6(55, 5) = 6(55, 4) 
= 6(54, 10) = 6(54, 9) = 6(54, 8) = 6(54, 7) = 6(54, 6) = 6(54, 5) = 6(54, 4) 
= 6(53, 10) = 6(53, 9) = 6(53, 8) = 6(53, 7) = 6(53, 6) = 6(53, 5) = 6(53, 4) 
= 6(52,10) = 6(52,9) =6(52,8) 
= 6(50,7) = 6(50,6) 

= 6(49, 10) = 6(49, 9) = 6(49, 8) = 6(49, 7) = 6(49, 6) = 6(49, 5) = 6(49, 4) 
= 6(1,7)=0. 

Subsequently, we checked that this recurrence is a left-multiple of the certified recurrence 
of order 42, thereby establishing validity of this candidate recurrence of order 60. This 
last recurrence was then converted into a linear differential equation with polynomial 
coefficients for the series 

H{z) := E^"^ = E |Hom(r3(3),S„)|^. 

n=0 ■ n=0 

Finally, this last mentioned differential equation can be translated into a Riccati-type 
differential equation for the generating function 

Siz) = J2sn+i{Tsi3))z^ (12.1) 

n>0 

for the subgroup numbers of r3(3) in the same way as we obtained (I11.5P in the previous 
section. It turns out that this differential equation has integral coefficients, so that it is 
amenable to our method from Section |H The differential equation cannot be displayed 
here since this would require about 100 pagesO Its reduction modulo 16 is written out 
in (112. 3p . By applying our method from Section H] with the minimal polynomial for the 
modulus 16 (!) in place of the polynomial on the left-hand side of (14. 4p to (112.31) . we 
obtain the following theorem. It refines the parity result [201 E<1- (6.3) with \H\ = 1, 
g = 3, m = 3]. 

Theorem 30. Let ^{z) = ^„>o-2^"- Then we have 

J2sn+l{T3{3))z^ 
n>0 

= 4/2 + 4^53 ^ 4^44 ^ 4^35 ^ g^26 ^ 4^20 ^ 4^14 ^ 4^12 ^ 4^11 ^ 4^10 



2 



'~7 T O fi 

+ Az'' + Az'^ + 6z'' + Az^ + Az^ + Az + Q + — + — + — + — 

z^ z-^ z° z° 

/^. 3 . 2 4 6 6 6 2\ , , 

+ [ Az^ + Az^ + - + — + — + — + — ] $ z 
\ z z^ z"^ z^ z^ ) ^ ' 

+ ( 42;*^ + 4z^ + 42;^ + 6^2 + 4 + - + — + — + — + — + — + ^ 1 ^''{z) 

z z^ z-^ z^ z^ z^ z' 
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The integers appearing as coefficients have up to 320 digits. 
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+ (^4z' + ^ + 1 + 1 + 1 + 1 + 1) *'(^) modulo 8. (12.2) 

Proof. The Riccati-type differential equation for S{z) (as defined in (112. ip ) modulo 16 
iS 

qo{z) + qi{z)S{z) + q,{z)S{z)S\z) + q,{z) S (z) S' {zf + q,{z) S (z) S' {zf 

+ q5iz)Siz)S'izf + qe{z)Siz)S"iz) + q7iz)Siz)S"izf + qsiz)Siz)S"'iz) 

+ q,{z)S{z)S"'{zf + q,o{z)S{z)S""'{z) + q^r{z) S (z) S' (z) S" {z) 

+ qu{z)Siz)S'{z)S"'{z) + qi3{z)S{z)S'{zfS"'iz) + q,,{z)S{zf + q^,{z) S {zf S' {z) 

+ q,,{z)S{zfS'{zf + q,-j{z)S{zfS\zf + q,,{z)S{zfS\zt + q,,{z)S{zfS"{z) 

+ q^,{z)S{zfS"'{z) + q,,{z)S{zfS"'{zf + q^,{z)S{zfS""{z) + q,,{z)S{zfS'{z)S"{z) 

+ q,,{z)S{zfS'{zfS"{z) + q,,{z)S{zfS'{z)S"'{z) + q,,{z)S{zfS'{zf S"'{z) 

+ q27{z)S{zf + q2^{z)S{zfS\z) + g29(^)5(^)'5'(^)' + q,,{z)S{zfS\zf 

+ g3i(^)5(^)'5"(^) + g32(^)5(^)'5'"(^) + q,,{z)S{zfS'{z)S"\z) + g34(^)5(^)' 

+ q;,{z)S{zfS'{z) + g36(z)5(^)^5'(^)' + g37(^)5(^)'5'(^)' + q,,{z)S{zYS"{z) 

+ g39(^)5(^)'5'"(z) + q,o{z)S{zYS\z)S'\z) + g4i(^)5(^)'5'(^)5'"(^) + g42(^)5(z)^ 

+ q,,{z)S{zfS\z) + g44(^)5(^)'5'(^)' + q,,{z)S{zfS"{z) + g46(^)5(^)'5'"(;2) 

+ g47(^)5(z)6 + q,s{z)S{zfS'{z) + g49(^)5(z)65'(z)' + q^,{z)S{zfS"{z) 

+ g5i(;2)5(;2)^5'"(^) + g52(^)5(^)' + q,:i{z)S{zyS\z) + g54(^)5(z)*^ + q,^{z)S{zfS\z) 

+ 556(^)5(^)9 ^ g57(^)5(z)i° + q,^{z)S\z) + g59(^)5'(z)S"(^) + q^^{z)S\z)S"\z) 

+ g6i(^)5'(;.)S'"(^)' + g62(^)^'(^)5""(^) + q,,{z)S\zf + g64(^)5'(^)'5"(^) 

+ q,,{z)S\zfS"\z) + g66(^)^'(^)' + g67(^)5'(^)'5"(^) + q,s{z)S\zfS"\z) 

+ g69(^)^'(^)' + q7o{z)S\zf + gn(^)5"(z) + g72(^)5"(^)5'"(^) + q7z{z)S"{zf 

+ g74(^)^'"(^) + q7,{z)S"\zf + g76(^)^""(^) + g77(^)5'""(^) = 

modulo 16, (12.3) 
with coefficients 5^(2;), j = 0, 1, ... , 77 as displayed in Appendix O 
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We display the differential equation modulo 16 in order to prepare for Remark [ 
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The differential equation f ll2.3p has a unique solution since comparison of coefficients 
of z^ fixes the initial values, and yields a recurrence for the sequence (s„(r3(3)))^^j^ 
which computes s„+i(r3(3)) from terms involving only Sj(r3(3)) with i < n. 

Now we apply the method from Section H] with the polynomial in flll.6p in place of 
the polynomial on the left-hand side of ( 14. 4 p to the differential equation (112.31) . This 
yields the above result by means of a straightforward computer calculationo D 

Also here, if we want to know criteria in terms of n when a subgroup number s„(r3(3)) 
is congruent to a particular value modulo 8, then we must first apply the algorithm from 
Section [3] to the right-hand side of (112. 2p . This leads to the identity 



n>0 ^ ^ ^ ^ 



2 



4 4 6 4 4 6 



4=' + ;i + 73 + ii + ;s + pj + irl^»(- 



I A R .A .•^^9 . 4 6 6 5 6 2 5.^^,, 

+ [ Az^ + Az^ + Az^ + Qz^ + A+ - + — + — + — + — + — + — ] Hi{z) 
* z z^ z^ z^ z^ z^ z^ ' 



.62 , .53 , .^44 , .35 , « 26 , .20 



+ 2z^ + 2z^ + Az'^ + 2 + - + ^ + ^ + ^ + ^ + ^ modulo 8, (12.4) 



+ Az""' + Az°' + Az"" + 4^^° + Qz''' + Az'" + Az"" + Az'' + Az'' + Az' 

6 14 6 13 
z z^ z^ z^ z^ z^ 

from which we can extract the following explicit description of the behaviour of the 

subgroup numbers of r3(3) modulo 8. 

Theorem 31. The subgroup numbers s„(r3(3)) obey the following congruences modulo 
8: 

(i) s„(r3(3)) = 1 (mod 8) if, and only if, n = 1,2, 10, or if n is of the form 2°" — 3 
for some cr > 4; 

(ii) Sn(r3(3)) = 2 (mod 8) if, and only if, n = 7,9, 17, 18,27,42, or if n is of one 
of the forms 

3 • 2"" - 3, 3 ■ 2"" - 6, for some a > 5; 

(iii) Sn(r3(3)) =4 (mod 8) if and only if, n = 3,12,22,23,36,38,39,43,46,49,50, 
51, 53, 54, 63, or if n is of one of the forms 

2" + ^, 2" + !, 2" \\\, 2" \\A, 2'" + 17, 2^^ + 18, 2^^ + 19, 2^^ + 21, 

for some cr > 6, (12.5) 

2^^ + 2^ - 2, 2^" + 2^ + 1, 2^" + 2^ + 2, 2^" + 2^ + 3, 2'" + 2^ + 5, 2^" + 2^ + 10, 

2^^ + 2^ + 13, for some a, r with a > Q and 5 < t < a - 1, (12.6) 

2^^ + 2^ + 2^ - 6, 2^" + 2^ + 2*^ - 3, 

for some a, t, v with a >1 , Q < v < a — 1, and 5 < r < z/ — 1; (12.7) 



The calculation being straightforward, it nevertheless required a machine with substantial amount 
of memory (we had available 32 gigabytes of memory, of which almost 50% were used). 
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(iv) Sji(r3(3)) = 5 (mod 8) if, and only if, n = 5, or if n is of the form 2*^ — 6 for 
some cr > 5; 

(v) Sn(r3(3)) = 6 (mod 8) if, and only if, n = 6, 11, or if n is of one of the forms 
2<^-2, 2'' + 3, 2'" + 4, for some a > 4, (12.8) 

2'" + l, 2'' + 2, 2'' + 13, for some a > 5, (12.9) 

2'' + 10, for some a >Q, (12.10) 

2^" + 2^ - 6, 2'" + 2^ - 3, 

for some a, t with cr > 7 and 5 < r < a — 2; (12.11) 

(vi) in the cases not covered by items (i)-(v), s„,(r3(3)) is divisible by 8; in partic- 
ular, s„(r3(3)) ^ 3, 7 (mod 8) for all n. 

Remark 32. In the application of the method from Section H] in the proof of Theorem [311 
when we arrive at the mod-8-level, we obtain 



Y^ 5^+1(13(3)) Z"" = 4/2 + 4z^3 + 4^44 ^ 4^35 ^ g^26 ^ 4^2 



+ 4/^ + 4/2 ^4^1^ +4/ + 4^*^ + 4/ + - + 4 + 4 + 4 



z z"^ z^ z^ 



4 2 4 2 

Z^ Z'' Z'' J \ z z* z^ z^ 



+ (4^' + 4 + ^ + 3 + iF)*(^)+(4-' + 7 + T^ + 7I + -)*'(^) 



/ . 2 4 4 2 4 4 2 > 3, , 

+ (4^^ + ;i + ;^ + ^ + ;F + ;a + ^)^^W 

+ ( 4/ + 4^4 + 4/ + 6^2 + 12 + - + ^ + 4 + 4 + ^ + 4 + ^ 1 '^"('2) 

z z^ z"^ z^ z° z° z' 



, . 4 6 4 6 

+ ( 42;2 + - + - + - 

z^ z^ z^ z 



\\^^{z) modulo 8. (12.12) 



However, the system of equations for the next level, the mod-16-level, has no polynomial 
solutionso Even the enhancement of our method described in Appendix[D]fails. (There 
are actually several problems arising. It turns out that, due to the reduction modulo 2, 
the variables hi{z^ expressed in (ID.lOp do not solve the system (ID. 90 unless one puts 
further restrictions on a^°^{^z\ a^^\z\ . . . , d^°^[^z\ d^^\z\ But even if we ignore that and 
continue to follow the procedure described in Appendix [Dl then a contradiction arises 
at a later point: one of the factors of the polynomial Pi^z^ turns out to be (1 + 2;)^°, and 
the congruence (ID. lip with P- ^(2;) = (1 + z)^" has no solution.) This proves that it is 
impossible to find a polynomial in $(z) with coefficients that are Laurent polynomials in 
z over the integers which agrees with the generating function for the subgroup numbers 
of r3(3) modulo 16. 



The corresponding computation took almost 5 hours, using 94% of the 32 gigabytes of memory 
of the machine on which the computation was performed. 
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13. A VARIATION I: FREE SUBGROUP NUMBERS FOR LIFTS OF HeCKE GROUPS 

In this section, we consider the functional equation 

zf\z)-f{z) + 1 = 0, (13.1) 

which generahses the functional equation (15. ip for the generating function of Catalan 
numbers. It is easy to see that this equation has a unique formal power series solution. 
The coefficients of this uniquely determined series can be calculated explicitly by means 
of the Lagrange inversion formula, the result being 

{zlfiz) = -(^\), (13.2) 

n\n — 1/ 

but this will not be relevant hereo Again, the numbers in (113. 2p are special instances 
of numbers that are now commonly known as Fufi-Catalan numbers (cf. the paragraph 
containing (16. ip ). 

Our aim is to determine the coefficients of f{z) modulo powers of 2. Our solution 
of this problem is that, again, the series f{z) can be expressed as a polynomial in a 
"basic" series. Here, this basic series is 

$,(^) = ^^2-/{2'^-i). (13.3) 

n>0 

It will turn out (see Corollary [MP that an adaptation of the proof of the theorem 
below will allow us to treat as well the behaviour, modulo powers of 2, of free subgroup 
numbers of lifts of Hecke groups S){q), with q a Fermat prime. 

The theorem below, in a certain sense, extends Theorem [131 It does not, however, 
reduce to it for /i = 1, due to the choice that, in the proof below, the reductions in our 
algorithm are based on the polynomial relation (113. 7p for the basic series $h(-2), which, 
for /i = 1, is "weaker" than the relation (14. 4p which is used in the proof of Theorem [T3l 

Theorem 33. For a positive integer h, let $/i,(-2) = X]n>0'^^" ^''^ "^^ ^'"'^ ^^^ a be a 
further positive integer. Then the unique solution f{z) to (113. ip . reduced modulo 2^" , 
can be expressed as a polynomial in ^h{z) of degree at most 2^"+^^^ — 1 with coefficients 
that are Laurent polynomials in z^'^^ ^^' over the integers. 

Proof. For ease of notation, we replace z by z^"^^ in (113. ip . thereby obtaining the 
equation 

z^'-^f'\z)-f(z) + 1 = 0, (13.4) 

with f{z) = f{z'^ ~^). We now have to prove that, modulo 2^" , the series f{z) can be 
expressed as a polynomial in 

oo 

$,(.) = 5^ z^"^ (13.5) 

n=0 

of degree at most 2*^""*"^^^ — 1 with coefficients that are Laurent polynomials in z. 
It is readily verified that 

l>f {z) + ^h{z) + z = modulo 2, (13.6) 



i^See Footnote [31 
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whence 

(^^t{z) + ^h{z) + zy =0 modulo 22"'. (13.7) 

We modify our Ansatz (14.21) to 

2(Q+l)h_i 

f{z) = J2 <^Whiz) modulo 2^"\ (13.8) 

i=0 

where the ai{z)^s are (at this point) undetermined Laurent polynomials in z. 

Next, we gradually find approximations ai^p{z) to ai{z) such that (113.41) holds modulo 
2^, for f3 = 1,2, ... , 2"^. To start the procedure, we consider the differential equation 
in^ modulo 2, with 

f{z)= Yl «mW^U^) modulo 2. (13.9) 

We substitute the Ansatz (113. 9p in (113. 4p . reduce high powers of ^h{z) by using Re- 
lation (113. 7p . reduce the resulting expression modulo 2, thereby taking advantage of 
the elementary fact that ^'h{z) = 1 modulo 2, and we finally see that the left-hand 
side of (I13.4P becomes a polynomial in ^h{z) of degree at most 2'^"+^^^ — 1 with coeffi- 
cients that are Laurent polynomials in z. Now we compare coefficients of powers $^(;z), 
k = 0,1, . . . , 2^"+^)'^ — 1. This yields a system of 2^°^^^^^ equations (modulo 2) for the 
unknown Laurent polynomials ai^i{z), i = 0,1, . . . , 2(°+^)'* — 1. Since we have already 
done similar computations several times before, we content ourselves with stating the 
result: all Laurent polynomials ai^i{z) must be zero, except for ao,i(-2) and a2cyhi{z), 
which are given by 

a-l 

«o,i(^) = 2^2 

fc=0 

a2ahx{z) = z~^. (13.10) 

After we have completed the "base step," we now proceed with the iterative steps 
described in Section HI Our Ansatz here (replacing the corresponding one in (I4.6p - (l4.8p ) 
is 

]{z)= Yl «^/3+iW^U-2) modulo 2'^+\ (13.11) 

1=0 

with 

ai,^+i{z):=a,A^)+2''bi,p+,{z), z = 0, 1, . . . , 2("+i)'^ - 1, (13.12) 

where the coefficients ai^jj^z) are supposed to provide a solution 

2(c,+i)h_i 
i=0 

to (113. 4p modulo 2^. This Ansatz, substituted in (I13.4p . produces the congruence 

z'''-'fJ{z)-fp{z)+2^ Y 6,^+i(^)l>U^) + l = modulo 2/^+1. (13.13) 

i=0 
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By our assumption on fi3{z), we may divide by 2^. Comparison of powers of ^h{z) then 
yields a system of congruences of the form 

bi,(s+iiz) + Poli(z) = modulo 2, i = 0, 1, . . . , 2^'^+^'^^ - 1, (13.14) 

where Polj(2;), i = 0, 1, . . . ^2^""*"^-"* — 1, are certain Laurent polynomials with integer 
coefficients. This system being trivially uniquely solvable, we have proved that, for an 
arbitrary positive integer a, the modified algorithm that we have presented here will 
produce a solution f2<^h{z) to (113. 4p modulo 2^° which is a polynomial in ^h{z) with 
coefficients that are Laurent polynomials in z. D 

It has been shown in [27] that the parity pattern of free subgroup numbers in Hecke 
groups Sj{q), q a Fermat prime, coincides with the parity pattern of (special) FuB- 
Catalan numbers. More precisely, let f[^ denote the number of free subgroups of index 
2qX in the Hecke group S^{q). (For indices not divisible by 2g, no free subgroups exist 
in Sj{q).) Then (see [271 Eq. (37)]) 



/w^ = - r \ ; modulo 2. 



" ~ X[ A-1 

The reader should keep in mind that g — 1 is a 2-power. Theorem |33] says that the 
generating function for the Fufi-Catalan numbers (113.21) . when reduced modulo a given 
power of 2, can be expressed as a polynomial in ^h{z)- We are now going to show that 
the same is true for the generating function for free subgroup numbers in the Hecke 
group i5(g), although the equation it satisfies is different from the functional equation 
(113. ip for the generating function of Fufi-Catalan numbers. In the corollary below, we 
present actually a more general result: even the generating function for free subgroup 
numbers of the lift Tmiq), when reduced modulo a given power of 2, can be expressed 
as a polynomial in ^h{z) in the case where g is a Fermat prime. In a certain sense, this 
extends Theorem [ini although it does not reduce to it for h = 1. Again, the reason 
lies in the choice that, in the proof below, the reductions in our algorithm are based on 
the polynomial relation (I13.7P for the basic series ^h{z), which, for h = 1, is "weaker" 
than the relation (14.41) which is used in the proof of Theorem [121 On the other hand, 
the corollary does largely extend the parity result [271 Cor. A']. 

Corollary 34. Let g = 2^ + 1 be a Fermat prime, and let 7 he some positive in- 
teger. Then, for every positive integer m, the generating function Fm{q', z) = 1 + 
Ylx>i fx i''^)^'^ 0/ /ree subgroup numbers ofTmiq), when reduced modulo 2^^, can be 
expressed as a polynomial in $2/ (-2) 0/ degree at most 2^ 7 — 1 with coefficients that are 
Laurent polynomials in z^^^'^~'^\ where the series ^h{z) is defined as in (113. 3p . 

Proof. In view of Proposition [181 the assertion is trivially true for even m, the poly- 
nomial in $2/(^) being a polynomial of degree zero in this case. We may thus assume 
from now on that m is odd. 

Equation (17. 3p provides a Riccati-type differential equation for Fm{q; z) = 1 + 
zGm{(l]z). Moreover, this equation, considered modulo 2, is the same for every odd 
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m. Namely, we have 

„-i „-i „, „ -,n V/^1' • • • '/^J^/ 



MlH \-IJ.iy=fl 



(m,-i) 



l[(^l{Fmiq;z)-l)y' modulo 2 



X 

Moreover, it is shown in f27\ Prop. 2] that, modulo 2, this differential equation reduces 
to 

zF^-^q; z) - F„,{q; z) + 1 = modulo 2. (13.15) 

The latter statement means that reduction of coefficients modulo 2 and usage of the 
simple fact that 

F^{q; z) = modulo 2 (13.16) 

leads from the original differential equation 07. 3p for Fm{q; z) = 1 + Gm{q',z) to the 
congruence (113. ISp . With q — I being a power of 2 by assumption, we observe that, 
disregarding the restriction to modulus 2, Equation (113. 15p is the special case of (113. ip 
where h = 2^ . In particular, if, for the moment, we assume that 7 = 2"^ , for some 
positive integer a, then we see that the base step of the Ansatz outlined (and applied) 
in the proof of Theorem [33] (with h = 2^) can be successfully performed here: it would 
yield exactly the same result as there, namely (113.90 with the Laurent polynomials 
ai^i{z) being given in the paragraph containing (113. lOp . 

However, also the subsequent iterative steps would just work in the same way as 
in the preceding proof! Indeed, transform the Riccati-type differential equation (17. 3p 
for Fm{q', z) by the substitution z \-^ z'^~^ (in analogy with the substitution leading to 
(113. 4p ). This yields a Riccati-type differential equation for Fm{q', z'^~'^) . A fine point 
to be observed here is that, in this equation, the coefficients will not necessarily be 
integral; due to the substitution rule for differentials, denominators that are powers of 
{q — 2) may occur. As in the proof of Theorem [331 we now continue with the Ansatz 

F„(g;z'?-2)= Yl «.,/3+i(^)'^2/(^) modulo 2^+\ (13.17) 

with 

ai,p+iiz):=a,,p{z)+2''bi,p+iiz), t = 0,1, . . . ,2^''+'^^' - 1 (13.18) 

(which is analogous to (I13.1ip - (ll3.12p ). where the coefficients ai^i3{z) are supposed to 
provide a solution 

to the differential equation for Fm{q; z'^^^) modulo 2^. The fact that reduction modulo 2 
and usage of (113.160 leads from the original differential equation for Fm(g; z) to (113. ISp 
implies that substitution of the Ansatz (113.170 - 013.180 in the differential equation for 
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Fm{q', z'^^'^) yields an equation completely analogous to (113. 13p . namely 

m,/3 ' 



■^^-'^'^-'{q;z)-F„,AT,z) + 2^ J2 h^+ii^Whi^) 



i=0 

+ 1 + T{z, Fm,p{q; z)) = Q modulo 2^^+^ 

Here, T{z, Fm^p^q; z)) consists only of terms that may depend on Fm^/siq', z) but do not 
depend on the 6j,g+i(2;)'s. The rest of the procedure is then as in the preceding proof: 
we divide by 2^, compare powers of $/i(z), and obtain a system of congruences of the 
form fll3.14l) . which is trivially solvable. The powers of (g — 2) that may appear in the 
denominators of the coefficients in the polynomials involved here are disposed of by 
interpreting them appropriately as elements of 'L/2'^'L. 

f 

Finally, if we are able to express Fm{q', z) as a polynomial in $2/ (-2) modulo 2^ = 2^" 
for all a, then the same assertion must hold for every 7. D 

In order to illustrate the algorithm described in the last proof, let us consider the 
case of the Hecke group i?(5), that is, the case of Corollary |3H where / = 1. The 

Riccati-type differential equation for the series Gm(^) := Gm{5;z) = ^'^Lo fx+ii''^)^^ 
that one obtains from (17. 3p in this special case reads 

Graiz) = 189m^ + AGOOm^zGmiz) + U30m^ z^ Gl^{z) + 80mz^Gl^{z) 

+ z^Gm{z)^ + U300mh'^G'^{z) + 2400m^z^Gn,{z)G'^{z) + QOmz^Gl,{z)G'^{z) 

+ 300mh%G'^{z)Y + 8000mh^G'^{z) + AOOmh^GUz)Gl{z) + 1000mh^GZ{z). 

Since Gm{z) = Gmi^; z) = ^(^^(5; z) — Ij, we obtain the differential equation 

1 + (189m^ - 300m3 + UOm"^ - 20m + l)z 

+ ((300m^ - 260m^ + 60m - A)z - 1)F™(5; z) + {130m^ - 60m + 6)zF^{5; z) 

+ (20m - 4)2F^(5; z) + zF^{h] z) + (4300m3 - lOOOm^ + 60m)2^F4(5; z) 

+ (lOOOm^ - l20m)z^F^{h] z)F'^{5; z) + mmz^F^{^; zfF'J^^; z) + 300m^ z^ F'J^^; zf 

+ (SOOOm^ - 400m2)2^F^(5; z) + 400m2;z^F„(5; ^)F^(5; z) + 1000mVF^'(5; z) = 

(13.19) 

for -Fm(5; z\ We have implemented the algorithm described in the proof of Theorem IMl 
for this differential equation. For the modulus 16, it produces the following result. (It 
is independent of m because of the high divisibility of the coefficients in the differential 
equation (113. 19p by powers of 2. The parameter vn will show up for 2-powers higher 
than 2^ = 16.) 

Theorem 35. Let $2(-2) = X]n>o'^^"''^' ^^ before. Then, for all positive odd integers m, 
the generating function Fm{5; z) = 1 + J2\>i fx im)z^ for the free subgroup numbers 
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of Tm{5) satisfies 
F„(5; z)=4z + l + 12z2/^$2(;z) + lOz^^^^Kz) + U^z) + (4^^/^ + 7z-^/^) <l>^(z) 
+ 4z^/^^l{z) + 4<l>l{z) + 12z-^/^<!>l{z) + 8z^/^<l>l{z) + 4$^(z) 

+ 2z-^/^(l>l°{z) + 12(l>l^{z) + 12z-^/^(l>l^{z) modulo 16. (13.20) 

Clearly, coefficient extraction from powers of $2(-2) (and, more generally, from powers 
of ^h{z)) can be accomplished by appropriately adapting the results in Section [3l 

14. A VARIATION II: SUBGROUP NUMBERS FOR HeCKE GROUPS 

In Section [9l we proved that the generating function for the subgroup numbers of 
the inhomogeneous modular group PSL2{'L) = i)(3), when reduced modulo a power 
of 2, can always be expressed as a polynomial in the basic series ^{z) with coefficients 
that are Laurent polynomials in z. Here, we discuss possible extensions of this result to 
Hecke groups S){q), where g is a Fermat prime. Again, we have to modify the original 
method from Section |4] by using the series ^h{z) defined in ( 113. Sp (for suitable h) instead 
of ^{z). We conjecture (see Conjecture [3H1) that this variation of our method will be 
successful for arbitrary Fermat primes q. If g = 5, we are actually able to demonstrate 
this conjecture, thereby largely refining the g = 5 case of [271 Theorem B]. 

Theorem 36. With notation from the previous section, let ^2{z) = '^n>o^^"^^^ '^^^ 
let a be a positive integer. Then the generating function S{z) = S^(^5){z) {see the 
first paragraph of Section^ for the definition), reduced modulo 2*^", can be expressed 
as a polynomial in $2(-2) of degree at most 4"+-*^ — 1 with coefficients that are Laurent 
polynomials in z^l"^ over the integers. 

Proof. Let 

h{n) = —hc2{n)hc,{n). 

Using the routine to compute recurrences for the Hadamard product of recursive se- 
quences, implemented in gfun [31] and GeneratingFunctions |23j (cf. [32l Theo- 
rem 6.4.12] for the theoretical background), one obtains the recurrence 

n(16 - 72n + 17471^ - 155n^ + 65n^ - 13n^ + n^)h{n) 

- (184 - 620n + 854n^ - 555n^ + 177n^ - 25n^ + n^)h{n - 1) 

- (856 - 1636n + 125071^ - 479^^ + lOln^ - 13n^ + n^)h{n - 2) 

- 4(n - 6)(n - ?,f{-7 + ?,n)h{n - 3) + 8(n -7){n- 4)(3n - 7)h{n - 4) 

- (1136 - 856n + 129271^ - 2930n^ + 3115n^ - 1718n^ + 516n^ - 80n^ + 5n®)/i(n - 5) 

- 2(1856 - 5376n + 6828^^ - 4868n^ + 2174n^ - 651n^ + 133n^ - 17n^ + n^)h{n - 6) 
-4(r2-6)(r2-5)(r2-3)2(r2-2)(3r2-7)/i(n-7)+8(n-7)(n-6)(n-4)(n-3)(3n-7)/i(n-8) 

- 16(n - 8)(n - 7)(n - 5)(-7 + ?>n)h{n - 9) 
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-{n-9){n- 8){n -6){n- 3)(16 + 12n - 16n^ - 5n^ + 15n^ - 7n^ + n^)h{n - 10) = 

(14.1) 

for the sequence (^h{n))^^^. Since the leading coefficient (i.e., the coefficient of h{n)) is 
not n, this recurrence is not suitable for being translated into a Riccati-type differential 
equation with integral coefficients via (lll.ip . to which we can apply our method from 
Section HJ Using Euclidean division of difference operators, one can see that we also 
have 



n 



h{n) - h{n - 1) - h{n - 2) - {5n^ - lln - 44) h{n - 5) - 2{n - 4)(n - 2)h{n - 6) 
+ 12/i(n - 7) - 4/i(r2 - 8) - [n^ - 20n^ + 95^^ + 260n - 2000) h{n - 10) 
+ 4(9n - 85)/i(n - 11) - 8 [n^ - 19n + 89) h{n - 12) 
+ 4(n - 14)(n - 13)(n - ll)(r2 - l)h{n - 15) 

- 4(n - 15)(n - 14)(n - 12)(n - 9)/i(n - 16) = 0. (14.2) 

If we now apply the procedure of converting such a recurrence for homomorphism num- 
bers (divided by n!) into a Riccati-type differential equation for the generating function 
of the corresponding subgroup numbers as explained in the paragraph containing 011.2p . 
then we obtain the differential equation 

224z^^ - 256^^^ + 40z" - 56^^° + lOO/ + 4z^ - Vlz^ + \<oz^ + 1<6z^ + z + 1 

+ (736^^^ - 824;z^^ + 48z^2 - 36^^^ + 276z^° + 14^^ + 44z^ - l) S{z) 

+ (448^^^ - 488z^^ + 8^^^ + 162^" + Iz'' + 5z^) S'^{z) 

+ [mz^^ - Mz^' + 262^2^ S^{z) + (42^9 - 4^1^ + z^'') S\z) 

+ (448z^^ - 488^^^ + 8;z^^ + 162;z" + 2;z^ + 5z^) S'{z) 

+ (I2z^9 - 12;2^^ + 'iz^'') {S'f{z) + (240^^^ - 252^1^ + 78^^') S{z)S\z) 

+ (24^19 - 24^1^ + 62^3) 52(;2)5'(2) + (80^1^ - 84^1^ + 26^^^) ^"(2) 

-F (16;2^9 - 16^^^ + 4^1=^) S{z)S"{z) + (4^19 - 4^i« + z^"") S"\z) = 0. (14.3) 

For convenience, we replace z by z^ in (114. 3p . Writing S{z) = S{z^), the above differ- 
ential equation translates into 

224z^^ - 256^^2 + 40z^^ - 56z^^ + lOOz^'^ + \z^^ - I2z^^ + IQz^^ + 26^^^ + z^ + 1 
+ (736^^^ - 824/^ + 48;z=^^ - 36^^^ + 276^^° + 14^^^ + 44z^^ - l) S{z) 
+ (448^^^ - 488;z^^ + 8^^^ + 162^3^ + 2z^^ + 5z^^) S'^{z) 
-I- (80/^ - 84^51 ^ 26^3'') S\z) + {Az'''' - Az''^ -|- z^"") S\z) 

- —z^^ (96^^^ - 3688z^^ + 3928z^° - 72^^! + 2Az^^ - 1312^^^ - 18^^ - 45) S'{z) 
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+ 1^35 (4^18 _ 4^15 ^ ^) (^/)2(^) _^ ^gQ^52 _ g4^49 ^ ^Qz'"') S{z)S'{z) 
o 

+ (SZ^S - 8Z'' + 2Z^') S'{Z)S'{Z) + ^^32 (4^21 ^ ^4^18 _ ^g^lS + ^3 ^ g) ^.(^) 

9 

+ 4=^'' (4^^' - 4^15 + 1) 5'"(2) = 0. (14.4) 

^ ( 

We have to prove that, modulo 2^", the series S{z) can be expressed as a polynomial in 
$2(2;) (as defined in (113. 5p ) of degree at most 4""^^ — 1 with coefficients that are Laurent 
polynomials in z. 
We make the Ansatz 

S{z) = J2 «i(^)*2(^) modulo 2^^ (14.5) 

i=0 

where the aj(z)'s are (at this point) undetermined Laurent polynomials in z. 

Next we gradually find approximations 04^/3(2) to ai{z) such that (114. 4p holds modulo 
2^, for /3 = 1,2, .. . ,4". To start the procedure, we consider the differential equation 
(ITOl) modulo 2, with 

S{z) = Y^ ai^i{z)^{z) modulo 2. (14.6) 

i=0 

We substitute the Ansatz (114.61) in (114.41) . reduce high powers of $2 (-2) by using the 
relation (113. 7p with h = 2, and reduce the resulting expression modulo 2, thereby taking 
advantage of the elementary fact that $2 (-2) — 1 modulo 2. The powers of 3 that appear 
in the denominators of the coefficients in the polynomials involved here are disposed 
of by interpreting them appropriately as elements of Z/2^'^Z. We finally see that the 
left-hand side of (114.41) becomes a polynomial in $2(2;) of degree at most 4°+^ — 1 with 
coefficients that are Laurent polynomials in z. Now we compare coefficients of powers 
$2 (-2) for k = 0,1, . . . , 4"+^ — 1. This yields a system of 4""*"^ equations (modulo 2) for 
the unknown Laurent polynomials ai^i{z), i = 0,1, . . . , 4""*"^ — 1. Since we have already 
done similar computations several times before, we content ourselves with stating the 
result: all Laurent polynomials ai^i{z) must be zero, except for 

a,^,{z) = Z-' + z-'' J2 ^'' + '~' E ''''' ' 
fc=i fc=i 

040 1(2;) = z , 

02.4^1(2:) = z-\ (14.7) 

After we have completed the "base step," we now proceed with the iterative steps 
described in Section |H Our Ansatz here (replacing the corresponding one in (I4.6l) - (l4.8p ) 
is 

S{z) = J2 «i,/3+i('2)^2('2) modulo 2^+\ (14.8) 

1=0 
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with 

a^,p+l{z) := a^^z) + 2\p+^{z), z = 0, 1, . . .,4^+1 - 1, (14.9) 

where the coefficients ai,i3{z) are supposed to provide a solution 

4Q+l_l 
1=0 

to (114. 4p modulo 2^^. This Ansatz, substituted in (114. 4p . produces a congruence of the 
form 

4"+l-l 

ns^iz)) +2^ Y. (^^./3+i(^) + ^'%p+ii^) + (^ + m+i,p+i{z))mz) + 1 = 



i=0 



modulo 2^+\ (14.10) 



where T{Sp{z)) represents terms that only depend on Si3{z). Inductively, we have al- 
ready computed Spi^z), and we know that T{Sp{z)) must be divisible by 2^. Comparison 
of powers of $2(2;) then yields a system of congruences that is equivalent to a system 
of the form 

&^,/3+l(^) + ^''&•,;3+l(^) + PoM2) = modulo 2, z = 0,l,...,4"+i-l, 

where Polj(2;), i = 0,1,..., 4"+^ — 1, are certain Laurent polynomials with integer 
coefficients. By Lemma [T2l these equations are solvable for any polynomials Polj(2;). 
Thus, we have proved that, for an arbitrary positive integer a, the modified algorithm 
that we have presented here will produce a solution S^a (z) to (114.41) modulo 2^° which 
is a polynomial in $2(^) of degree at most 4"+^ — 1 with coefficients that are Laurent 
polynomials in z. D 

Again, we have implemented the algorithm described in the above proof. For a = 1, 
that is, for the modulus 16, we obtain the following result. 

Theorem 37. Let ^2{z) = '^n>o^'^"^^^ ^^ before. Then, for the generating function 
Ss){5){z) = Yl'^=o^n+i{^{5))z"' for the subgroup numbers of S){5), we have 

7 



Sf,i5){z) = 8z^^ + Az^ + 82^ + 8z^ + 2z^ + 8z^ + Az + U+ — 






+ f 8^32/3 + 8^20/3 ^ 8U1/3 ^ 8^8/3 _^ 8^5/3 _^ ^ ^ ^ ^ ^ ^ _!_] $4(^) 



12 

y4/3 ' ^7/3 ' ^10j3 ' ^13/3 

(8 19 4 
8z' + 8z + - + - + -]^tiz) 
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+ (8^22/3 ^ g^l9/3 ^ g^l6/3 ^ g^lO/3 ^ ^2^7/3 + Sz^/^ 

^ _8_ ^ J2_\ s 

+^273 + ^573 + ^sTs + px^ ] 2(2;) 



^ ^ z4 ^ ^2 y -2V-; -r ^ ^4/3 ^ ^i3/3y -2 V-; ^ ^3/3 
+ ( 8^6 + 8;,^ + 8^ + 8 + - + 4 + 4 + 4^) <l>lHz) 



2; Z"^ Z'^ z^ ' 
+ (^ + ;|3+;^)^2^W+(;^ + ;i?7?)^2^W modulo 16. (14.11) 

We conjecture that Theorems [2T] and [5B1 extend to any Hecke group S){q), where q is 
a Fermat prime (note that PSL2{Z) = C2 * C3 = Sj{3)). 

Conjecture 38. For a positive integer h, let $h(-2) = X]n.>o '^^" /(2 -1), j^q^ a be a 
further positive integer, and let q = 2'^ + 1 be a Fermat prime. Then the generating 
function Ss^(^q-j{z) {see the first paragraph of Section^ for the definition)., reduced modulo 

2^" , can be expressed as a polynomial in $2/ (-2) of degree at most 2("+-^)2 — 1 with 
coefficients that are Laurent polynomials in ^^/(5-2) ^^g^ ^/^g integers. 

Note that Theorems 121] and [36] are the special cases corresponding to / = and / = 1, 
respectively. In particular, we conjecture that the obvious extension of the algorithm 
described in the proofs of the two theorems would be successful modulo any 2-power. 
In more detail, given a Fermat prime q, the first step consists in deriving a recurrence 
relation for the Hadamard product of the sequences (^C2(^))„>n ^^^ (^c'9(^)/^0n>o" 
By the procedure explained in the paragraph containing (lll.2p (where we now use (IlLip 
with F = io ((?)), this leads to a Riccati-type differential equation for the generating 
function X]n>o"^"+i(-^(^)) ^" ^°^ ^^^ subgroup numbers of S){q). The open questions 
are whether it will be possible to complete the base step, and whether it will always 
be possible to carry out the subsequent iterative steps in (the variation of) our method 
or, if necessary, its enhancement outlined in Appendix [D] Given the description of the 
parity pattern of the subgroup numbers proved in [27] Theorem B] , it is highly probable 
that the first question has a positive answer. What the answer to the second question 
is, remains entirely open. (The reader should recall that, in Section [121 we met a case 
where our method worked initially, but then stopped to work for modulus 16). 

Appendix A. Expansions of powers of the 2-power series $(z) 
Recall the notation 

U (.y\ - ST ^ai2"l+a22"2+-+ar2"r 

J-^ai,a2,...,ar\'^) — / _, ^ 

ni>n2>->nr>0 

from Section [3] In this appendix we list the expansions of $^(z) for i^ = 5, 6, 7, 8 in 
terms of the series Hai^a2,...,ar{z) where all a^'s are odd, obtained by using the algorithm 
described in Section [3] (see (13. ip and the proof of Lemma [9]). Namely, we have 
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$^(^) = l6H5{z) - 40iJ3,i,i(^) - 40iJi,3,i(^) - 40iJi,i,3(^) + 120iJi,i,i,i,i(^) 
- 80/73,1(2) - SOHi^siz) + 240/Ji,i,i,i(2) + {20z - 90)H3{z) 
- {60z - 270)ifi,i,i(z) - {120z - l90)Hi^i{z) 

+ {25z^ - 125^ + 75) Hi{z) + 50z^ - 75z, 

<^%z) = 96H5,iiz) + 96/7i,5(2) + 807^3,3(2) - 240/J3,i,i,i(z) - 240iJi,3,i,i(2) 

- 240i7i,i,3,i(2) - 240ifi,i,i,3(2) + 720/7i,i,i,i,i,i(^) + 2A0H,{z) 

- 6001/3,1,1(2) - 600iJi,3,i(2) - 600iJi,i,3(2) + 1800iJi,i,i,i,i(2) 

+ (120;^ - 840)i73,i(2) + (I2O2 - 840)/7i,3(2) - (360^ - 2520)/Ji,i,i,i(2) 
+ (300^ - 76A)H3{z) - (900^ - 2340)iJi,i,i(2) + (ISO^^ - 1200z + 1470)iJi,i(^) 

+ (375^2 - IO2O2 + 525)Hi{z) - 61z^ + 495^^ - 525z, 

$^(2) = -272H^{z) + Q72H5,i,i{z) + 672iJi,5,i(2) + 672iJi,i,5(2) 

+ 560/J3,3,i(2) + 5607/3,1,3(2) + 560iJi,3,3(2) + 2016i75,i(2) + 201QH,^5{z) + 16801/3,3(2) 

-1680iJ3,i,i,i,i(2)-1680iJi,3,i,i,i(2)-1680iJi,i,3,i,i(2)-1680/Ji,i,i,3,i(2)-1680i7i,i,i,i,3(2) 

+ 5040iJi,i,i,i,i,i,i(^) -5040/73,1,1,1(2) -5040/Ji,3,i,i(2)-5040iJi,i,3,i(2)-5040iJi,i,i,3(2) 

+ 15120iJi,i,i,i,i,i(2)-(3362-3360)iJ5(2) + (8402-8400)iJ3,i,i(2) + (8402-8400)iJi,3,i(2) 

+ (840z - 8400)i:fi,i,3(2) - (2520^ - 25200)/Ji,i,i,i,i(2) + (2520^ - 9408)1/3,1(2) 

+ (2520^ - 9408)/Ji,3(2) - (7560^ - 28560)/Ji,i,i,i(2) - (350^^ _ 40502 + 7434)iJ3(2) 

+ (1050^2 - 12180Z + 23310)i7i,i,i(2) + (3150^^ - 13020z + 13230)i7i,i(2) 

- (427z^ - 5040^2 + 9555z - 4347)i7i(2) - 1281/ + 5208^ - 4347^, 

^^{z) = -21761/7,1(2) - 2176ifi,7(2) - 17921/5,3(2) - 17921/3,5(2) 

+ 5376115,1,1,1 (2) + 537611i,5,i,i(2) + 537611i,i,5,i(2) + 537611i,i,i,5(2) + 4480113,3,1,1(2) 

+ 4480113,1,3,1(2) + 4480113,1,1,3(2) + 448011i,3,3,i(2) + 448011i,3,i,3(2) + 4480111,1,3,3(2) 

- 13440113,1,1,1,1,1(2) - 13440111,3,1,1,1,1(2) - 13440111,1,3,1,1,1(2) - 13440111,1,1,3,1,1(2) 

- 13440111,1,1,1,3,1(2) - 13440111,1,1,1,1,3(2) +4032011i,i,i,i,i,i,i,i(2) - 7616117(2) 

+ 188161l5,i,i(2) + 1881611i,5,i(2) + 1881611i,i,5(2) 

+ 15680113,3,1(2) + 15680113,1,3(2) + 15680111,3,3(2) - 47040113,1,1,1,1(2) 

-4704011i,3,i,i,i(2) -4704011i,i,3,i,i(2) -4704011i,i,i,3,i(2) -47040111,1,1,1,3(2) 
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+ 141120/7i,i,i,i,i,i,i(2) - (2688^ - 36288)^^5,1 (z) - (26882 - 36288)Hi^5{z) 

- (2240^ - 30240)if3,3(z) + (6720^ - 90720)i73,i,i,i(z) + (6720^ - 90720)i7i,3,i,i(2) 
+ (67202-90720)i7i,i,3,i(2) + (6720z-90720)ifi,i,i,3(2)-(201602-272160)iyi,i,i,i,i,i(z) 

- (9408z - 4:726A)H5{z) + (23520z - 119504)i73,i,i(2) + (23520z - 119504)i7i,3,i(^) 

+ (23520^ - 119504)ifi,i,3(;z) - (70560z - 361200)/7i,i,i,i,i(2) 

- (2800^^ - 44240^ + 115304) Hs^i^z) - (2800^^ - 44240z + 115304) Hi^^^z) 

+ (8400^2 - 132720^ + 355320)iJi,i,i,i(z) - (9800^^ - 55832z + 80892)/73(z) 

+ (29400^=^ - 1688402 + 260820)ifi,i,i(2) - (3416^^ -55020;z2 + 154980^- 135982)/7i,i(2) 

-(11956/-684742^ + 1012062-41245)iJi(2) + 13852^-223582^ + 59961;z^-41245z. 

Appendix B. The coefficients in the differential equation (111.51) 

Here we provide explicit expressions for the coefficients in the Riccati-type differential 
equation (111.51) . when reduced modulo 16: 

Po{z) = 8z^' + 8^46 + I2z45 + 4^43 ^ ^2z« + Uz^o ^ 4^39 ^ ^^z's ^ 3^37 ^ 4^36 ^ 4^34 ^ 2^33 

+ ll2^^ + 62^" + Uz^^ + Uz^^ + 13z^^ + 6z^® + 9z^^ + \\z^^ + \z^^ + 7z^^ + 9z^^ + z^° 
+ 72^^ + 15zis + 142^^ + 122^^ + llz^^ + 2" + IO2" + 52^2 ^ 22" + 82^° + 92^ + 152^ 
+ 52^ + 132^ + 42* + 142^ + 2^ + II2 + 15, 
Pl(z) = 82*8 + 122*6 ^ 122*5 ^ 4^43 ^ 122*2 ^ 4^41 ^ 4^40 ^ 4^39 ^ 8^38 ^ 4^36 ^ ^3^34 ^ g^33 
+ 102^1 + z30 + 4^29 ^ 11^28 ^ 8^27 ^ ^3^25 ^ ^24 ^ g^23 ^ ^22 ^ ^21 ^ ^^20 ^ 14^19 



- 92*8 + 62*^ + 142*6 _^ 8^15 _^ 13^14 _^ g^l3 _^ 22^2 + 9^10 + 11^:9 + 62'^ + 826 + 6 



2 



52* + 32^+42+1, 



P2(z) = 122^1 + 22*9 ^ i22*« + 42*7 ^ 4^46 ^ 122** + 102*3 + 22*^ + 42*1 ^ 8^40 ^ 14^39 

+ 52^7 + 22'''6 + Iqz^'^ + 102^* + 82^3 + 102^2 + 52^* + Tz-"^" + 122^9 + 72^8 + 42^7 

+ 122^6 + 92^5 + 122^* + 22^3 + 142^2 + 102^1 + 122^" + 102*^ + 32**^ + 62*6 + 92!'^ 

+ 122** + 152*3 + 142^2 _^ 102*1 + 22^" + 122^ + 142^ + 132^ + 7z6 + 82^ + 2-3 + 22^, 
P3(z) = 6z52 ^ 8^50 ^ 12^48 ^ 8z*7 + 142*6 ^ 4^45 ^ 8z** + 6z*3 + 82*^ + 42** + 52*° + 623^ 

+ 122^*^ + 22^7 + 22^6 + 72^* + 62^3 + 72^1 + 142^0 + 82^9 + 52^^ + 14^27 + 8;j26 

+ 142^5 + 142^* + 22^2 + 8z2i + 122^" + 142^9 + 62*^ + 122^7 + 92^6 ^ 4^15 ^ 22*3 
+ 2z*2 + 152*° + 142^ 
V^{z) = 4z53 + 152^* + 4z*7 + 2z*6 + 2*5 + 92** + 152*^ + 42** + lOz*" + 4239 ^ g^38 ^ 4^37 

+ 22^6 + 6235 + 1023* + 11233 _^ 8232 ^ 8^31 + 10^29 + 12226 + 1222^ + 2223 + 4222, 
V^{z) = 13z5* + 152*» + 62*7 ^ 12^46 ^ 3^45 ^ 10^44 ^ 4^43 ^ 12^42 ^ iQ^41 ^ 2^39 + 2238 

+ 15^36 + I4z35, 

V^{z) = 122^1 + 22*9 ^ 12^48 ^ 4^47 ^ 4^46 ^ 12^44 ^ 10^43 ^ 22*2 + 42*^ + 82*" + 14239 
+ 5237 _,_ 2236 + 6235 + 1023* + 8233 + 10232 + 523* ^ ^^30 ^ 12^29 _^ 7^28 j^ 4^27 
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+ 12z26 + 9z25 ^ i2z24 + 2z23 + 14^22 + 10^21 + I2z2" + 10^1^ + 3z^^ + 6z^^ 
+ 9z^^ + I2z^^ + 15z^^ + Uz^^ + lOz" + 2z^° + Ylz^ + 14z^ + 13z^ 

+ 7z6 + 8z^ + z^ + 2z^ 

Vl{z) = 12z'" ^ 13^51 ^ ^2z47 + 6z*6 _,. 3^45 ^ ^^^44 ^ ^3^42 ^ ^22^1 + 14^40 ^ y^^i-A ^ 2^38 

+ 12z3^ + 62^*^ + 2z3''^ + \^Z^^ + Z^3 + 82^2 ^ g^31 ^ ^4^29 _^ 4^26 ^ 4^25 ^ g^23 ^ ^2^22^ 
P8(^) = 2z52 + Sz^O + 4z48 + 8z47 + 10^46 + 12^45 + 8^44 ^ 2^43 + 8^42 + 12^41 + 15^40 ^ 2^39 

+ 4z38 + 6z3^ + 6z36 + Sz^"' + 2z33 + 5^31 + lOz^^" + Sz^^ + iSz^* + lOz^^ + Sz^e 
+ lOz^s + 10z24 + 6z22 + Sz^i + 4z2" + lOz^^ + 2z^^ + 4zi^ + \\z^^ + 12zi'"^ 
+ 6z" + 6zi2 ^ 13zi" + lOz^ 
P9(z) = 8z53 _^ loz^i + 8z47 + I2z46 + 6z45 + Gz^^ + 10z42 + gz^i + 12z4" + gz^^ + 4z38 + Sz^^ 

+ 12z3^ + 4z^'^ + 12z^'' + 2z^^ + 12z^^ + Sz^® + Sz^'^ + Ylz^"^ + 8z2^ 
pio(z) = 2z54 + ez^s + 12z47 + 8z46 + I4z45 ^ 4^44 ^ 3^43 ^ 3^42 ^ 4^41 ^ 4^39 ^ 4^38 

+ 6z^^ + 12z^^ 
pil(z) = 3z54 ^ ^48 _^ 10^47 ^ 4^46 ^ ^3^45 ^ g^44 ^ ^2z43 + 4z42 + Gz^i + Mz^^ 

+ 14z^** + z^*^ + 2z^^ 
Pi2(^) - 6z52 ^ g^50 ^ i2z48 + 8z47 + I4z46 + 4z45 + gz^^ + 6z43 + 8z42 + 4z" + Sz^" + Gz^^ + 12z38 

+ 2z^^ + 2z^'^ + Tz^-* + %7?^ + Iz^^ + 14z^" + Sz^^ + hz^^ + 14z2^ + Sz^'^ + 14z2^ + 14z^'^ 

+ 2z22 ^ 8^21 _^ j^2z2° + Wz^s + 6zi* + 12zi^ + Oz^^ + ^z^^ + 2z" + 2zi2 ^ ;l5^io _^ ^4^9^ 
Pi3(2) - 12z'^i + Sz^s + 4z45 + 4z44 + 12z42 ^ 3^40 ^ 3^38 ^ 3^36 ^ 3^35 ^ 3^34 ^ ^2z33 

+ 8z29 + 8z23, 
Pi4(z) = 2z54 + 6z48 + 12z4^ + Sz^e _^ 14^45 ^ 4^44 ^ 3^43 ^ 3^42 ^ 4^41 ^ 4^39 

+ 4z38 + 6z36 + I2z35, 
Pi5(z) = 2z54 + 6z*« + 12z4^ + 8z46 + 14z45 + ^z^^ + gz^^ + 8z42 + ^z^^ + 4z39 

+ 4z38 + 6z36 + 12z35, 
Pi6(z) = 4z53 + 15z5i + 4z47 + 2z46 + z45 ^ 9^44 ^ ^5^42 ^ 4^41 ^ ^^^m ^ 4^39 _^ g^38 ^ 4^37 _^ 2z36 

+ Gz^^ + lOz^'^ + llz-''^ + 8z32 _^ g^3i ^ ^Q^29 _^ j^2z2'^ + 12z2-^ + 2z^^ + 4z2^ 
pi7(z) = z54 ^ ^^^48 ^ 14^47 ^ ^2z46 + ISz^^ ^ 2z44 + 4z43 + 12z42 + 2z4i + lOz^*^ + 10z38 

+ llz3'^ + 6z3^ 
Pi8(z) = 13z54 + 15z4« + 6z47 + 12z46 + 3z45 + lOz^^ + 4z43 + I2z42 + lOz^i + 2z39 + 2z38 

+ ISz^^ + 14z3^ 

Appendix C. The coefficients in the differential equation (112. 3p 

Here we provide explicit expressions for the coefficients in the Riccati-type differential 
equation f ll2.3p . when reduced modulo 16: 

q^{z) = 8z50 + 8z49 + 4z47 + I2z46 + gz^^ + I0z43 + 2z42 + 2z« + 12z40 + 12z38 + Tz^^ + llz^e 
+ lOz^^ + llz^-* + 13z33 + ISz^^ + 14z^i + ISz^o + lOz^^ + 14z2** + %z^'^ + lOz^*^ + Oz^^ 
+ z23 + 4z22 + 8z2" + Mz^^ + 12zi** + Sz^^ + Oz^*^ + ^z^^ + ^z^^ + Sz^^ + z^^ + 14^11 
+ 9zi° + 2z^ + 14z* + lOz^ + 4z^ + llz^ + 4z'' + 4z3 + z^ + 8z + 15, 

qi(z) = 8z5i ^ 8^50 ^ 8^49 ^ i2z48 + 4z47 + I2z46 + 4z45 + Mz^s + I0z42 + Sz^o + 4z39 + 4z38 
+ 9z3^ + 4z36 + 14z^^ + hz^^ + 14z^^ + 13z^^ + 5z^° + ^z^"^ + 14z28 + 6z^^ + Sz^^ 
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+ lOz^s + 12z23 + 9z22 + 12z2i + 13z2" + 6z^^ + z^^ + Tz^^ + Ibz'^^ + llz^"^ + 7z^^ 
+ 5zi2 + ;j" + 6;j^° + lOz^ + 2z^ + Az'^ + Sz^ + 2z'^ + llz'^ + Az^ + Sz^ + 7z + 1, 
<Z2(z) = 8z'" + 82^6 + 82^5 + 8;^^^ + 8;2^3 + 12z''^i + 10z^° + 122^9 + 82^*^ + lOz^^ + ISz^*^ + 42^^ 

+ 10z43 + z42 + 11^41 + ^40 ^ 2z39 + z3« + Sz^^ + Bz^S + 7z35 + 8z34 ^ ^33 ^ g^31 ^ 9^29 

+ 14^28 + Uz^^ + 5z26 + 2z24 + 6z23 + 4z22 + 6z2i + 13z20 + ^i^ + lOz^s + z" + lOz^^ 
+ 12zi4 + z^-'' + 9zi2 + 9z" + 7z^" + 4z^ + 12z® + 3z^ + 2z^ + 9z^ + Sz'^ + 9z^, 

q^{z) = 14z56 + 2z55 + 14z54 + llz'^^ ^ g^51 ^ ^q^50 ^ ^^^49 ^ 2^48 + 15^47 ^ g^46 ^ ^Q^45 ^ g^44 
+ 4z43 + 5z42 + 11^41 + 5^40 ^ ^0^39 ^ 4^38 ^ ^2z37 + 2z36 + hz^^ + z^^ + 6z33 + 5z32 

+ ISz^i + 12z3" + Sz^^ + lOz^s + ISz^^ + 14z26 + Az"^^ + 12z24 ^ 3^23 ^ 9^22 _^ g^2i 
+ Wz^o + 7zi^ + Sz^* + 4zi^ + 4zi6 + 8zi4 _^ 2^5^13 ^ ^^3^11 ^ 9^9 ^ y^s + ^^ + 4^5^ 
(74(2) = Sz^*^ + ISz'"^* + 14z" + 12z''^*^ + 9z^4 ^ 7^53 _^ 9^52 _^ j^q^si ^ ;^2z^° + z49 + 10z4S + lAz'^^ 

+ 6z45 ^ g^44 ^ g^43 ^ ^4^42 ^ ^^^41 ^ ^3^40 ^ ^o^39 ^ ^38 ^ ^3^37 ^ 7^36 ^ 2^34 

+ 12z^^ + 4z^2 + 13z^^ + 13z-''° + 4z2^ + llz^** + Tz^^ + ^z"^^ + 2z'^^ + Bz^'^ + ISz^^ 

+ llz^i + 2z2" + 9zi9 + 6zi» + 2zi^ + 15zi6 + Sz^^ + 3zi4 + Mz^^ + 2zi2 + lOz" 

+ 9zl° + 6z^ + 12z^ 
(Z5(z) = 5z59 + 13z5*^ + 12z56 ^ ^55 ^ 8^54 ^ 5^53 ^ ^2z52 ^ ^5^51 ^ g^so ^ ^4^49 ^ 2^48 + 4z47 

+ 4z46 + 10^45 ^ ^2z44 + 8z43 + 4^41 ^ 2^4" + 8z38 + 12z37 + 2z36 + 12z35 + 8z34, 
q^z) = 12z49 + 4z48 + 12z46 + 12z45 + 12z44 + 4z4--' + 12z42 + I2z4i + 12z39 + 4z38 + 4z37 

+ 4z-''^ + Az^"^ + 4z-''^ + 12z2^ + 12z2'' + Uz"^^ + 12z2^ + 12z2" + 12z^'^ + 12z^^ + 4z^^ 

+ Az^^ + 12zi2 + 12z" + 12zi° + 4z^ + 4z*^ + 4z'^ + 12z^ + 4z^ 

qji^z) = lOz'^*^ + ez'^'^ + lOz^^ _^ g^52 ^ g^49 ^ 2^41 + 14;,40 _^ g^38 _^ ^4^37 ^ -^q^36 _^ ^4^31 

+ Wz^o + 2z2« + lOz^^ + 14z23 + 2z2i + Bz^^ + lOz^^ + lAz^^ + 2z^4 + g^io^ 
<Z8(z) = lOz^e + 6z55 + lOz^^ + 9z52 ^ 2^51 + Mz^o + 9z49 + Bz^^ + Sz^^ + 2z46 + ^4^45 ^ 3^44 
+ 12z43 + 7z42 + 9z4i + 7z4" + 14z39 + 12z3S + 4z-''^ + Gz^'^ + 7z'^ + Uz^^ + 2z'^ 
+ 7z^^ + ISz^i + Az^" + z^^ + 14z2* + 15z2^ + lOz^^ + \2z^'^ + 4z24 _^ 3^22 ^ c^^2\ 
+ lOz^o + 13zi^ + 7zi* + 12zi^ + 12z^*^ + Sz^^ ^ ^b^w ^ 3^9 ^ ^^3^8 _^ j^j^^e ^ J^2z^ 
(79(z) ^ 7z''^^ + Voz^^ + llz^'"^ + 7z''^-'' + Sz'"^! + lOz'^^ + 6z4S + 14z45 + ^z^^ + 6z■''^ 
gio(z) = z5« + 7z54 ^ 9^53 ^ 7^52 ^ ^5^49 ^ 5^41 ^ 3^40 ^ ^5^38 ^ 3^37 ^ 9^36 ^ 3^31 ^ 3^30 

+ 5z2« + 9z27 + 3z23 + 5z2i + 7zi9 + z^^ + H^is + ^3^14 ^ 7^10^ 
gii(z) = 12z55 + 4z5i + 12z5" + 4z47 + 12z46 + 4z45 + 12z43 + 4z42 + Az'^^ + Az^^ + 12z34 

+ Az^^ + 4z-''2 _^ 22z29 + 4z26 + 12z25 + \2z^'^ + 12z^® + 12z^^ + 4z^'^ + 12z^3 

+ 12z^ + 4z*+4z6, 
gi2(z) = 15z'^*^ + 14z'" + 12z^6 + 9z''^'^ + 7z^3 ^ 9^52 _^ j^q^si ^ j^2z'™ + z40 + 10z4S + 14z46 + %z^^ 

+ 6z44 + 6z43 + I4z42 + 11^41 + 13z40 + 10z39 + z3« + 13z-''^ + 7z36 + 2z34 + 12z33 + 4z32 

+ 13z3i + 13z3" + 4z29 + llz^s + 7z2^ + 2z25 + 6z24 + 13z23 + Uz^i + 2z20 + 9zi^ + Bz^* 

+ 2zi^ + ISz^s + Sz^^ + Zz^^ + 14zi3 ^ 2zi2 ^ j^Qz" + 9zi° + 6z^ + 12z*, 
gi3(z) = lOz'^^ + lOz'"^* + 2z^''^ + Wiz^"^ + Wz'^^ + 12z49 + 4z4S + 4z45 + 4z40 + Az^^ , 
quiz) = 8z'" + 8z56 ^ 8^55 ^ ^2z53 ^ 8^51 ^ 82^" + 8z49 + 8z48 + 10z47 + 14z46 + lOz^s + 8z44 

+ 9z43 + 12z42 + 2z4i + 6z4" + llz39 + 3z38 + 12z37 + 8z-''6 + 7^35 ^ 2^34 + ^33 + z32 
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+ Tz^l + Vlz^^ + 12z29 + z28 + ;j26 ^ ;^5^25 _^ 9^24 _^ j^4^23 ^ ^22 ^ 3^21 ^ 3^20 _^ 4^19 

+ 7z^^ + 7z^^ + 6z^^ + Sz^'^ + z" + 3z^2 + 2z" + 14z^ + 2z* + Gz'^ + 2z'^ + Mz"* 

+ 6z^ + 14z^ 
gi5(z) - 12z55 + 12z54 ^ g^50 ^ 2z49 + Mz^* + gz^" + lOz^e + lOz^^ + lOz^^ + Mz^^ + lOz^^ + lOz^^ 

+ 2z-''^ + 14z3^ + Gz^^ + 14z3'^ + 12z3'* + \^z^'^ + \^z^^ + 4z3" + Sz^*^ + lOz^^ + 12z2^ 

+ 10z24 + 12z23 + 2z22 + 2z2i + 2z20 + ^z^"^ + lOz^*^ + Sz^^ + 4zi6 + 6z^^ + Sz^'^ + Wz^^ 

+ 2zi2 + 2z" + 2zi° + 6z^ + 14z*^ + 14z'^ + 2z^ + 14z^ + 12z'^, 
gi6(z) = 4z6^ + lOz" + 6z^6 + Qz^'"^ + 8z''^3 ^ q^52 ^ -^-^^51 ^ ^^3^50 _^ 2z'*9 + 12z''^ + Sz"*^ + ISz"*^ 

+ 3z45 + 10z44 + 13z43 + 3z42 + 3z4i + 2z4" + 4z^^ + 7z'^ + Az^' + Gz^^ + Gz^^ + 5z34 

+ 15z33 + 3z32 _^ j^Q^Sl ^ ^4^30 _^ ^29 _^ 2^28 _^ 2^4^27 ^ ^^^26 _^ g^25 _^ 22^^ + 9^22 ^ g^21 

+ 2z20 + 8zi« + 10zi« + 12z" + 9zi6 + Sz^^ + Uz^^ + Sz^^ + gz^^ + lOz^" + z^ + 7z« 
+ Qz'^ + 7z^ 

gi7(z) = 8z'^« + 12z'^5 + 8z54 ^ g^53 ^ g^52 ^ 4^51 ^ ^2z5° + Sz^^ + g^48 ^ 4^47 ^ ^2z45 ^ g^43 

+ 8z42 + 12z4i + 8z4" + 12z39 + 8z38 + 4z36 + Az^^ + 8z-''4 + 4^33 ^ i2z32 + 8z3i ^ 4^29 

+ 8z2^ + 8z2® + 8z22 + 12z2^ + 8z^^ + 12z^^ + 8z^^ 
gi8(z) = 6z6" + Uz'"^^ + 4z" + lOz^^ + llz"^^ + 2z^'' + Az"^^ , 
gi9(z) = 4z'^6 + 12z55 + 4z'^4 + 2z'^2 + 4^51 ^ ^2z5" + 2z^'> + I2z^'' + lOz^^ + Az^"" + I2z^'> + Gz^^ 

+ lAz'^'^ + 2z'*i + 14z4" + 12z39 + 12z36 + Mz^^ + Gz^^ + Az^^ + 14z32 + lAz^^ + 2z29 

+ 12z2* + 14z2^ + 4z^*^ + 6z^^ + Az^^ + 4z^" + lOz^^ + 14z^^ + lOz^'"' + 14z" + 6z^ 

+ 10z* + 6z^ 

q2o{z) = 12z65 + 14z'" + 2z56 + 3^55 ^ ^^52 ^ 9^51 ^ ^5^50 ^ g^49 ^ 4^48 ^ ^47 ^ ^^^46 ^ ^45 

+ 14z44 + 15z43 + z42 + z^i + 6z4" + 12z-''s + ISz^* + I2z^' + 2z^'' + 2z^'> + Tz^'' + 5z'^ 
+ z^^ + Wz^i + lOz-'^o + Uz^^ + 6z28 + lOz^^ + 9z26 + Sz^^ + Gz^^ + Sz^^ ^ gz^" 
+ Wz^^ + 4z" + Sz^s + 7zi^ + 9z^^ + 7zi3 _^ 14^10 ^ ;l;l^9 ^ ;^3^8 ^ 2z'^ + 13z^ 
921 (z) = 2z6" + Oz'"^^ + lAz^^ + 9z^'' + 6z^^, 

q22{z) = 3z5« + 5z54 ^ ^^^53 ^ 5^52 ^ ^3^49 ^ 15^41 ^ g^40 ^ ^3^38 ^ 9^37 ^ ^^^36 ^ g^31 

+ 9z3° + 15z28 + Uz^^ + 9z23 + 15z2i + Sz^^ + Sz^^ + ^is ^ 7^14 ^ 5^10^ 
923 (z) = Wz'"^*^ + 12z'" + 2z54 + 14^53 ^ 2z52 + 4^51 ^ 2z49 + 4z48 + 12z46 + I2z45 ^ ^2z^4 ^ ^2^43 

+ 12z42 + 6z4i + 10z4" + 4z-''9 + 2z-''« + lOz^^ + Uz'^ + Az^^ + lOz^i + lOz^o + 6z^' 

+ 14z27 + 4z25 + 12z24 + 10z23 + 6z2i + Az^" + 2zJ^ + I2z^^ + 4z" + lAz^^ + lOz^^ 

+ 6zi4 + 12zi3 + 4zi2 + 4z" + 2zi° + 12z^ 
924(2) = 4z''^^ + 4z''^*^ + 4z''^^ + 4z^3 ^ ;l2z5\ 

925 (z) = 12z'^^ + 4z'^i + 12z^" + 4z''^ + 12z''-^ + 12z'^^ + 12z3'^ + 4z-''^ + 4z-'''^ + Az^^ + I2z^^ + Az^'^ 

+ 12z2i + 12z", 

926 (z) = 12z6" + ez'"^^ + 4z^6 + Gz^'^ + 4z5'', 

8z'" + 8z56 + g^55 ^ g^54 ^ g^53 ^ 4^51 ^ ^4^50 ^ 4^49 ^ g^48 ^ ^4^47 ^ 5^46 ^ i2z45 
+ 14z43 ^ ^^^42 ^ g^41 ^ ^^^40 ^ g^39 ^ ^^^38 ^ ^37 ^ 2^36 + 13^35 ^ g^34 ^ ^^^33 ^ 2, 

+ 3z29 + lOz^s + 9^27 + 7z2^ + 6z24 + 2z23 + Uz"^"^ + 2z2i + ISz^o + llz^^ + 14zi^ 
+ llz^^ + lAz^^ + Az^'^ + llz^^ + 3z^2 _^ 3^11 ^ ;l3^io _^ j^2z^ + 4z® + z'^ + Bz'^ 



(727(2 



;z3i 



66 M. KAUERS, C. KRATTENTHALER, AND T. W. MULLER 

+ 3z5 + z^ + 3z3, 

q2s{z) = 42^6 + i2z55 ^ 4^54 ^ 2z52 ^ 4^51 ^ ^2z''° + 22*9 + 12^48 + lOz^^ + 42^*^ + 12^45 + Qz*" 
+ 8z*3 + 14z42 + 2z4i + Wz^o + I2z39 + gz^^ + 8z'' + 12^36 + u^^^ ^ g^34 ^ 4^33 
+ Uz^^ + 142^1 + 8z^" + 2z29 + 12z2« + 14z2^ + 42^6 + 8z^^ + Sz^"^ + 6z^^ + iz^^ 
+ 4:z^° + lOz^^ + 142^^ + Sz^'' + 8z^^ + 10z^3 + 142^^ + 6z^ + lOz* + 6z^ + 8z^ 

529(2) = Sz^^ + 15z^^ + Uz^'' + 122'"^*^ + 'dz^'^ + 7z^3 + 9z''^2 _^ j^Qz'"^! + 12z-^° + z^^ + lOz"** + Wz'^^ 

+ 6z45 + 6z44 + 6z43 ^ ^4^42 ^ ^^^41 ^ ^3^40 ^ ^Q^39 ^ ^38 ^ ^3^37 ^ ^^36 ^ 2z34 
+ 12z33 + 4z32 _^ 2^3^31 ^ ;^3^30 _^ 4^29 _^ j^j^^28 ^ y^27 _^ g^26 _^ 2z25 + Qz'^'^ + I3z23 

+ llz^i + 2z2" + 9zi9 + 6zi« + 2zi^ + 15zi6 + Sz^^ + Sz^^ + Mz^^ + 2zi2 + lOz" 
+ 9zl° + 6z^ + 12z^ 
930(2) = 12z59 + 12z5« + 12z-^5 ^ ;^2z53 ^ 4^51 ^ 3^49 ^ 3^48 ^ 3^45 ^ g^40 ^ g^36^ 

931(2) = 4z55 + 12z5i + 4z5° + I2z^' + 4z46 + 12z45 + ^z^^ + 12z42 + I2z« + 12z38 + Az^^ + 12z33 
+ 12z^2 + Az'^'^ + 12z2*^ + Az^"^ + 4z22 + 4z^*^ + 4z^^ + 12z^'' + 4z^^ + 4z^ + 12z*^ + 12z^ 

'732(2) = Sz'"^* + lOz'" + 4z^6 + Sz^"* + 13z'^3 ^ 3^52 ^ -^4^51 _^ 4^50 _^ j^j^^49 ^ ^^4^48 _^ j^q^46 ^ 2z'^5 
+ 2z44 + 2z« + 10z42 + 9z4i + iSz^o + Mz^^ + llz^^ + ISz^^ + ISz^^ + Gz^^ + 4z33 
+ 12z32 + ISz^i + ISz^" + 12z29 + Oz^*^ + ISz^^ + Bz^^ + 2z2'^ + ISz^'"' + Oz^^ + Gz^" + Sz^^ 
+ 2zi« + 6zi^ + 5zi6 + 7zi'^ + zi4 + lOz^-^ + Gz^^ + 14z" + Sz^" + 2z'' + 4z«, 

g33(z) = 12z'^^ + 12z''^* + 12z-^'^ + 12z'^3 ^ 4^51^ 

934(2) = 2z'^5 + lOz^^ + 4z5° + 3z4" + 13z4« + 4z4^ + Tz^^ + 7^45 ^ ^5^44 ^ 5^43 ^ ^5^42 ^ ^5^41 

+ Sz^O + llz39 + 13z3S + z37 + 5z36 + 8^35 ^ 2z34 + 5z32 ^ 5^31 ^ 14^30 ^ 8^29 ^ 4^28 

+ ISz^^ + lOz^^ + ISz^'^ + lOz^^ + llz^^ + llz^^ + llz^" + 3z^^ + 15z^® + 4z" + 14z^^ 

935(2) = 12z65 + 14z'" + 2z56 + 3^55 ^ 8^53 ^ 2z52 + Qz^i + 15z5° + Gz^^ + 4z4« + z^' + ISz^^ 
+ z^s + I4z44 + 15z43 + z42 + z^i + 6z4" + 12z39 + 13z38 + Uz^'^ + 2z36 + 2z^'> + 7z^^ 
+ 5z^^ + z^^ + lAz^^ + lOz^" + llz^s + 6z2® + lOz^^ + 9z^^ + Sz^'^ + 6z^^ + 3z^^ + 8z^^ 
+ 6z20 + 8zi9 + 14zi8 + 4z" + 3zi6 + 7z^^ + 9z^^ + 7z^' + 8z'^ + lAz^° + Hz" + 13z8 
+ 2z■^ + 13z^ 

536 (z) = 8z67 + 12z5*^ + Wz^s ^ 4^54 ^ 4^53 ^ 4^52 ^ 2z5i + Gz^o + 12z49 + 12z48 + lOz^^ + 6z^^ 
+ 8z44 + i2z43 + 12z42 + Gz^i + 12z40 + Bz^^ + 4z38 + 2z36 + lOz^^ + 12z34 + 2z33 
+ 6z32 ^ -^2z3i + lOz^^ + 12z2s + 12z2*^ + 8z25 + Uz^^ + lAz^^ + 8z^^ + Az?^ + 14z" 
+ 4zl6 + 8zl^ 

937(z) == 4z6" + 2z''^^ + 8z'" + 12z^6 + 2z''^^ + 12z''^'^ + 8z^^ , 

q^siz) = 5z5*^ + lOz'" + 4z56 ^ 3^54 ^ ^3^53 ^ 3^52 ^ ^4^51 ^ 4^50 ^ ^^^49 ^ ^4^48 ^ ^oz^e + 2z^'' 

+ 2z44 + 2z« + 10z42 + 9z4i + I5z40 + Mz^^ + Hz"''® + ISz^^ + 13z36 + 6z'^ + Az^' 
+ 12z^^ + ISz^i + ISz^" + 12z29 + Oz^*^ + 13z2^ + 6z^^ + 2z^^ + 15z^^ + 9z^^ + 6z^" 
+ 3zi9 + 2zi« + 6zi7 + 5zi6 + 7zi5 + z^^ + lOz^^ + 6z^^ + 14z" + 3zi° + 2z'> + Az\ 
qsgiz) = 4z58 + lOz^s + 12z54 ^ ^2z53 ^ ^2z-^2 ^ g^5i ^ 2z5" + 4z49 + 4z48 + Mz^^ + 2z''^ + Az^^ 
+ Az^^ + 2z4i + 4z« + 2z39 + 12z38 + 6z'^ + lAz^" + Az'^ + 6z'^ + 2z'^ + Az^^ + lAz^^ 
+ Az^^ + Az^^ + Az^^ + lOz^i + 12zi® + lOz^^ + 12zi^, 
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q4o{z) = 12z^^ + 12z^^ + 12z^^ + 12z^^ + 4z^\ 

quiz) = 4z^" + 2z''^ + 12z''^ + 2z^'' + Uz''"^, 

qi2{z) = 22''^*^ + 14z55 + 2z54 ^ 5^52 ^ ^Q^51 ^ g^50 ^ 5^49 ^ ^4^48 ^ ^47 ^ ^Qz^^ + Qz^^ + 72^4 

+ 12z43 + llz42 + 5z4i + llz« + 6z39 + 12z3« + Az^' + Uz^^ + Uz^^ + l^z^^ + lOz^a 
+ llz^^ + Sz^i + 4z3" + ISz^^ + 6z28 + Sz^^ + 22^^ + 122^5 + Az"^^ + Sz^'"' + 72^^ + 10;?^^ 
+ 2z20 + 9^19 + 11^18 + I2z" + 12zi6 + Sz^'^ + z^^ + a^^n + y^^a + g^s + 15^6 + 12z^ 
943(2) = 8z66 + 3z5« + 6z'" + 12z56 ^ 5^54 ^ ^^^53 ^ 5^52 ^ 3^51 + 12^50 ^ ^3^49 ^ ^z^b ^ g^46 

+ 14^45 + 14z44 + 14^43 ^ g^42 ^ ^5^41 ^ 9^40 ^ 2^39 + 13^38 ^ g^37 ^ ^^^36 ^ ^Qz^^ 
+ I2z^^ + 4z32 _^ 9^31 ^ 9^30 ^ 4^29 ^ -^5^28 _^ j^j^^27 ^ g^26 ^ ^q^25 _^ j^4^24 ^ 9^23 

+ 15z2i + lOz^" + Sz^^ + Wz^^ + lOz" + Sz^s + z^'' + 7zi4 + Gz^^ _^ j^Qz^^ ^ 2z" 
+ 5z^° + 14z^ + 12z^ 
944(2) = 2z59 + 2z5« + 8z56 + loz^s ^ 2z53 ^ 3^52 ^ g^5i ^ i2z49 + 4z48 + gz^^ + gz^e + 4z45 
+ gz44 + 8^41 ^ 4^40 ^ 8^37 ^ 4^36 ^ 8^35^ 

(745(2) = 8z^\ 

946(2) = 6z59 + 6z58 + 14z55 + 6z53 + 2z5i + 4z49 + 12z48 + I2z^^ + 12z4" + 12z36, 

947(2) = 4z6^ + 2z'" + lAz"^^ + ISz^'"^ + gz'^-'^ + 14z^2 _^ 7^51 _^ ^50 _^ j^q^49 ^ ;^2z'^^ + ISz^^ + z^^^ 

+ 15z45 ^ 2z44 + ^43 ^ ^5^42 ^ ^5^41 ^ ^^^40 ^ 4^39 ^ 3^38 ^ 4^37 ^ ^4^36 ^ ^4^35 ^ 9^34 

+ llz^^ + 15z32 _^ 2z3i + 6z-''° + 5z29 + lOz^* + Gz^^ + Iz"^^ + ISz^^ + lOz^'^ + ISz^^ + gz^^ 
+ lOz^o + gzi9 + 2z^^ + 12z" + ISz^s + 9zi5 + Tz^^ + 92^3 + gz^^ + 2zi" + Sz^ + 3z« 
+ 14z■^ + 3z^ 

948(2) = 8z5« + 4z'^5 + gz'^4 + gz53 + 8z52 ^ ^2z'^i + 4z'^° + gz^^ + 8z48 + 12z4^ + Az^"" + gz^^ + ^^^2 
+ 4z'^i + gz'^o + 4z-''^ + gz^s + 12z36 + I2z35 + gz34 + Uz^^ + Az^"^ + gz^i + 12z29 + gz^^ 
+ gz^e + gz22 + 4z2i + gz^® + 4z" + gz^^ 

949(z) = 4z6" + lOz'"^^ + gz" + 12z^6 + lOz'"^-^ + 12z-^'* + gz'^^^ 

95o(z) = 4z''^^ + 4z''^* + 4z''^-^ + 4z^3 ^ 22z^\ 

951 (z) = 12z6" + lAz"^^ + 4z^6 + lAz"^^ + 4z^^ 

952(2) = gz^e + 7z5« + 14z'" + 12z56 ^ ^54 ^ ^5^53 ^ ^52 ^ ^g^si ^ ^2z5" + 9z49 + lOz^s + lAz^^ 
+ 6z45 + 6z44 + g^43 ^ ^4^42 ^ 3^41 ^ 5^40 ^ ^Oz39 + 9z-''« + bz^' + 15z36 + 2z34 + 12z33 
+ 4z32 + bz^^ + 5z-''° + 4z29 + 3z2« + 15z2^ + gz^^ + 2z'^^ + Bz^^ + Sz^^ + Sz^^ + 2z2" 
+ zi« + 6zi8 + 2zi7 + 7zi6 + 13zi-^ + llzi4 + 14zi-'' + 2z^^ + lOz" + z^" + 6z» + 12z«, 

953(2) = 4z59 + 4z5« + 4z55 + 4z53 + i2z5i + gz^^ + gz^s + gz^^ + gz^o + gz^e, 

954(2) = 4z6^ + lAz"^^ + gz" + ISz^'"^ + lOz''^'^ + 10z^3 _^ 2z''^2 _^ ^5i _^ ^^^w ^ ^^4^49 _^ g^48 ^ 5^47 

+ gz^e + 3^45 ^ i2z44 ^ g^43 ^ g^42 ^ 3^41 ^ g^40 ^ 3^39 ^ 2z3S + Qz^^ + Sz^^ + Gz^^ 

+ 9z33 + 3z-''2 _^ g^3i _^ 8z^''" + 13z29 + Gz^s + Bz^*^ + Az"^^ + gz^'^ + Wz^^ + ISz^^ + 4zi^ 

+ 2zi^ + ISz^^ + lOz^*^ + 4zl^ 
955 (z) = 6z6" + 3z''^^ + 4z'" + lOz^s + 3z''^^ + 2z''^'' + 4z^^ 
956(2) = 5z59 + 13z58 + 12z56 + z^s + gz54 ^ 5^53 ^ i2z52 ^ ^5^51 _^ 8z5" + 14z49 + 2z^^ + Az^^ 

+ Az^'' + lOz^s + 12z44 + gz43 + 4z4i + 2z4" + gz^* + 12z37 + 2z-''6 + ^^z""'' + gz^^, 
957(z) = 14z6" + ISz'"^^ + 4z" + 2z''^'^ + ISz'"^^ + lOz^^ + Az^"^ , 
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958 (z) = 8z^^ + 8z^s + 8z^'^ + 12z53 j^ 8^51 _^ g^so _^ §^49 ^ g^48 ^ ;^0z''^ + 14z'*s + lOz"^^ + Sz^"^ 

+ 9z43 + I2z^2 ^ 2^41 _^ g^40 ^ ^^^39 ^ 3^38 ^ ^22^7 + 8z36 + 7z35 + 22^4 + ^33 ^ ^32 
+ TZ^I + 12^30 ^ ^2z29 + ^28 ^ ^26 ^ ^5^25 _^ g^24 ^ ^^^23 ^ ^22 ^ 3^21 ^ 3^20 ^ 4^19 

+ 7z^^ + 7z^^ + Gz^'"^ + Sz^'^ + z" + 3z^2 + 2z^^ + Uz^ + 2z* + 6z^ + 2z^ + \Az^ 
+ 6z^ + 14z^ 
959 (^) = 4z56 ^ ^2z55 + 4z-^4 ^ 2z-^2 ^ 4^51 ^ ^22^" + 22^9 + 12^48 + 10^47 ^ 4^46 ^ ^22^5 + Qz^^ 
+ Wz'^^ ^ 2z'*i + 142^" + 12z39 + 122^6 + I4z35 + gz^^ + 42^3 + 142^2 + Mz^i + 22^9 
+ 12z28 + 142^^ + 42^*^ + 6^22 + 4z2i + 42^" + lOz^^ + 142^^ + lOz^''' + 14z" + 62^ 

960 (^) = 12z^' + 14z'" + 2z56 + 3z55 ^ 2z52 ^ 9^51 ^ ^5^50 ^ g^49 ^ 4^48 ^ ^47 ^ ^5^46 ^ ^45 

+ 14;244 ^ ^5^43 ^ ^42 ^ ^41 ^ g^40 ^ ^2^39 + 13^38 ^ ^2^37 + 2^36 + 2z^'> + 72^4 + 5^33 

+ z32 + Mz^i + lOz^o + 112^9 + 6z28 + lOz^^ + 9z26 + Sz^^ + Gz^^ + Sz^^ + 62^" + Mz^^ 
+ Az^'' + Sz^^ + Tz^'^ + 9z" + 7zi3 ^ ;l4^io _^ ^^^9 _^ ^3^8 ^ ^^7 ^ -^3^6^ 
961 (z) = 2z6" + Oz'"^^ + 14z'^*^ + 9z-^'^ + 6z''^'^, 

962(^) = 3z5« + 5z54 ^ ^^^53 ^ 5^52 ^ ^3^49 ^ ^5^41 ^ g^40 ^ ^3^38 ^ 9^37 ^ ^^^36 ^ g^31 ^ 9^30 

+ 15z28 + llz27 + 9z23 + 15z2i + 5zi9 + 3zi6 + ^i5 ^ 7^14 ^ 5^10^ 
963 (^) - 6z55 + 14z'^4 + 12z5° + 9z49 + Tz^* + 12z4^ + Sz^*^ + Sz^^ + ISz^^ + ISz^^ ^ ^3^42 
+ ISz^i + 8z4" + z39 + 7z38 + 3z37 + 15z36 + Sz^^ + 6z34 ^ ^5^32 ^ ^5^31 ^ ^Qz^o 
+ Sz^'^ + 12z2^ + ISz^^ + 14z2''^ + Vdz"^^ + 14z23 _^ ^22 _^ ^21 ^ ^20 ^ 9^19 _^ ^^3^18 
+ 12zi7 + lOzie + 3zi5 + 4zi4 + Tz^^ + ^12 + ^11 + ^10 ^ 3^9 ^ ^5^8 ^ 7^7 ^ ^6 
+ 15z^ + 6z'^, 
964(z) = 15z'^*^ + 14z'" + 12z-^6 + %z^'^ + 7z53 ^ 9^52 _^ j^q^si _^ j^2z'^" + z''^ + lOz^^ + 14z''^ + Qz^^ 
+ 6z44 + 6z« + 14z42 + Uz^i + 13z« + lOz^^ + z^® + ISz^^ + 7z^^ + 2z^^ + 12z33 + 4z32 
+ ISz^i + ISz^" + 4z29 + 11^28 + 7^27 _^ 2z25 + 6z24 + ISz^^ + Uz^i + 2z20 + 9zi^ + Bz^* 
+ 2zi^ + ISz^s + 5z^^ + iz^'^ + 14zi3 _^ 2zi2 _^ j^Qz" + 9zi° + 6z^ + 12z*, 
965 (z) = 12z'^* + lAz"^^ + 4z^4 + 4z'^3 _^ 4^52 ^ 2^51 + 6z^° + 12z'^^ + 12z4*^ + lOz'^^ + Qz'^^ + 12z'^3 

+ 12z42 + 6z41 + 12z4" + 6z39 + 4z38 + 2z36 + lOz^^ + 12z34 + 2z33 + 6z32 ^ ^^^Zl ^ ^q^29 

+ 12z28 + 12z2*^ + 12z22 ^ ;l4^21 _^ 4^18 ^ ^4^17 _^ 4^16^ 

966 (z) = 12z6^ + 14z'" + 2z^6 + 3z'^^ + Sz'^^ ^ 2z^2 ^ 9^51 ^ ^^5^50 _^ g^49 ^ 4^48 ^ ^47 ^ ^^5^46 _^ ^45 
+ 14z44 ^ ^5^43 ^ ^42 ^ ^41 ^ g^40 ^ ^2z39 + 13z3S + I2z^' + 2z36 + 2z35 + 7z34 ^ 5^33 
+ z32 + Mz^i + lOz^o + llz29 + 6z28 + 10z27 + 9z26 + Sz^^ + Gz^^ + 3z22 + gz^^ + Qz^° 
+ 8zi9 + 14zi8 + 4z" + 3zi6 + 7zi''^ + %z^^ + 7z^^ + 8z^^ + Uz^" + llz^ + ISz^ 
+ 2z■^ + 13z^ 

967(z) = 12z'^^ + 12z''^* + 12z-^'^ + 12z'^3 ^ 4^51^ 

968 (z) == 4z6" + 2z''^^ + 12z''^'^ + 2z^''^ + 12z5'', 

969 (z) = 4z^^ + Uz^^ + 8z" + Tz^'"^ + lOz^^ + lOz^^ _^ 2z'^2 _^ g^si _^ 3^50 ^ -^4^49 _^ g^48 _^ ^^3^47 
+ 8z46 + llz^s + 12z44 + 6z43 + 6z42 + Uz^ + 6z^° + llz^'^ + 2z^' + z^e + 13z35 + 6z'* 

+ z33 + llz32 _^ g^Sl _^ g^30 _^ 5^29 ^ g^28 ^ g^26 _^ 4^25 _^ g^24 _^ ^^4^22 ^ ^^21 _^ 4^19 

+ 2zl^ + 7zl^ + 10zl6 + 4zl^ 
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<?7o(z) = 142^" + I5z^^ + 4z5^ + 2z^^ + I5z^^ + lOz^^ + Az^^, 

qniz) = 4z51 + 14z50 + 4z49 + Mz^^ + 5^46 ^ ^2^45 ^ ^^4^43 ^ ^^^42 _^ 9^41 ^ -^^^40 ^ g^39 ^ ^^^38 

+ z^^ + 2z36 + Uz^^ + llz^^ + 2z^^ + Sz^^ + lOz^*^ + gz^'^ + Tz^^ + Gz^" + 22^^ + 122^2 
+ 2z2i + 15z20 + llzi9 + 14zi« + Hz" + 14zi5 + 42^4 + Uz^^ + 3^12 + 3^11 + 13^10 
+ 12z^ + 4z* + z'^ + Gz*^ + 3z-^ + z^ + 3z^ 
g72(z) = Sz^* + lOz" + 4z^6 + Sz^'' + ISz^'"' + Sz^^ + Wz^^ + llz''^ + Uz^^ + lOz^^ + 2z'^5 + 2z'^'' 

+ 2z43 + 10z42 _^ g^41 ^ ^5^40 ^ ^4^39 ^ ^^^38 ^ 15^37 ^ ^3^36 ^ g^34 ^ ^5^31 ^ ^5^30 

+ 9z28 + 13z27 + 6z25 + 2z24 + 15z23 + Qz^^ + Gz^" + Sz^^ + 2zi« + Gz^^ + Sz^^ + Tz^^ 

+ z^'' + lOz^^ + 6zi2 + 14z" + Sz^" + 2z^ + 4z^ 
573 (z) = 2z5^ + Gz^i + 10z'^° + Gz-*^ + lOz''^ + ez'*^ + lOz''^ + Gz^^ ^ q^h ^ 24z3« + lOz^'' + Uz^^ 

+ 6z^^ + 2z2^ + ez^^ + 2z^^ + 2z^2 + 2z^® + lOz^^ + Qz^^ + lOz" + 2z^ + 14z^ + 14z^ 
g74(^) = 2z55 + 10z54 + 4z50 + 3z49 + 13z48 + 4z47 + 7z^^ + Tz^s + ISz^^ + Sz^^ ^ ^5^42 ^ ^5^41 

+ llz^^ + 13z^* + z^^ + Sz^^ + 2z^'' + 5z^2 + 5z^^ + 14z^° + Az^^ + ISz^'^ + lOz^^ 

+ 15z24 + 10z23 + llz22 + llz^i + llz^o + iz^'^ + ISz^* + 4z" + Mz^s + z^^ + 12zi'^ 

+ 13z^^ + llz^2 + Hz" + llz^° + z^ + 5z^ + 13z'^ + llz*^ + 5z^ + 2z'', 
q7z{z) = lOz'^^ + 13z5^ + Wz^^ + 14z'" + 6z52 _^ 3^51 _^ ^50 _^ j^q^49 ^ 2z*** + ISz"*^ + Qz'^^ + 2z'*3 

+ 2z42 + 9z« + 2z« + 9z39 + 6z38 + Hz^e + iSz^^ + 2z^^ + llz^^ + 9z32 ^ 2z3i + Jz^"" 

+ 2z28 + 2z26 + 10z22 + 13z2i + Gz^^ + 13z" + lAz^^ , 

qreiz) - 5z52 ^ 5^49 ^ ^47 ^ ^^44 ^ ^^^42 ^ 5^41 ^ ^^^40 ^ ^^^35 ^ ^5^34 ^ ^^^32 ^ 3^31 ^ ^3^29 

+ 3z27 + 7z22 + 9zi9 + llzi« + zi3 + 3z" + Tz^^ + 9z8 + ISz^, 

g77(^) = 132^' + 7z51 ^ ^50 ^ ^5^47 ^ ^46 ^ ^5^45 ^ ^43 ^ ^5^42 ^ ^^^41 _^ 3^38 ^ 9^34 ^ ^^^33 

+ 15z^^ + Sz^^ + 7z2^ + 13z2^ + 13z22 + 13zi^ + 9z" + 7z^^ + 9z^^ + Sz^ + 3z^ + 3z^ 

Appendix D. The method "reloaded" 

As the reader will recall from Section H] (cf. Remark ITT]) . our method described there 
is based on the "hope" that, if a polynomial in $(2) is zero modulo a 2-power 2^ (as a 
formal Laurent series), then already all coefficients of powers of ^{z) in this polynomial 
vanish modulo 2^. (This is manifest in each comparison of coefficients of powers of 
$(2;) in Section m) In general, however, this implication does not hold (see Lemma 139) 
below for the case of modulus 2^ = 16). It may consequently happen that the method 
from Section H] fails to find a solution modulo 2^ to a given differential equation in the 
form of a polynomial in $(2) with coefficients that are Laurent polynomials in z over 
the integers, while such a solution may in fact exist. As it turns out, this situation 
occurs when treating the subgroup numbers of SL2{1') and of r3(3) modulo 16, see the 
paragraph above Theorem |2S] and Remark [32J (In the former case, there is indeed a 
solution, while in the latter there is not.) 

Our aim here is to explain how the method from Section H] can be enhanced so that 
one can decide whether or not such a solution modulo a given 2-power 2^ exists; and, if 
it exists, how to find it. In principle, it should be possible to describe such an improved 
method for an arbitrary 2-power 2^ . Since, in the present paper, we need it only for 
the modulus 16, and since we are not able to rigorously establish the validity of the 
enhancement we have in mind in general (it would depend on Conjecture HI which at 
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present we are not able to prove), we content ourselves with describing the enhanced 
method for the modulus 16. From this description, the reader should have no difficulty 
to "extrapolate" to arbitrary 2-powers, assuming the truth of Conjecture ID 

We begin by characterising when a polynomial in $(;z) with coefficients that are 
Laurent polynomials in z over the integers vanishes modulo 16 as a Laurent series in z. 

Lemma 39. As before, let ^{z) = Yln>o^'^" ■ Furthermore, let P{z,^{z)) be a poly- 
nomial in ^{z) with coefficients that are Laurent polynomials in z over the integers. 
Then, as a Laurent series in z, 

P(z,$(2))=0 modulo 16 

if, and only if, the coefficients of powers of ^ in P{z,^{z)) agree modulo 16 with the 
corresponding ones in 

c,{z)M,,{z,<l>{z))+2{c2{z)<l>{z)+cs{z))Ms{z,<^{z)) 

+ 8{c^{zMz) + c,{z))M2{z, $(^)). (D.l) 

Here, M2{z,t),Ms{z,t),MiQ{z,t) are the minimal polynomials for the moduli 2,8,16, 
respectively, given in Proposition^ and ci{z) , C2{z) , cs^z) , C4^{z) , c^i^z) are suitable Lau- 
rent polynomials in z over the integers. 

Proof. We assume that P{z, $(;z)) = modulo 16. 

Recall that, by definition, Miq{z,^{z)) is a monic polynomial in $(2;). We use this 
fact to perform division of P{z,^{z)) by Miq{z,^{z)) (as polynomials in ^{z)), thus 
obtaining 

P{z, $(^)) = c,{z)Mie{z) + Pi{z, $(^)), (D.2) 

where Pi{z, ^{z)) is a polynomial in ^{z) of degree at most 5, with coefficients that are 
Laurent polynomials in z over the integers, say 

P,{z, ^z)) = 4 W$' W + d,{z)<l>\z) + d,{z)<l>'\z) 

+ d2{z)<^\z) + di{z)<l>{z) + do{z). (D.3) 

As Laurent series in z, both P{z, $(-2)) and Miq{z, $(-2)) vanish modulo 16. Using this 
observation in (]D.2p . we see that Pi{z,^{z)) vanishes modulo 16 as well. Now recall 
from (12. 5p and the proof of Lemma [1] that 

Pi{z,<!>{z))=d,{z)5\E,{z)+Q,{z) 

(with a suitable series Qi{z)), where D{Qi{z),16;n) has strictly smaller asymptotic 
growth (in n) than D{E^{z), 16; n). Since, as we already observed, Pi{z, ^{z)) vanishes 
modulo 16, it follows that bld^lz) must vanish modulo 16, that is, there exists a Laurent 
polynomial C2{z) over the integers such that 0^5(2;) = 2c2{z). We use this observation in 
(1D.3P to see that 

P{z, <t>{z)) = Ci{z)M,,{z) + 2c2{z)^z)Ms{z, <^{z)) + P2{z, $(^)), (D.4) 

where P2{z, ^{z)) is a polynomial in ^{z) of degree at most 4, with coefficients that are 
Laurent polynomials in z over the integers, say 

P2{z, <^{z)) = e,{z)<!>\z) + es{z)<^^{z) + e2{z)<l>\z) + e,{z)^z) + eo{z). 
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Applying the same kind of argument again, we further deduce that 

P{z, <^{z)) = ci{z)M,,{z) + 2(C2(Z)$(2) + C3{z))Ms{z, $(^)) + Psiz, <^{z)), (D.5) 

where 03(2;) is a Laurent polynomial in z over the integers and P^iz, ^{z)) is a polynomial 
in ^{z) of degree at most 3, with coefficients that are Laurent polynomials in z over 
the integers, say 

P,iz, ^z)) = Mz)^%z) + Mz)<^'{z) + MzMz) + foiz). 

As Laurent series in z, all of P{z,^{z)), Mie(z, $(z)), and 2M8(z, $(2;)) vanish mod- 
ulo 16. Using this observation in fID.Sp . we see that Psi^z, $(2;)) vanishes modulo 16 as 
well. Equation (12.51) and the proof of Lemma [T] give that 

P3{z,<t>{z))=d,{z)3\E3{z)+Q,{z) 

(with a suitable series Q3{z)), where D (Qsi^z), 16; n) has strictly smaller asymptotic 
growth (in n) than D{Es{z), 16; n). Since, as we already observed, ^3(2;, ^{z)) vanishes 
modulo 16, it follows that Sld^lz) must vanish modulo 16, that is, there exists a Laurent 
polynomial 04(2) over the integers such that ^3(2;) = 804(2;). By another application of 
the same kind of argument, this leads to 

P{z, $(^)) = c,{z)M,,{z) + 2{c2{z)<^{z) + Cs{z))Ms{z, ^z)) 

+ 8{c^{zMz) + c,{z))M2{z, $(^)) + P^{z, ^z)), (D.6) 

where 05(2) is a Laurent polynomial in z over the integers and ^4(2;, $(2;)) is a polynomial 
in ^{z) of degree at most 1, with coefficients that are Laurent polynomials in z over 
the integers, say 

P,{z,^z))=g^{zMz)+go{z). 

Since P(z,^(z)), Mie(z,^{z)), 2M8{z,^{z))), 8M2{z,^{z)) all vanish modulo 16, also 
P4(z, $(z)) must have this property; but this means that gi{z) and go{z) both vanish 
modulo 16. 

If we combine flD.2p . flD.4p . fID.Sp . flD.6p . then we obtain our claim. D 

Now we put ourselves in the situation that we want to describe the coefficients of 
the formal power series F{z) modulo 16, where F{z) solves a Riccati-type differential 
equation of the form (14. ip . and that we try to solve this problem by expressing F{z) in 
the form 

5 

F{z) = ^ai{z)¥{z) modulo 16, 

i=0 

where the aj(2;)'s are Laurent polynomials in z over the integers to be determined. Let 

us assume that, while following the approach outlined in Section H] (with Miq{z, $(2;)) in 

place of the polynomial in (14.40 ). we have already reached the level of modulus 8, that 

is, that we have found Laurent polynomials Gq^^^z), 01^3(2;), 02,3(2;), 03,3(2;), 04,3(2;), 05,3(2;) 

such that 

5 

E«^'3(^)$^(^) 
i=0 
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solves the differential equation (14.1 p modulo 8. According to the Ansatz fl4.6l) - fl4.8p 
with ^ = 3, we now substitute 



^(a,,3(^) + 86,,4(^))$*(^) (D.7) 

1=0 

(where the 64,4(2;) 's are at this point undetermined Laurent polynomials in z) instead of 
F{z) in (14. ip . For the sake of better readability, in the sequel we write bo^z) for bQ^^^z), 
etc. After simplification of the left-hand side of (14. ip modulo 16 as described below 
(14. 8p . and after reduction of the resulting expression modulo Miq{z, $(-2)) (which is a 
polynomial in $(-2) of degree 6), we obtain a polynomial of the form 



t> 
8J2{pii^) + G^{zMz),^'{^)))^'{^), (D-8) 

i=0 

where the Pi{z)^s are certain Laurent polynomials in z over the integers, and the 
Gi(^z, h(z), h'{z)ys are certain linear forms in 

bo{z), bi{z), b2{z), bsiz), b4{z),b5{z) and bQ{z),b[{z),b2{z),b'r^{z),b'4{z),b'^{z), 

with coefficients that are (known) Laurent polynomials in z over the integers. 

Our goal is to find Laurent polynomials bo{z), bi{z), 62 (-2)? b^i^z), 64(2:), b^i^z) such that 
the expression (1D.8P is zero modulo 16 as Laurent series in z. Lemma [391 combined 
with the explicit forms of M2{z,t) and Ms{z,t) given in Proposition [21 then says that 

8(po{z) + Go{z,h{z),h'{z))) = {Az + 10z'^)c3{z) + Szc^iz) modulo 16, 

8(^Pi{z) + Gi {z, h{z), h'{z))^ = {Az + 10z^)c2{z) + (12 + Az)c3{z) + Szc^iz) + 805(2) 

modulo 16, 

8(p2{z) + G2{z, h{z), b'(z))) = (12 + 4z)c2iz) + (6 + 4z)c3{z) + 804(2) + Sc^iz) 

modulo 16, 

8(p3{z) + G3{z,h{z),h'{z))^ =(6 + 4z)c2(z) + 12c3(z) + 8c4(z) modulo 16, 

8(pa{z) + G4{z,h{z),h'{z))) = 12c2{z) +203(2) modulo 16, 
8(p5{z) + G5{z,h{z),h'{z))) =202(2) modulo 16, 

for suitable Laurent polynomials 02(2), 03(2), 04(2), 05(2). From the last congruence one 
sees that 02(2) is actually zero modulo 4, and then the next-to-last congruence implies 
that the same holds for 03(2). Writing 4a(2) = 02(2), 46(2) = 03(2), c{z) = 04(2), 
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d{z) = c^iz), we see that the above system of congruences simphfies to 
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G'o(^,b(^),b'(^) 

G2{zMz),h'{z) 
G3(^,b(z),b'(z) 
G4(^,b(z),b'(z) 
G5(;2,b(z),b'(z) 



Po{z) + z'^b{z) + zd{z) modulo 2, 

Pi{z) + z'^a{z) + zc{z) + d{z) modulo 2, 

P2{z) + b{z) + c{z) + d{z) modulo 2, 

P3{z) + a{z) + c{z) modulo 2, 

P4,{z) + b{z) modulo 2, 

P5{z) + a{z) modulo 2. 



(D.9) 



This puts us in the situation of Lemma [T21 except that on the right-hand sides of 
the congruences (denoted by ri{z), i = 1,2, ... ,N, in Lemma [T2|) there appear the 
unknown Laurent polynomials a{z) , b{z) , c{z) , d{z) . Still, the idea of the proof of 
Lemma [12] may be applied: the system of congruences (JD.9P can be solved with re- 
spect to the "variables" bo{z), bi{z), 62 (^), b^i^z), 64(2;), b^i^z) by separating odd and even 
parts, and thereby converting the original system (14.91) of congruences into the system 
(I4.10p of linear congruences for the odd and even parts of the variables. We solve 
this last system over the field of rational functions over Z/2Z, where odd and even 
parts of the "auxiliary variables" a{z) , b{z) , c{z) , d{z) "sit" inside the odd and even 
parts of the "constants" ri{z). In the end, if odd and even parts of the variables 
bQ{z),bi{z),b2{z),b3{z),b4{z),b5{z) are put together, then we are able to express these 
variables in the form 



bAz) 



q,iz) + H,{z,a("\z),a('\z), . . . ,S"\z),d('\z)) 



P(z) 



modulo 2, i = 0,1, ... ,5, 



(D.IO) 
where the qi{z)^s are (known) Laurent polynomials in z over the integers, P{z) is a 
(known) polynomial in z over the integers, and the 

Hi{z, a^°\z), a^^\z), ..., d^"\z),d'^''\z))'s 

are linear forms in a'^°\z), a^^\z), . . . , d^°\z), d^^\z) with coefficients that are (known) 
Laurent polynomials in z over the integers. 

The task now is to choose a^°\z) , a^^\z) , . . . , d^°\z) , d^^\z) in such a way that in 
each of the fractions on the right-hand sides of ( ID.10|) the denominator P{z) cancels 
out. 

In order to carry out this task, we decompose P{z) into its prime factors (over Z/2Z), 
say 



P(z) = Yl Pp (z) modulo 2. 



i=i 
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Using a standard inductive procedure jij we find a^°\z) , a^^\z) , . . . , S°\z) , S^\z) (if 
tliere are") sucli tliat 



qi{z) + H,{z,a^"\z),a^-\z), . . . ,S''\z),S-\z)) = {mod Pp (z)) , 

z = 0,l,...,5, j = l,2,...,i, (D.ll) 

(again, over tfie field Z/2Z), and tlien put the particular results for each j together by 
means of the Chinese remainder theorem. We only discuss the generic case here, the 
discussion for other cases being completely analogous. Namely, generically, having to 
solve 6 equations in 8 variables, one will be able to express six of the variables in terms 
of two "free" variables. Let us say, b^"\z) , b^^\z) , c^^^z) , c^'^^z) , S"\z) , S''\z) can be 
expressed in terms of a^°''{z), a^'^\z), 

+ v,iz)a^'\z) 



b^"){z) = si{z) + ui{z)a^°\z 
b^^\z) = S2{z) + U2{z)a^"\z 
c^"\z) = s;iz)+u,iz)a^"\z 
c^^){z) = s^{z) + Ui{z)a^"\z 
S°\z) = S5{z) + U5{z)a^°\z 
S'\z) = seiz) + UG{z)a^°\z 



V2{z)a 
V2,{z)a 



+ v^{z)a^'\z) 



(e)( 
(e)( 



modulo 2, 
modulo 2, 
modulo 2, 
modulo 2, 



+ v^{z)a'''\z), modulo 2, 
+ VQ{z)a'''\z) modulo 2, 



(D.12) 



where the Sj(2;)'s, the Mj(2;)'s, and the fj(z)'s are certain (known) Laurent polynomials 
in z over the integers, and where we are free to choose a^"\z) and a^^\z). If this is 
substituted in (ID. 101) . then on the right-hand sides the denominator P{z) cancels out, 
and bQ{z),bi{z),b2{z),b3{z),bi{z),b^{z) will all be equal to Laurent polynomials in z 
over the integers. 

We are still not finished, though. In the "solution" (ID. 120 the Laurent polynomials 



Mi 



5(0) ( 



z),c 



(0)1 



,S-°\z) must be chosen as odd, while the Laurent polynomials 



a'^^\z)^b'^^\z)^c'^^\z)^S-^\z) must be chosen as even. In order to achieve this, we must 
(again) separate odd and even parts: doing so in (ID. 121) yields the system 



= 


s^\z)+u^^ 


{z)a<^"\z)+vi^\z)a^^\z) 


modulo 2, 


= 


S2 {z) + U2 


iz)a^°\z)+vi''\z)a^'\z), 


modulo 2, 


= 


sf{z)+uf 


{z)a^'^\z)+v,ie)iz)a^^\z), 


modulo 2, 


= 


s^\z)+uf 


iz)a^°\z)+vi"\z)a^^\z), 


modulo 2, 


= 


st\z)+uf 


iz)a('^\z)+vi^\z)a^^\z), 


modulo 2, 


= 


sf{z)+uf 


(^)a(°)(z)+4")(^)a(^)(^) 


modulo 2. 



(D.13) 

This is a system of six linear congruences with two variables, a^°\z) and a^^\z), where 
the first of these should be an odd Laurent polynomial and the second an even one. 



One first solves (JD.III) modulo Pjiz) (instead of P'"'(z)); this means solving a system of linear 
equations over a field. If one has solved (|D.11[) already modulo P,'*(z), for each variable var(z) one 
makes the Ansatz var(z) = varo(z) + vari(z)P''(z), where varo(2;) is the value of var(z) in the solution 
modulo Pj^iz). If this is substituted in (|D.11|) . after cancellation, solving (jD.lip modulo FJ'+^(z) boils 
again down to solving a system of linear equations modulo Pj{z), that is, over a field. 
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It is of the type of the system of congruences fl4.10p . How to solve such a system is 
explained in the paragraph below the proof of Lemma [T21 (One first solves over the 
field of rational functions in z over 'L/2'L, and then cancels denominators, if possible.) 
Moreover, the argument in the proof of Lemma [12] showing that, if fl4.10p has some 
solution in Laurent polynomials, then it also has a solution in which all /• (2;)'s are 



.•(2) 



even Laurent polynomials and all /• (-2)'s are odd Laurent polynomials, also applies to 



'3 

the system (ID. 131) to guarantee that, if one is able to find some solution a^''^(z), a^^\z)^ 
then one can also find one in which a^°\z) is an odd Laurent polynomial and a^^\z) is 
an even Laurent polynomial. 

If one is able to carry through this procedure, then one has found the unknowns 
bi,i{,z), -i = 0, 1, . . . , 5, so that (]D.7|) produces the desired description modulo 16 of the 
solution F{z) to the Riccati-type differential equation (14.11) . Conversely, if one of the 
systems of linear congruences which one has to solve along the way (these are (10.90 . 
(1D.11|) . and (JD.ISP ) has no solution, then one has proved that it is impossible to describe 
the series F{z) modulo 16 in terms of a polynomial in ^{z) with coefficients that are 
Laurent polynomials in z over the integers. 
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